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GLOSSARY  OF  SYMBOLS 


Notation 

Upper-case  is  used  to  denote  both  the  signal  inputs  and  outputs  of 
logic  blocks  (e.g.,  X3(t),  Y2(t),  Z2(t)).  Upper-case  x  is  reserved  for 
the  inputs  to  the  first  level  of  logic  (e.g.,  XI (t),  X2 (t) ... .Xn(t))  while 
upper-case  z  is  reserved  for  the  outputs  from  the  last  level  of  logic 
(e.g.,  Zl(t),  Z2(t) , . . ,Zm(t)) .  Signals  defined  at  points  internal  to  a 
logic  block  are  denoted  by  lower-case  (e.g.,  x(t),  y(t),  z4(t)).  When 
modeling  logic  blocks,  delay  elements  are  always  placed  at  the  output. 

The  lower-case  notation  denoting  the  input  to  the  delay  element  is  always 
identical  to  the  upper-case  notation  denoting  the  output  (e.g.,  if  the 
delay  element  input  is  denoted  by  z3(t)  the  delay  element  output  will 
be  denoted  by  Z3(t)). 

A  random  variable  (r.v.)  is  denoted  by  a  symbol  with  an  underbar 
(e.g.,  _T ,  xi*  uO  •  A  stochastic  process  is  denoted  as  a  time  function 
with  an  underbar  (e.g.,  X2(t),  £(t),  Z(t)).  The  time  variable  is  denoted 
by  t.  A  specific  time  instant  is  denoted  by  t^;  k  =  1,2,...  Time  delays 
are  denoted  by  i.  Time  derivatives  are  denoted  by  an  upper  dot  (e.g., 
y(t),  Xl(t)). 

Probabilistic  quantities  are  denoted  as  follows: 

F  (X)  =  Pr  {X  <  X)  =  Cumulative  Distribution  function  (c.d.f.) 

A 

fx(X)  =  —  F^(X)  =  probability  density  function  (p.d.f.) 

E[X]  =  /  X  fy  (X)dX  =  expected  value  of  the  r.v.  X. 


E[X(t)j  ■  expected  value  of  the  stochastic  process  X(t)  . 

• 

E[X(t)]  =  —  E[X(t)]  =  time  derivative  of  the  expected  value  of  the 
process  X(t)  . 


RXY  =  =  cross  correlation  function  of  the  random 

variables  X(tp  and  Y(t2)  • 

(n)  =  E[X(t)Y(t-n)]  =  cross  correlation  function  when  X(t)  and 
Y(t)  are  mutually  stationary  random  processes. 


X(t)>  = 


tf-t0 


X(t)dt  =  time  average  of  X(t)  over 


the  finite  observation  time  interval  (t  ,t,)  . 

o  f 


X(t)*Y(t)  = 


X(0)Y(t-0)d0  =  convolution  operation  between  the  functions 


X(t)  and  Y (t) . 

The  notation  associated  with  0,1  binary  signals  is  as  follows  (also, 
see  Fig.  1.5): 

u^  =  time  at  which  the  i—  rise  transition  occurs 

d^  =  time  at  which  the  i— —  fall  transition  occurs 

th 


u.  =  u.-d.  ,  =  duration  of  the  i —  gap 
lii-l  &  * 


th 


6.  =  d.-u.  =  duration  of  the  i— —  pulse 
ill  v 


U.  =  u.-u.  .  =  duration  of  the  i-^-  interval  between  successive  rise 
l  li-l 

transitions 

=  d^-d^_j  =  duration  of  i— —  interval  between  two  successive  fall 
transitions . 


X  (t)  =  counting  signal  for  the  rise  transitions  of  X(t) 


X  (t)  =  counting  signal  for  the  fall  transitions  of  X(t) 


-rT’  -fT 


-rG’  -fG 


E.[Z(t)] 


yt> 


Total  delays  for  rise  and  fall  output  transitions  in  the 
longest  path. 

Mean  path  delays  for  rise  and  fall  output  transitions. 

Weighted  mean  path  delays  for  rise  and  fall  output  transitions. 
Assigned  Gaussian  delays  for  rise  and  fall  output  transitions. 
Output  expected  value  calculated  using  i*''1  strategy  for 
delay  characterizations. 

Reference  waveform  for  E[Z(t)],  obtained  by  Monte  Carlo 
simulation. 


The  notation  utilized  in  the  discussion  of  switching  functions 
is  as  follows: 

X'  =  complement  of  switching  variable  X 

XvY=  OR  operation  between  X  and  Y 

XXY  =  X-Y  =  AND  operation  between  X  and  Y 

(X*Y) '  =  NAND  operation  between  X  and  Y 
(XvY) '  =  NOR  operation  between  X  and  Y 

X©Y  =  XOR  (exclusive  OR)  operation  between  X  and  Y 
X0Y  =  EXNOR  (exclusive  NOR)  operation  between  X  and  Y 


V  denotes  end  of  example 


viii 


Schematic  Symbols 


Buffer 


AND  gate 


OR  gate 


Inverter 


NAND  gate 


NOR  gate 


ill 


1 .  INTRODUCTION 


1.1  Motivation 

In  recent  years  the  trend  has  been  towards  the  fabrication  of 
microelectronic  circuits  capable  of  operating  at  increased  speeds. 

By  way  of  example,  the  goal  of  the  VHSIC  (Very  High  Speed  Integrated 
Circuits)  program  is  to  achieve  clock  rates  as  high  as  100  MHZ.  As 
integrated  circuits  become  faster,  the  occurrence  of  unintentional  prop¬ 
agation  delays  within  the  chip  can  seriously  affect  system  performance. 

Because  microelectronic  circuits  are  highly  complex,  they  are  ex¬ 
tremely  difficult  to  analyze.  The  situation  is  further  complicated 
when  the  circuits  are  exposed  to  electromagnetic  interference  (EMI) . 

In  an  effort  to  simplify  the  analysis  problem,  a  probabilistic  approach 
has  been  proposed  [1].^  Using  this  point  of  view,  the  desired  and  inter¬ 
fering  signals  are  treated  as  stochastic  processes  while  the  suscepti¬ 
bility  level  of  the  system  is  considered  as  a  random  variable.  Approxi¬ 
mations  to  probability  density  functions  (p.d.f.'s)  are  obtained  through 
measurements  and  statistical  methods.  These  are  then  used  to  evaluate 
the  probability  of  unacceptable  system  performance  due  to  EMI. 

The  probabilistic  nature  of  integrated  circuits  was  observed  in 
an  experiment  performed  at  Rome  Air  Development  Center  (RADC)  [2],  The 
measured  susceptibility  levels  of  252  "equivalent"  7400  TTL  NAND  gates 
were  found  to  vary  in  a  random  fashion  from  0.85  volts  to  a  level  in 
excess  of  3.0  volts.  Commonly  used  p.d.f's  were  shown  to  approximate 
the  data. 

^Numbers  in  brackets  refer  to  listings  in  the  list  of  references. 


This  dissertation  is  concerned  with  the  increase  in  propaga¬ 
tion  delay  caused  by  EMI,  as  observed  by  Alkalay  and  Weiner  [3] .  In 
their  computer  simulation  sinusoidal  EMI  was  injected  into  the  output 
of  a  7400  TTL  NAND  gate  which  was  loaded  by  a  fan  out  of  10  identical 
gates.  The  gate  under  test  was  driven  by  an  input  which  caused  the. 
output  to  switch  from  the  HIGH  to  LOW  state  in  the  absence  of  inter¬ 
ference.  As  expected,  for  a  large  enough  interfering  signal,  the 
gate  output  would  not  remain  in  the  LOW  state  after  switching.  However, 
for  smaller  interfering  signals,  even  though  the  gate  output  would 
remain  in  the  LOW  state  after  switching,  increases  in  the  propagation 
delay  were  noted.  The  increment  in  the  delay  was  referred  to  as  "excess 
propagation  delay."  A  plot  of  excess  propagation  delay  versus  frequency 
is  shown  in  Fig.  1.1  with  the  peak  amplitude  of  the  interfering  signal 
as  a  parameter.  Since  typical  propagation  delays  in  the  absence  of 
interference  are  in  the  order  of  10  nanoseconds,  some  of  the  observed 
excess  propagation  delays  were  considerably  larger  than  the  typical 
delays  without  interference. 

An  illustration  of  how  unintentional  delays  can  upset  system  per¬ 
formance  is  given  in  Example  1.1. 

Example  1 . 1  Consider  the  simple  combinational  logic  circuit  shown  in 
Figure  1.2.  Assume  a  HIGH  level  is  assigned  the  value  1  while  a  LOW 
level  is  assigned  the  value  0.  Input  B  is  assumed  to  be  LOW  for  all 
time.  The  inputs  at  terminals  A  and  C  are  denoted  by  Xl(t)  and  X2(t), 
respectively.  For  the  waveforms  shown,  the  output  Z(t)  is  identically  0. 

In  the  presence  of  EMI,  assume  the  output  of  gate  01  is  delayed  1.5  ns 
while  the  output  of  gate  02  is  delayed  by  0.75  ns.  The  resulting  waveforms 
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Fig.  1.2  Combinational  logic  circuit  and  signals  for  Example  1.1 


and  the  output  Z(t)  are  plotted  in  Figure  1.3.  Because  of  the  delays, 
a  pulse  is  now  present  in  Z(t).  This  is  referred  to  as  a  hazard  and 
could  cause  difficulty  depending  upon  the  use  of  the  circuit. 
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Fig.  1.3  Waveforms  for  the  circuit  of  Figure  1.2  illustrating 
propagation  delays  due  to  the  presence  of  EMI. 


1.2  Objective 

The  objective  of  this  work  is  to  evaluate  the  performance  of 
combinational  logic  circuits  in  which  the  signals  experience  random 


propagation  delays.  Except  for  delays,  the  logic  Is  assumed  to  perform 
In  an  ideal  manner.  Random  propagation  delays  are  present,  even  in  the 
absence  of  EMI,  due  to  parameter  fluctuations  in  the  fabrication  process 
of  integrated  circuits.  As  pointed  out  previously,  excess  random  pro¬ 
pagation  delays  may  also  occur  due  to  EMI.  In  this  dissertation,  both 
effects  are  combined  together  in  a  single  overall  random  propagation  delay. 
If  the  delays  associated  with  a  combinational  logic  circuit  are 

specified,  it  is  treated  as  deterministic  and  a  known  input  produces  a 

/\ 

known  output.  Denote  the  predicted  output  by  Z(t).  Specific  valifes  for 
the  delays  can  be  obtained  in  a  variety  of  ways.  For  example,  they  may 
be  measured,  they  may  be  chosen  as  the  average  of  the  maximum  and  mini¬ 
mum  delays  specified  by  the  manufacturer,  or,  in  a  worst  case  design, 
chosen  as  the  maximum  delay  quoted  in  the  data  sheets.  In  general, 
the  delays  are  unknown  and  the  actual  output  Z(t)  may  differ  from  the 

A 

predicted  output  Z(t).  One  possible  measure  of  system  performance  is  the 
probability  that  Z(t)  differs  from  Z(t). 

For  analytical  simplicity,  and  without  loss  of  generality,  it  is 
convenient  to  assign  the  value  1  to  a  HIGH  level  and  the  value  0  to 
a  LOW  level.  Therefore,  at  any  time  instant,  each  waveform  assumes 
either  of  the  values  0  or  1.  We  refer  to  such  waveforms  as  0,  1 
binary  signals.  Consider  a  combinational  logic  circuit  in  which  the 

A 

propagation  delays  have  been  specified.  Assume  Z(t^;  =  0.  The  proba- 

A 

bility  that  Z(tp  differs  from  Z(tp  equals  the  probability  that 

A 

Z(tp  =  1.  This  is  also  the  probability  of  error  given  that  Z(tp  =  0. 


We  write 


Pr  {error  |  Z^)  -  0}  -  Pr  {Z(tL)  -  1}  .  (] 

A 

Now  assume  Z(t2)  “1.  It  follows  that 

A 

Pr  (error  [ Z^Ct2>  -  l)  »  Pr  {Z(t,,)  -  0} 

-  1  -  Pr  (Z(t2)  -  1}.  (] 

In  each  case  the  conditional  probability  of  error  can  be  expressed  in 
terms  of  the  probability  that  the  output  equals  unity.  Note  that  the 
expected  value  of  the  output  is  given  by 

E[Z(t)]  =  0  •  Pr  lZ(t)  =0}  +  1  •  Pr  {Z(t)  -  1} 

-  Pr  (Z(t)  ■*  1}.  C 

As  a  result,  the  conditional  probabilities  of  error  can  be  expressed  as 

/\ 

Pr  {error  |  Z(t^  =  0}  =  EtZU^] 

Pr  {error  |  Z(t2)  =  1}  =  1  “  E[Z(t2>],  ( 

Because  of  the  above  considerations,  the  effort  in  this  dissertation 
is  focused  on  evaluating  the  output  expected  value.  Example  1.2  illus¬ 
trates  the  application  of  equations  (1.4). 

Example  1.2  For  a  given  input  and  a  specified  set  of  delays,  assume  the 

A 

output  Z(t)  is  derived  deterministically  to  be  the  waveform  shown  in 

Figure  1.4(a).  Taking  the  random  nature  of  the  propagation  delays  into 


account,  let  the  expected  value  of  the  output, E[Z(t) ] »be  as  depicted 

/\ 

in  Figure  1.4(b).  The  conditional  probability  of  error  given  Z(t)  is  de¬ 
termined  frota  Eq.  (1.4)  and  is  sketched  in  Figure  1.4(c).  In  general. 


Z(t) 


t 


Fig.  1.4  Waveforms  for  Example  1.4. 

(a)  derived  signal  Z(t),  (b)  expected  value  E[Z(t)], 

(c)  conditional  error  Pr  {error |Z(t)}  as  a  function  of  time. 


the  conditional  probability  of  error  will  be  strongly  dependent  on  the  pre 

/s 

dieted  output  Z(t)  which,  of  course,  depends  upon  the  specified  delays 
and  the  given  input.  V 

1.3  Problem  Description 

In  this  section  the  three  basic  considerations  which  enter  into  the 
problem  description  are  discussed.  They  are:  1)  the  signals  involved, 

2)  the  subdivision  of  combinational  logic  circuits  into  logic  blocks, 
and  3)  the  interconnection  of  the  various  logic  blocks. 


As  mentioned  in  Sec.  1.2,  all  of  the  digital  signals  are  assumed  to 
be  0,1  binary  (i.e.,  take  on  the  values  0  and  1  with  zero  transition  Limes 


from  one  state  to  the  other).  A  typical  signal  is  shown  in  Fig.  1^.5. 

The  following  notation  is  used: 

*  time  at  which  the  i^  rise  transition  occurs  . 

(A  rise  transition  is  transition  from  0  to  1.) 

=  time  at  which  the  i^  fall  transition  occurs  . 

(A  fall  transition  is  a  transition  from  1  to  0.) 

=  u^  -  d^_^  ■=  duration  of  the  i*"*1  gap.  (A  gap  corresponds  to 
the  interval  between  transitions  during  which  the  signal  is  0) 
6^  =  d^  -  u^  =  duration  of  the  iC^  pulse  . 

(A  pulse  corresponds  to  the  interval  between  transitions 
during  which  the  signal  is  1  ) . 

U.  =  u.  -  u.  .  *  duration  of  ifc^  interval  between  two  successive 
1  x  x-1 

rise  transitions. 

=  d^  -  d^  ^  «■  duration  of  i^  interval  between  two  successive 
fall  transitions. 

By  convention,  the  first  fall  transition  occurring  after  t  =  0  is 


denoted  by  d^.  u^  is  the  rise  transition  preceding  d^  and  may  or 
may  not  be  positive.  Both  deterministic  and  stochastic  0,1  binary  signals 
are  included  in  this  work.  For  deterministic  signals  the  rise  and  fall 
transition  times  are  assumed  to  be  specified.  On  the  other  hand,  u^  and 
d^  are  viewed  as  random  variables  for  stochastic  processes. 

Each  0,1  binary  signal  x(t)  can  be  decomposed  into  two  counting 
signals  as  illustrated  in  Fig.  1.6.  A  counting  signal  is  a  nondecreasing 
integer  valued  step-like  waveform  in  which  every  step  has  unit  height. 
x+(t)  accounts  for  the  rise  transitions  of  x(t)  while  x_(t)  accounts 
for  the  fall  transitions.  By  definition, 

x ( t)  =  x+(t)  -  x  (t)  .  (1.5) 

Observe  that  a  rise  transition  always  precedes  a  fall  transition,  and 
vice  versa.  Symbolically, 


Ui  '  di  "  Ui+1 


di-l  <  Ui  <  di 


(1.6) 


Only  combinational  logic  circuits  are  considered  in  this  work. 
The  basic  building  blocks,  which  correspond  to  the  smallest  entities 
needed  for  the  analysis,  are  gates,  buffers  and  inverters.  These  are 
modeled  as  ideal  logic  circuits  followed  by  ideal  delay  elements  as 
shown,  in  Fig.  1.7. 


In  practice,  the  delays  experienced  by  rising  transitions  differ  from 
the  delays  x^  of  the  falling  transitions  (See  Figs.  1.7(b)  and  (c)). 
Consequently,  the  ideal  delay  circuit  is  referred  to  as  a  discriminating 
delay  circuit.  The  delays  are  assumed  to  be  random  variables  whose 
randomness  is  caused  by  variation  in  the  fabrication  process  plus  the 
effects  of  EMI.  The  probability  density  function  (p.d.f.'s)  of  x^  and 
x^  are  assumed  to  be  known. 

For  complicated  combinational  circuits  it  is  convenient  to  sub¬ 
divide  the  network  into  larger  logic  blocks.  These  are  also  modeled  as 
ideal  logic  circuits  followed  by  ideal  delay  elements.  It  is  shown 

that  the  same  analytical  techniques  used  with  the  single  logic  gates 
can  be  extended  to  the  more  complicated  logic  blocks. 

Various  interconnections  may  arise  in  combinational  circuits.  The 
network  is  said  to  be  tree-like  when  no  more  than  one  path  exists  from 
each  input  to  every  output  (See  Fig.  1.8(a)).  The  network  is  said  to 
contain  a  reconvergent  fanout  when  two  or  more  branches,  originating 
from  a  single  point,  merge  together  at  some  other  point  of  the  network. 
(See  Fig.  1.8(b)).  By  definition,  networks  containing  reconvergent 
fanouts  are  not  tree-like.  A  feedback  loop  is  said  to  exist  in  a 
network  when  the  output  of  a  logic  block  is  fed  back  to  a  previous 
block  whose  output  is  connected  in  a  path  leading  to  an  input  of  the 
original  block.  (See  Fig.  1.8(c)).  In  this  work  only  tree-like  networks 
and  networks  containing  reconvergent  fanouts  are  considered. 

Given  a  description  of  the  input  signals  and  the  p.d.f's  of 
the  delays,  the  problem  is  to  evaluate  the  expected  value  at  the 


Fig.  1.8  Examples  of  (a)  tree-like  network,  (b)  network  with 
reconvergent  fanout,  (c)  network  with  feedback. 

network  output.  As  pointed  out  previously,  this  can  be  used  as  a 
measure  of  the  output  error.  Since  the  output  may  be  used  as  the  input 
to  another  combinational  network,  it  is  desirable  that  the  analytical 
technique  developed  be  capable  of  generating  the  input  information 
required  for  analysis  of  the  second  circuit. 

1.4  Problem  Classification 

A  variety  of  problems  may  arise  depending  upon  the  signals,  logic 
blocks,  and  networks.  For  the  purpose  of  organization  and  clarity,  a 
classification  scheme  is  proposed  in  this  section. 

As  mentioned  previously,  all  signals  are  0,1  binary.  However, 
during  a  specified  observation  time  they  may  be  either 

1)  time- invariant  (i.e.,  no  transitions  occur) 
or 

2)  time  variant  (i.e.,  transitions  do  occur). 

In  addition,  signals  may  be  either 

a)  deterministic  (i.o.,  all  transition  times  are  specified) 
or 

b)  stochastic  (i.e.  transition  times  are  considered  as  random 


variables  with  known  probability  density  functions). 


In  this  work  the  signals  at  the  network  inputs,  referred  to  as 
primary  inputs,  are  assumed  to  be  deterministic.  However,  because 
deterministic  signals  become  random  after  being  processed  by  a  logic 
block  containing  random  delay,  the  analysis  of  logic  blocks  is  performed 
assuming  stochastic  inputs. 

In  this  work  only  zero  memory  logic  blocks,  such  as  gates, 
buffers,  and  inverters,  are  considered.  They  may  be  further  classi¬ 
fied  according  to  whether  the  delays  associated  with  rise  and  fall  tran¬ 
sitions  are 

1)  equal  (i.e.,  signals  delayed  without  distortion) 
or 

2)  unequal  (i.e.,  signals  delayed  with  distortion)  . 
and  whether  they  are 

a)  zero  (i.e.  no  delay) 

b)  specified  (i.e.  deterministic) 

c)  statistically  independent  random  variables 

d)  statistically  dependent  for  a  given  logic  block  but  statisti¬ 
cally  independent  from  one  logic  block  to  another 

e)  statistically  dependent  . 

With  regard  to  the  networks,  only  combinational  networks  are  con¬ 
sidered  (i.e.,  no  feedback  loops  exist).  However,  the  networks  may  either 
have 

or  1)  no  reconvergence  (i.e.,  be  tree-like) 

2)  one  or  more  reconvergences. 

In  the  ensuing  discussion  the  proposed  classification  procedure 
will  be  used  to  characterize  the  various  problems  analyzed. 


1.5  Literature  Review 

The  earliest  attempts  at  probabilistic  analysis  of  digital  logic 
circuits  dealt  with  evaluation  of  the  probability  that  a  given  switching 
expression  (i.e.,  binary  Boolean  function)  takes  on  the  value  1.  Fratta 
and  Montanari  [4]  introduced  algorithms  for  evaluating  this  probability 
using  an  orthogonal  form  of  the  switching  expression.  Parker  and 
McCluskey  [5,6]  formalized  the  theory  by  presenting  several  theorems 
which  show  the  relationship  between  Boolean  operations  and  algebraic 
operations  upon  probabilities.  Kumar  and  Breuer  [7]  derived  additional 
results  and  showed  how  their  approach  is  related  to  the  Walsh  coeffi¬ 
cients  and  Reed-Muller  canonic  form  of  the  switching  expression. 
Murchland  [8]  considered  multi-state  systems  and  evaluated  the  probabil¬ 
ity  of  being  in  a  state  as  well  as  the  average  transition  rates  between 
states.  In  particular,  he  showed  how  to  use  the  Karnaugh  map  to  deter¬ 
mine  the  transition  rates.  Schneeweiss  [9]  presented  a  conceptually 
simple  method  which  lends  itself  to  efficient  hand  calculations  for 
some  problems.  Debany  [10]  developed  a  computationally  efficient  com¬ 
puter  program  for  evaluating  the  probability.  The  above  investigators 
were  motivated  by  a  variety  of  applications.  Fratta  and  Montanari 
were  interested  in  computing  the  terminal  reliability  in  a  communica¬ 
tion  network.  Parker  and  McClusky  applied  their  results  to  the  analysis 
of  faulty  logic  circuits.  Kumar  and  Breuer  were  interested  in  the 
random  testing  of  digital  circuits.  Murchland  was  motivated  by  the 
reliability  analysis  of  systems  with  numerous  components.  Schneeweiss 


applied  his  results  to  the  evaluation  of  fault  trees  used  in  reliability 
analysis  while  Debany  was  also  interested  in  the  random  testing  of 
digital  networks.  All  of  the  work  cited  above  is  applicable  to  evalu¬ 
ation  of  the  output  expected  value  of  a  digital  circuit  without  delay. 

A  first  attempt  at  introducing  delay  was  made  by  Coraluppi  [11]. 

His  approach  is  a  graphical  one  in  which  deterministic  pulses  are  mapped 
onto  the  complex  plane  by  associating  the  time  instants  at  which  the 
pulses  rise  with  the  imaginary  axis  and  the  time  instants  at  which 
the  pulses  fall  with  the  real  axis.  A  pulse,  therefore,  is  represented 
by  a  vector  in  the  complex  plane.  A  discriminating  delay,  as  described 
in  Fig.  1.7,  is  accounted  for  by  adding  the  vector  T ^  to  the 

vector  representing  the  pulse.  Coraluppi  gives  an  interesting  graphical 
interpretation  to  the  effects  of  delays  associated  with  basic  gates 
in  the  output  of  a  digital  system. 

In  recent  years  significant  effort  has  been  devoted  to  studying 
random  delays  in  digital  circuits  by  means  of  digital  simulation. 
Magnhagen  and  Flisberg  [12]  characterize  the  delays  by  assigning  Gaussian 
probability  density  functions  (p.d.f.'s)  to  each  delay.  The  means  and 
variances  are  determined  from  the  maximum  and  minimum  delays  specified 
by  the  manufacturer.  Applications  of  their  simulation  program,  called 
DIGSIM,  to  digital  design  verification  of  race  conditions  are  given 
in  [13].  It  is  important  to  recognize  that  DIGSIM  does  not  perform 
Monte  Carlo  simulations.  Instead  probabilistic  results  are  obtained  by 
evaluating  the  mean  and  variance  of  the  delay  at  the  output  and  cal¬ 
culating  desired  probabilities  from  the  resulting  Gaussian  p.d.f. 


Another  simulation  program  of  interest  is  LDPS  (LSI  Delay  Path 
Specification)  which  was  recently  developed  by  Al-Hussein  and 
Dutton  [14].  LDPS  computes  the  propagation  delays  of  both  rise 
and  fall  transitions  for  all  paths  in  an  integrated  circuit.  It 
takes  into  account  various  factors  which  affect  the  delays  such  as 
loading,  fabrication,  and  temperature.  Empirical  equations  involv¬ 
ing  these  factors  were  obtained  using  SPICE  simulations.  Other  simu¬ 
lation  programs  are  discussed  by  Grundmann  [15]  and  Bass  and  Grundmann 
[16]. 

Major  contributions  to  the  probabilistic  analysis  of  digital  cir¬ 
cuits  containing  random  delays  have  been  made  by  Grundmann  [15],  Bass 
and  Grundmann  [16],  and  Grundmann  and  Bass  [17].  They  discuss  equal 
and  unequal  random  propagation  delays  for  rise  and  fall  transitions 
in  both  combinational  and  sequential  digital  circuits.  They  have  also 
written  a  computer  program,  called  SLAP,  which  implements  the  theory. 
SLAP  evaluates  as  a  function  of  time  the  probability  that  the  output 
of  a  digital  circuit  takes  on  the  value  1  without  the  need  to  perform 
a  Monte  Carlo  simulation. 

The  work  in  this  dissertation  expands  on  the  excellent  contribu¬ 
tions  made  in  [15]  and  [16].  First,  many  of  the  results  in  [15]  and 
[16]  are  derived  from  a  different  point  of  view.  The  new  derivations 
are  more  general  and,  it  is  believed,  more  easily  understood.  In 
addition,  some  of  the  approximations  used  by  Grundmann  and  Bass  when 
analyzing  networks  containing  reconvergent  fanouts  have  been  removed 
by  using  higher  order  statistics.  Finally,  strategies  are  investigated 
for  significantly  reducing  the  amount  of  computation  by  subdividing  a 
complex  circuit  into  large  blocks. 


1.6  Dissertation  Outline 

Chapter  2  is  devoted  to  the  analysis  of  tree-like  networks  having 
statistically  independent  inputs.  The  case  of  pure  delay  (i.e.,Tr  =  x^) 
is  discussed  first.  It  is  shown  that  E[ZJ ,  the  expected  value  at  the 
output  of  a  logic  block,  is  obtained  by  convolving  E[z_]  ,  the  expected 
value  of  the  output  from  the  ideal  logic  circuit,  with  f^x),  the  p.d.f. 
of  the  pure  delay  element.  The  probability  expression  theorem  is  then 
introduced  as  a  simple  way  for  evaluating  E[_z]  in  terms  of  the  expected 
values  of  the  inputs. 

The  case  of  discriminating  delay  (i.e.,  x^  ^  i^)  is  discussed  next. 
After  decomposing  the  0,1  binary  signals  into  the  differences  between 
their  counting  signals,  an  expression  for  E[Zj  is  derived  which  involves 

•  -f-  •  _ 

convolutions  of  the  time  derivatives  E[£  ]  and  E[jz  ]  with  ^Ir(Tr)  and 
fTf  (Tf)»  respectively.  It  is  shown  how  to  obtain  arithmetic  expressions 
for  z.  and  £  in  terms  of  time  derivatives  of  the  input  signals  by 
employing  the  arithmetic  expression  theorem. 

Finally,  a  systematic  method  is  developed  for  analyzing  tree-like 
networks  consisting  of  several  logic  blocks. 

Chapter  3  discusses  the  analysis  of  combinational  circuits  con¬ 
taining  reconvergent  fanouts.  The  analysis  is  considerably  more  com¬ 
plicated  because  output  higher-order  moments  are  now  required  even 
though  the  primary  inputs  are  statistically  independent.  As  in  Ch.  2, 
the  pure  delay  and  discriminating  delay  cases  are  discussed  separately. 
For  the  pure  delay  case,  the  arithmetic  expression  theorem  simplifies 
determination  of  the  higher-order  moments  at  the  output  of  the  ideal 


logic  circuit.  For  the  case  of  discriminating  delay,  higher-order 

*+  •  — 

joint  moments  involving  the  time  derivatives,  z  and  z_  ,  are  needed. 

The  same  procedure  used  in  Chapter  2  is  employed  for  obtaining  arithme- 
tic  expressions  for  z_  and  z^  .  Having  derived  the  basic  approach,  a 
method  for  analyzing  combinational  circuits  with  reconvergent  fanouts 
is  developed.  Finally,  an  approximation  utilized  without  discussion  by 
a  previous  investigator  [15]  is  analyzed  and  its  limitations  pointed  out. 

To  simplify  the  analysis  of  very  large  circuits,  it  is  proposed 
in  Chapter  4  that  the  circuits  be  subdivided  into  large  logic  blocks 
containing  several  levels  of  logic.  Guidelines  for  appropriate  sub¬ 
division  of  the  network  are  given.  Possible  approaches  for  placing 
the  delay  elements  in  the  model  of  the  large  logic  block  are  considered. 
The  chosen  approach,  referred  to  as  the  simple  output  delay  model,  is 
justified  in  terms  of  simplicity  and  application  of  the  model  to  situa¬ 
tions  involving  EMI.  Various  strategies  for  characterizing  the  output 
random  delays  are  presented  and  applied  to  a  specific  example.  The 
results  are  compared  with  each  other  as  well  as  with  a  computer  genera¬ 
ted  Monte  Carlo  simulation. 

Finally,  a  summary  and  suggestions  for  future  work  are  given  in 


2.  ANALYSIS  OF  TREE-LIKE  NETWORKS 


2. 1  Introduction 

As  mentioned  in  Chapter  1,  only  combinational  logic  circuits  are 
considered  in  this  work.  When  complicated  combinational  circuits  are 
encountered,  it  will  be  convenient  to  subdivide  the  network  into  basic 
logic  blocks.  The  notation  to  be  used  for  the  various  signals  is  illus¬ 
trated  in  Fig.  2.1.  Following  the  practice  of  other  investigators 
[14-17],  basic  logic  blocks  are  modeled  as  ideal  logic  elements  cascaded 
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Fig.  2.1  A  combinational  circuit  which  has  been  subdivided  into 


with  ideal  delay  circuits.  The  ideal  logic  circuits  are  characterized 
by  switching  functions  and  the  logic  is  performed  without  delay.  The 
ideal  delay  circuit  delays  rise  transitions  by  Tr  and  fall  translations 
by  t^.  In  general,  Tr  ^  T£  and  the  ideal  delay  circuit  is  referred 
to  as  a  discriminating  delay.  Observe  that  upper-case  is  used  to  denote 
both  the  signal  inputs  and  outputs  of  logic  blocks.  In  addition,  upper 
case  x  is  reserved  for  the  inputs  to  the  first  level  of  logic  while 
upper-case  z  is  reserved  for  the  outputs  from  the  last  level  of  logic. 
Signals  defined  at  points  internal  to  a  logic  block  are  denoted  by  lower¬ 
case.  With  respect  to  the  discriminating  delay  elements,  the  lower  case 
notation  denoting  their  inputs  is  always  identical  to  the  upper-case  nota¬ 
tion  denoting  their  outputs. 

A  typical  input-output  pair  for  a  discriminating  delay  circuit  is 

shown  in  Fig.  2.2.  Note  that  the  rise  and  fall  transition  instants  of 
k  z(t) 


|Z(t) 

i ! 


Fig.  2.2  A  typical  input-output  signal  pair  for  a  discriminating 
delay  circuit  , 

the  output  are  related  to  those  of  the  input  by 


u„,=u.+t  ;  i  =  1,2,. 

Zi  zi  r 


(2.1) 


dZi  =  dzi  +  lf  5  1  "  1,2,1 
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A  relatively  simple  case  to  analyze  occurs  when  Tr  “  Tf  Then  the  delay 
circuit  output  is  a  delayed,  hut  undistorted,  version  of  the  delay  circuit 
input.  This  type  of  delay  is  referred  to  as  a  pure  delay  and  is  dis¬ 
cussed  first. 

2.2  Pure  Delay  (t^  ■  T^) 

A  typical  logic  block  when  lr=  Tf  M  is  illustrated  in  Fig.  2.3. 

For  ease  of  discussion,  a  single  output  2(t)  is  considered.  The  input 


Fig.  2.3.  A  typical  logic  block  with  Tr  =  Tf  =  T- 

signals  Xi (t) .... ,Xn(t)  are  assumed  to  be  stochastic  processes.  Note 
that  a  stochastic  process  is  denoted  as  a  time  function  with  an  under¬ 
bar.  Similarly,  a  random  variable  (r.v.)  is  denoted  as  a  symbol  with  an 
underbar.  As  discussed  in  section  1.2,  system  performance  can  be  evalua¬ 
ted  in  terms  of  E[Z^(t)]. 

Because  J_r  =  J.£  =  J_»  the  output  of  the  delay  circuit  is  a  delayed 

replica  of  its  input.  Specifically, 


It  follows  that 


E[Z(t)]  -  E[z(t  -  T) ] 


(2.3) 


The  randomness  associated  with  .z(t-jr)  is  due  to  the  randomness  of 
1)  the  input  signals  XI  (t) «...  ,Xn(t)  and  2)  the  delay  t_.  Assume  the 
probability  density  function  (p.d.f.)  of  jt  is  known  to  be  f^(i). 

A  convenient  way  to  evaluate  the  expectation  in  (2.3)  is  to  employ  the 
conditional"  expectation  theorem  [18,  p.208).  Let  Ea  [Yj  denote  the 
expected  value  of  Y_  with  respect  to  the  r.v.  ji.  Ec.(2.3)  can  then  be 
written  as 


EXl(t)...,*l(t),T  EX1(C) . Xn(t),X,-<t"-  1 


EXl(t).....Xn(t)[i(t-  ‘  T)lf1<I>dT 


e  W) . Xn(t)  -T))fT<T)  dT' 


(2.4) 


The  latter  result  is  recognized  as  the  convolution  of  E  „  ,  . [z(t) 

2LL*  •  •  •  >- v^i  y  t ) 

with  f  (t).  If  *  denotes  convolution,  Eq.(2.4)  becomes 


EXl(t),...,Xn(t),  T[-(t)]  "  EXl(t) . Xn(t)[-(t)]*fT(t)-  (2>5) 


'  *J 


For  simplicity,  subscripts  on  the  expectation  operator  will  be  dropped 
when  the  meaning  of  the  expectations  is  clearly  understood.  In  its 
simplified  form,  Eq.(2.5)  is 

E[Z(t)]  =  E[z(t)]  *  fx(t)  .  (2.6) 

Example  2 ■ 1  A  simple  graphical  example  is  shown  in  Fig.  2.4.  Observe 
that  the  random  delay  alters  system  performance  as  reflected  by  the  fact 
that  E[Z^(t)]  differs  from  E[z^t)]. 


(c) 

Fig.  2.4  Graphical  example,  (a)  Expected  value  of  input  to  delay 

circuit,  (b)  p.d.f.  of  random  delay,  (c)  Expected  value  of 
output  from  delay  circuit.  V 

A  major  step  in  evaluating  E[Z^t)]  according  to  Eq.(2.6)  is  the 
determination  of  E [_z^(t )  ] .  With  reference  to  Fig.  2.3.  the  ideal  logic 
circuit  is  characterized  by  its  switching  expression  [ 19)— [21] 


z/t)  =  g(Xl(t) , . . . ,Xn(t) ) . 


(2.7) 


Three  basic  operations  of  switching  algebra  are: 


NOT  (Inversion,  complement) 


z  =  X'  =  1 


AND 


z  =  XI  a  X2  =  XI  •  X2  =  < 


OR 


z  =  XI  v  X2 


0,  X  =  1 
1,  X  =  0 

1,  Xl=l  and  X2=l 
0,  otherwise 

0,  X1=0  and  X2=0 
1 ,  otherwise 


(2.8) 


(2.9) 


(2.10) 


It  is  well  known  that  the  NOT  and  AND  operations  form  a  complete  set 
as  do  the  NOT  and  OR  operations.  For  convenience,  all  three  operations 
are  used  in  this  work. 

An  implicant  is  defined  to  be  a  collection  of  switching  variables  or 
their  complements  connected  by  the  AND  operation  (e.g.,  Xl  A  X2'  A  X3 
is  an  implicant) .  Any  switching  function  can  be  expressed  as  a  set 
of  implicants  connected  by  the  OR  operation.  This  form  is  referred 
to  as  the  "sum  of  products"  form. 


Example  2,2  By  definition,  z  =  g(Xl,  X2,  X3,  X4)  =  (XI '  *  X2  •  X4)  •  X3 ' 

is  not  in  a  sum  of  products  form.  However,  using  the  distributive  law  • 

[19],  the  equivalent  sum  of  products  form  is 

z  =  XI '  •  X3’  v  X2  •  X4  •  X3'  V  (2.11) 

Two  implicants  1^  and  are  said  to  be  nonoverlapping  provided  ■  0 

for  all  possible  values  chosen  for  the  switching  variables. 
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When  the  switching  variables  are  random,  it  is  neceusary  to 
evaluate  E[jz].  Specifically,  it  is  required  to  express  IJ[z]  in  terms 

of  the  input  expectations  E[X1] . E[Xn] ■  This  expression  is  referred 

to  as  the  probability  expression.  Other  investigators  [5]-[7],  [9], 
[15],  [16]  have  demonstrated  how  this  can  be  accomplished  for  specific 
situations.  A  general  procedure  is  derived  below. 

Given  the  switching  expression 


£  =  gQ(Xl,...,Xn), 


(2.12) 


the  first  step  in  evaluating  E[_z]  is  to  express  Eq.(2.12)  in  a  sum  of 
products  form.  The  second  step  is  to  modify  the  sum  of  products  form 
into  a  sum  of  nonoverlapping  implicants.  This  procedure  is  referred 
to  as  a  Boolean  orthogonalization  [4]  and  is  based  on  the  relation 


X  vY  =  Xv  X'*  Y 


(2.13) 


Since  X  •  (X' •  Y)  =  0,  the  right  side  of  £q.(2.13)  consists  of  nonover¬ 
lapping  implicants  whereas  the  left  side  does  not.  The  Boolean  orthogon¬ 
alization  procedure  is  carried  out  as  follows: 

1)  Given  z  =  gQ(Xl , . . . ,Xn)  in  a  sum  of  products  form,  rearrange  the 
product  terms  such  that  successive  terms  have  an  equal  or  larger  number 
of  variables.  If  several  terms  have  the  same  number  of  variables,  arrange 
them  such  that  those  having  the  smallest  number  of  complemented  variables 
come  first.  Represent  the  rearrangement  by 


8o^X1 . Xn^  110  V  I20  V  '*•  V  Xk0,0 


(2.14) 
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where  kO  equals  the  number  of  implicants  (product  terms)  in  Eq.(2.14). 

2)  Define 


Z1  “  X10 


(2.15) 


3)  Construct 


gL (X1 .... >Xn)  =  I|Q  •  gQ (XI .... ,Xn) 


I'  •  I  V  T  '  *1  V 

10  10  10  20 


I10  ’  Ik0,0  *  (2-16^ 


Use  DeMorgan's  law  (i.e.  (X  •  Y)  '  =  X'  VY')  and  the  distributive 
law  (i.e.  (X  v  Y)  •  Z  =  X  •  Z  v  Y  •  Z)  to  expand  the  expression  into  a 
sum  of  products  form.  Where  possible  simplify  the  expression  using  the 
idempotency  law  (i.e.,  X  v  X  =  X  and  X  •  X  =  X)  and  the  null  law  (i.e, 

X  •  X’  =  0) .  Finally,  rearrange  the  terms  according  to  step  1 .  Represent 
the  rearrangement  by 


g  ^  (X 1 , . . , ,  Xn)  1  ^  ^  v  I2 1  v  •  •  •  Y  Ij^j  ^ 


(2.17) 


where  kl  equals  the  number  of  implicants  in  Eq.(2.17). 


4)  Define 


z2  xio  v  In* 


(2.18) 


5)  Construct 

g2(Xl , . . . ,Xn)  =  I[L  •  g1(XL,...,Xn) 


=  1  ’  •  1  vT  *T  v  . . .  v  I '  •  I  .  (2.191 

LU  11  11  21  11  kl  ,1  J) 
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Simplify  and  rearrange  (XI , . . . ,Xn)  as  in  step  3  Represent  the  rearrange¬ 
ment  by 


g2(Xl . Xn)  *  I  2  v  I  2  v 


v  I 

k2,2 


(2.20) 


where  k2  equals  the  number  of  implicants  in  Eq.(2.20) 
6)  Define 

•Z3  =  I10  V  I11  V  I12‘ 


(2.21) 


7)  The  procedure  is  repeated  £  times  until  (Xl,...,Xn)  =  0. 
The  desired  sum  of  nonoverlapping  implicants  is  given  by 


z  =  g  (xl, . . . ,Xn) 
as  proven  below. 


z  =  I  v  T  v 
2  10  11 


v  I 


1(2-1)’ 


(2.22) 


Proof 

The  proof  consists  of  showing  that  a)Iln,  I,,,...  I...  are 

I U  11  1  ( it- 1 )  1 

nonoverlapping  implicants  and  b)z  =  z 


(a)  In  general, 


H  <X1 . x"!  ■  IUVI21V  •••  v  w 


a  (Xl>  .  .  ,Xn) 

1 1  (i-l)  A  Il(i-2)  A  •'*  A  I10  A  8o(X1, • • • >Xn> i 


i=l,2, 


(2.23) 


It  follows  that 


g.(xi,...,xn)  a  ilk  =  0,  k  =  0 . (i-1). 

Since  I  is  an  implicant  of  (XI,..., Xn), 

xli  a  Iik  =  0  ;  k  =  0 .  (i-1)  • 

Because  Bq.  (2.25)  holds  for  all  i,  all  the  implicants  of  z^ 
are  nonoverlapping. 

b)  Mathematical  induction  is  used  to  prove  that  z  ■  z 
reference  to  equations  (2 . 14) —  (2 . 16) ,  note  that 

z  =  8Q(X1» . . • »Xn)  =  I1Q  v  (I20  v  ...  v  Itooj0) 


110  V  I10  '  (I20  v  *  *  *  v  Iko,0) 


110  v  g1(Xl, . . .  ,Xn)  ■  zl  v  g^xl, . . .  ,Xn)  . 


Assume  for  j  an  integer  between  2  and  £.  that 


z  =  Zj  v  gj  (XI,..., Xn). 


Observe  that 


It  follows  that  Eq.(2.28)  holds  for  all  j,  j»l . 1.  In  particular, 

z  *  Zj  v  (XI,.. .Xn) .  (2 

Since  g^(Xl , . . . ,Xn)  =  0,  it  is  concluded  that 

2  =  z£  (2 

Example  2.3  Given  the  expression  z  *  g(Xl,  X2,  X3,  X4)  -  X1*X2'  v  X3*X4 
write  the  expression  in  a  sum  of  nonoverlapping  implicants  form. 

Step  1.  The  expression  z  is  already  given  in  a  "sum  of  products"  form 
and  only  rearrangement  is  required. 

z  =  g0(Xl,  X2 ,  X3,  X4)  =  X3  •  X4  v  XI  •  X2 '  . 

Step  2.  zx  =  X3  •  X4  . 

Step  3.  gl(XI,  X2 ,  X3,  X4)  *  (X3  •  X4) ’ •  (X3  •  X4)  v  (X3  •  X4) ' (XI  •  X2 

*  (X3  •  X4)’.(X1  •  X2 ' ) 

=  (X3  '  v  X4  ’  )*(X1  •  X2 ' ) 

=  X3'*  XI  •  X2'  v  X4'  •  XI  •  X2'  . 

Step  4.  z2  =  X3  •  X4  v  X3'*  XI  •  X2' 

Step  5.  g2(Xl,  X2,  X3,  X4)  *  (X3'  •  XI  •  X2')'.(X3’*  XI  •  X2’) 

v  (X3 ’  •  XI  •  X2’)’  * (X4 ’  •  XI  •  X2') 

*  (X3 ' •  XI  •  X2')'  *(X4'*  XI  •  X2') 

-  (X3  v  XI'  VX2)  •  (X4'«  XI  •  X2') 


Step  6.  z3  -  X3  •  X4  v  X3'*  XI  •  X2' vx3  •  X4'*X1  •  X2'  . 

Step  7.  g  (XI,  X2,  X3,  X4)  «  (X3  •  X4’*  XI  •  X2’)*  (X3  •  X4’*  XI  •  X2')  - 


Thus  z  =  z,  -  X3  •  X4  v  X3'*  XI  •  X2'  v  X3  •  X4’*  XI  •  X2'  . 


Note  that  the  implicants  are  nonoverlapping.  V 

Having  demonstrated  how  to  express  the  switching  expression  as 
a  sum  of  nonoverlapping  implicants,  the  case  in  which  the  switching 
variables  are  random  is  now  considered.  Given 


—  •  •  • » Xn)  =  JI  ^ Q  v  _I ^ l  v  ...  v  I j ^ £_ j ^  ,  (2.3 

it  is  necessary  to  evaluate  E[z].  Since  the  implicants  in  Eq.(2.31) 
are  nonoverlapping.,  for  any  choice  of  values ^giver  cne  switching  variables 
only  one  implicant  can  assume  the  value  1.  Therefore, 


E[,l  -  Pr(z  -  1)  -  Pr[I1Qv  ...  vl  -  I] 


.Pr[il0-ll  +  ...  >1 


'  E[I10]  +  ••• 

The  expected  value  of  _z  is  seen  to  equal  the  arithmetic  sum  of  the 
expected  values  of  the  nonoverlapping  implicants. 


Each  implicant  of  z_„  is  of  the  form 


-ik =  -ik  *  -2k’  •••■  4k  ;  k "  0”-*,a  1}- 


where 


Aj ^  £  (^  ,  X  ' »  1)  ;  j  ■  1 , . . .  ,n. 


(2.34) 


Note  that 


EHik>  ’  PtHik  -1' 


»r(Alk  -  1,  A2k  -  i . A,k  -  n 


(2.35) 


Assuming  the  switching  variables  to  be  statistically  independent,  it 
follows  that 


=  Prt-lk  =  ^  Pr*-A2k  “  ^  •••  Pr  = 


*  E[A,k)  E[A2kl  ...  E[4k]. 


Substitution  of  £*q.(2.36)  into  Eq.(2.32)  yields 


E[z]  =  E[A1Q]  E[A2Q]  ...  ElA^] 


+  ...  +  E  tA2(£_1)  1  E[4i(£-1)] 


(2.36) 


(2.37) 


The  previous  development  suggests  the  following  theorem  for  evaluation 
of  the  probability  expression,  E[zJ . 


The  Probability  Expression  Theorem 

Let  z_  =  gp (XI .... ,Xn)  be  a  switching  expression  of  the  n  statis¬ 
tically  independent  switching  variables  2LL»  •  *  •  »2i£.*  If  is  expressed 
as  a  sum  of  nonoverlapping  implicants,  then  E[z]  is  readily  obtained 
by  interchanging  each 

1)  OR  operation,  v  ,  by  the  arithmetic  operation  of  addition, 

2)  AND  operation,  ♦,  by  the  arithmetic  operation  of  multiplication, 

3)  switching  variable,  Xk,  by  its  expected  value. 


A*  ,**  .  ”  »  *  .  *_»*_•  *.«  *. 


“.  %■  *.  \  % 
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Straightforward  application  of  the  probability  expression  theorem 
yields  an  expression  for  E[zJ  in  terms  of  E[Xk]  and/or  E[Xk* ] 

where  k  ■  1,2 . n.  Expectations  involving  the  complement  of  a  - 

switching  variable  can  be  eliminated  by  noting  that 

E[Xk'  ]  »  Pr[Xk'  -  1]  -  Pr[Xk  =  0] 

=  1  -  PrtjOc  =  1]  =  1  -  E[Xk] .  (2.38) 

Example  2.4  Consider  an  ideal  logic  circuit  with  statistically  inde-  * 

pendent  inputs  XI,  j(2,  X3,  X4.  Given  that  the  switching  expression 
for  the  output  is 

z  -  XI'  v  X2  •  X3  v  X4,  (2.39) 

evaluate  the  probability  expression,  E[zJ . 

Since  the  given  switching  expression  is  not  a  sum  of  nonoverlapping 
implicants,  the  first  step  is  to  perform  the  Boolean  orthogonal iza- 
tion.  Following  the  procedure  outlined  previously, 

z_  =  X4  v  )0_ '  v  X2  •  X3 

z~  X4,  gL  =  X4'  •  z_  =  X4 '  *  XI/  v  X4'  •  X2  •  5(3 

£?=  X4  v  X4 '  *  XI  ’ ,  Zl2  =  QM '  •  XI ' ) '  •  (X4 '  •  X2  •  X3) 

=  ( X4  v  5U_ )  •  (X4 *  •  X2  •  X3)  =  }Q  •  X2  •  X3  •  X4' 

£3  =  )C4  V  )a '  •  X4'  v  30_  •  X2  •  X3  •  X4' ,  ^  =  0  '-.v 

"  m  12  °  ~  W  —  (2-40)  ’ .  / 
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Direct  application  of  the  probability  expression  theorem  to 
Eq.(2.40)  results  in 


E[z]  -  E[X4]  +  E[X1_' ]  E  [X4']  +  E[Xl]  E  [X2]  E[X3]  E  [X4‘] 

-  1  -  E[X1]  +  E[3Q]  E[X4]  +  E[5U]  E[X2]  E[X3] 

-  E[X1]  E[X2]  E[X3]  E[X4] .  V  (2.41) 

In  summary,  provided  the  inputs  to  a  logic  block  are  statis¬ 
tically  independent  and  the  signals  experience  a  pure  delay  (i.e., 
x^  ~  J.f  ■  T.) »  the  expected  value  of  the  output  from  the  logic  block 
in  Fig.  2.3,  E[Z(t)],  is  readily  obtained.  First,  E[jz(t)],  the 
output  from  the  ideal  logic  circuit  is  evaluated  by  expressing  z^(t) 
as  a  sum  of  nonoverlapping  imp licants  and  then  using  the  probability 
expression  theorem.  Finally,  E[Z(t)]  is  obtained  from  E[z(t)]  by 
convolving  E[z/t)]  with  the  p.d.f.  of  T,  f  (t) ,  as  derived  in  Eq.(2.6) 

2.3  Discriminating  Delay  (x^  f  x^) 

As  pointed  out  in  Chapter  1,  the  delays  x^  experienced  by  the  rise 
transitions  typically  do  not  equal  the  delays  of  the  fall  transi¬ 
tions.  This  situation  can  be  analyzed  by  decomposing  each  0,1  binary 
signal  into  two  counting  signals,  as  illustrated  in  Fig. 1.6.  x+(t) 

is  referred  to  as  the  rise  counting  signal  while  x  (t)  is  called  the 
fall  counting  signal. 

A  typical  logic  block  with  x  ^  x^  is  illustrated  in  Fig.  2.5. 
Again,  for  ease  of  discussion,  a  single  output  Z^t)  is  considered. 

The  input  to  the  discriminating  delay  element  is 


Circuit 


Delay  element 


n(t) 


Fig.  2.5  A  typical  logic  block  with  Tr  ^  Tf. 
z(t)  =  z+(t)  -  z  (t)  . 


(2.42) 


By  definition,  the  output  of  the  discriminating  delay  element  is 


Z(t)  =  _z+  (t-xr)  -  £  (t-  !f) . 


(2.43) 


Because  the  delay  element  only  delays  the  rising  and  falling  transitions 
of  z(t),  it  follows  that 


Z^(t)  =  _z+  (t-Tr) 
I  (t)  =  £ (t-xf)  . 


(2.44) 


Therefore,  Eq.(2.43)  may  also  be  expressed  as 
Z(t)  =  Z+(t)  -  Z“(t). 


(2.45) 


Of  interest  is  E[Z_(t)].  Unfortunately,  a  straightforward  method 
for  directly  obtaining  z+(t)  and  z  (t)  given  the  counting  signals  foi 


¥  \ 


XI (t) .... ,Xn(t)  is  not  known.  Difficulties  arise  because  the 
counting  signals  are  not  binary.  The  problem  is  solved  by  con¬ 
sidering  time  derivatives  of  the  counting  signal  expectations.  Later 
in  this  section  a  method  is  derived  for  obtaining  E [£  (t)]  and 

E[jz-(t)]  where  the  dot  denotes  differentiation  with  respect  to  time. 

• 

For  now,  it  is  shown  how  to  obtain  E[£(t)]  from  E[£  (t)]  and 
E[z~(t)]. 

By  examination  of  Eq.(2.43),  observe  that 


E[Z (t) ]  =  E[z  (t-O]  -  E[z  (t-T -) ] 


(2.46) 


where,  by  the  conditional  expectation  theorem  [18,  p.208]. 


E[£+(t-Tr)l  =  Et  >z+(c)  [/(t-x^.)] 


(2.47) 


=  E  [E  +.  .[z+(t-T  )lx  -t  ]] 
T  z  (t)  —  — r  I— r  f 


E[£~(t-Tf) ] 


Vito1* 


yVwii'-vly'fH' 


(2.48) 


Taking  the  time  derivative  of  Eq.(2.46)  and  making  use  of  Eqs.(2.47)  and 
(2.48),  it  follows  that 

hjW)  -  ~£  E,  [Ez+(t)[£+(t-Tt)|lI-'rl] 


iEz-<t)^  If"  Tf!1 

— t  — 


(2.49) 


Let  f  (t  )  and  f  (x,)  denote  the  p.d.f.'s  of  T  and  X  _ ,  respec- 
— T  — t 

tively.  Carrying  out  the  expectations  with  respect  to  x  and  x,» 

— r  —r 

Eq.(2.49)  becomes 


E[Z(t))  =  ^ 


E[z  (t-T  )]f  (x  )dl 
0  —  r  x  r  r 


_d_ 

dt 


E[£_(t-Xf)lf  (xf)dx 
0  r  If  1  1 


E[z  (t-x  )]f  (x  )dx 
JO  —  r  x  r  r 


•  — 


EU  (t-xf)]fT  (xf)dxf 


=  E[z+(t)]*  fx  (t)  -  E[z~(t)]*fT  (t), 

-T  — f 


(2.50) 


For  convenience,  let  t  =  0  be  the  time  reference.  Integrating 
Eq.(2.50)  yields 


E[Z(t)  ]  =  E[Z(0)1  + 


{E[z  (0)]*f  (h)  -  E[z  ( •  ) ] * f  (  ) 

0  -r  -f 


(2.51) 


Note  that  Eq.(2.51)  expresses  E[Z^(t)]  in  terms  of  the  time  derivatives 
of  the  counting  signal  expectations  at  the  delay  element  input. 

Having  derived  Eq.(2.51),  attention  is  now  focused  on  obtaining 
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•  ^  •  •, 

E[*.  (t)3  and  E[z  (t)].  In  the  derivation  of  the  probability  expression 

theorem,  recall  that  Boolean  orthogonalization  was  used  to  express  the  • 

switching  function  for  z(t)  in  the  sum  of  nonoverlapping  implicants  form. 

This  form  leads  to  another  useful  result.  After  Boolean  orthogonali¬ 
zation,  the  switching  function  is  of  the  form  • 

—  ~  (XI » •  •  • » Xn)  =  _z^  m  _I^q  v  ^  v .  •  •  v  .1.^  ( £_i )  (2 . 52) 

Because' the  terms  are  nonoverlapping,  only  one  implicant  is  nonzero  • 

at  any  instant.  As  a  result,  an  equivalent  expression  for  the  switch¬ 
ing  function  is  obtained  by  replacing  each  OR  operation  in  Eq.(2.52) 
by  arithmetic  addition  as  shown  below:  • 


z 


-10  +  -11  +  •**  +  -l(Jl-l)* 


(2.53) 


In  addition,  each  implicant  is  of  the  form 


T  ac  A  •  A  • 

-lk  -lk  -2k  *•' 


•A  , 
— nk 


k  -  0,  1, 


(i-1) 


(2.54) 


where 

Ajk  £  (Xj,  Xj',1);  (2.55) 

Consequently,  an  equivalent  expression  for  I_  results  when  each 
AND  operation  in  Eq.(2.54)  is  replaced  by  arithmetic  multiplication 
as  follows: 


_ljk  =  (Ajk)  (^2lc^ ' "  *  nk^  ’  k  =  0 , 1 , . . . ,  (2,-1 ) . 


(2.56) 
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The  previous  observations  are  summarized  in  the  following  theorem: 

The  Arithmetic  Expression  Theorem 

Let  z_  =  gQ (X 1 , . . . , Xn)  be  a  switching  expression  of  the  n  switch¬ 
ing  variables  3U, . . .  ,Xn.  If  z_  is  expressed  as  a  sum  of  nonoverlapping 
implicants,  then  an  equivalent  arithmetic  expression  for  _z  is  obtained 
by  interchanging  each 

1)  OR  operation,  v ,  by  the  arithmetic  operation  of  addition, 

2)  AND  operation,  •  ,  by  the  arithmetic  operation  of  multiplica¬ 
tion. 

It  should  be  emphasized  that  performing  the  arithmetic  operations 
in  the  equivalent  arithmetic  expression  always  yields  the  correct 
binary  values  (i.e.  0  or  1)  for  z^. 


Example  2.5  Consider  the  switching  expression  z  =  go(Xl,X2,X3)  = 
(XI  v  X2')-X3  whose  truth  table  is  given  below. 


XI 

X2 

X3 

z 

0 

0 

0 

0 

0 

0 

1 

1 

0 

1 

0 

0 

0 

1 

1 

0 

1 

0 

0 

0 

1 

0 

1 

1 

1 

1 

0 

0 

1 

1 

1 

1 

A  sum  of  nonoverlapping  implicants  form  for  z  is 


z  =  XI  *X3  v  XI'  *X2 '  •  X3 .  (2.57) 

By  the  arithmetic  expression  theorem,  z  can  be  written  as 

z  =  (XI) (X3)  +  (X1')(X2')(X3).  (2.58) 

This  arithmetic  expression  can  be  shown  to  be  valid  by  checking 
each  row  of  the  truth  table.  E.g.,  for  the  sixth  row  in  which 
Xl=l,  X2=0,  X3=l ,  application  of  Eq.(2.58)  yields  z=(l) (l)+(0) (1) (1)  - 
as  predicted  by  the  truth  table.  V 

Application  of  the  arithmetic  expression  theorem  results  in  an 
arithmetic  form  for  z(t)  .  In  particular, 

l-l 

“  I  Iik  <t) 

k=0  1K 

'  jo  jll  V  (°  (2-59) 

The  time  derivative  of  _z(t), assuming  £(t)  is  mean-square  differentiable 
is  given  by 

.  X'_1  n  .  n 

£  (t)  =  ^  l  A  (t)  tt  Aft)  .  (2.60) 

k=0  i=l  1K  j=l 


Recall  from  Eq.  (2.55)  that 
A,,  e(xi,  xi',  1). 


(2.61) 


When  A,  =  Xi  ,  A,  =  Xi.  When  A,.  -  Xi’, 
— ik  —  — ik  —  — lk  — 


iuc  '  di  M  « 


(2.62) 


The  latter  result  follows  because  a  rise  transition  in  Xi_'  corresponds 

to  a  fall  transition  in  Xi  and  vice  versa.  Finally,  when  A^=  1, 

•  • 

_A^  =  Utilizing  the  above  and  collecting  terms  in  )^k  ,  Eq.(2.60) 

can  be  rewritten  as 


z(t)  =  l  Xk(t)Bk(t) 
k=l 


(2.63) 


Example  2.6  Given  the  switching  function 


z ( t )  =  Xl(t)  •  X2(t)  *  X3(t)  v  (XI (t)  •  X2  (t) ) '  •  X3'(t), 


(2.64) 


express  z  (t)  in  the  form  of  Eq.(2.63). 


For  convenience  the  time  symbol  t  and  the  underbar  symbol 
denoting  random  process  are  suppressed  in  the  ensuing  discussion 
whenever  the  meaning  is  clear.  The  first  step  is  to  express  z  as 
a  sum  of  nonoverlapping  implicants.  This  results  in 


z  =  XI  •  X2  •  X3  v  xi '  •  X2  •  X3 ' v  X2 '  •  X3 ' . 


(2.65) 


Next  the  arithmetic  expression  theorem  is  applied  to  yield 


z  =  (Xl)(X2)(X3)  +  (X1’)(X2)(X3')  +  (X2 ' ) (X3 ' ) . 


(2.66) 
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The  time  derivative  of  Eq.  (2.63)  is 

z  =  (XI)  (X2)  (X3)  +  (XI)  (X2)  (X3)  +  (XI)  (X2)  (X3) 

+  (-j-  [XI '  ] )  (X2)  (X3 ' )  +  (XI ' )  (X2)  (X3 ' ) 
at 

+  (Xl')(X2)(^[X3'])  +  (^[X2’])(X3’) 

+  (X2')  (—  [X3']). 


a 

a 


Substitution  of  Eq.(2.62)  into  Eq.  (2.67)  yields 

z  =  (X1)(X2)(X3)  +  (X1)(X2)(X3)  +  (X1)(X2)(X3) 

-  (XI) (X2)  (X3')  +  (XI ' ) (X2) (X3 ' )  -  (XI ' ) (X2) (X3) 

-  (X2)  (X3 '  )  -  (X2 '  )  (X3)  . 

Collecting  terms  result  in 

z  =  (X1)(X2)(X3  -  X3')  +  (X2)[(Xl)(X3)+(X3’)(xr-l)] 
+  (X3)  [ (X2) (Xl-Xl ' )  -  (X2 ' ) ] 

=  (XI) (Bl)  +  (X2)(B2)  +  (X3)(B3) 

where 

Bl  =  (X2) (X3-X3 ' ) 

B2  =  (XI)  (X3)  +  (X3')(xr-1) 

B3  =  (X2) (Xl-Xl ' )-(X2 ' ) 
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(2.67) 


(2.68) 


(2.69) 


(2.70) 


In  the  analysis  to  follow  it  is  necessary  to  have  analytical 
•+  •  — 

expressions  for  z_  (t)  and  z_  (t)  in  terms  of  the  switching  variables 

Xk  ,  k  «  l,...,n,  and  their  time  derivatives.  For  this  purpose,  an 

equivalent  arithmetic  expression  for  Bk(t)  must  be  obtained  such  that 

only  one  term  in  its  equivalent  expression  is  nonzero  at  each  instan.; 

and  the  value  of  that  nonzero  term  equals  the  value  of  Bk(t)  whenever 

it  is  nonzero.  This  form  is  referred  to  as  the  arithmetic  canonical  form. 

The  arithmetic  canonical  form  is  obtained  by  tabulating  as  in  a  truth 

table,  all  possible  combinations  of  the  values  for  the  variables  in  Bk 

(i.e.;  0,1),  along  with  the  corresponding  value  of  Bk  (i.e.;  -1,0,1). 

Having  this  table,  a  product  term  is  formed  for  each  row  by  writing 

Xi  as  a  factor  if  Xi  =  1  and  XI'  if  Xi  =  0*  Denote  the  product  term  for 

the  mth  row  in  the  table  by  Pk(m) .  In  forming  the  product  terms  in  this 

manner,  one  and  only  one  can  take  on  the  value  1  for  a  specified  set 

of  values  for  Xk;  k  =  l,...,n.  Specifically,  Pk(m)  =  1  only  for  the 

til 

values  of  Xk  specified  in  the  m  row  of  the  table.  It  follows  that 
the  expression  obtained  by  multiplying  Pk(m)  by  the  value  of  Bk  corres¬ 
ponding  to  the  row  and  summing  all  such  products  results  in  the 

desired  equivalent  expression  for  Bk.  Denote  the  value  of  Bk  corres¬ 
ponding  to  the  mt'"1  row  of  the  table  by  Bk(m)  .  The  arithmetic  canonical 
form  for  Bk  is  then  given  by 
2Y 

Bk  =  \  Bk(m)  Pk(m)  (2.72) 

m=  1 

where  ,  <  n-1  is  the  number  of  input  variables  Xk  upon  which  Bk  depends. 
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Example  2.6  (cont.)  Determine  the  arithmetic  canonical  forms 
for  Bl,  B2,  B3  in  Eq.  (2.70) . 

Following  the  procedure  outlined  above,  the  following  tables 
result  where  an  additional  column  has  been  included  for  tabulating  Pk(m) . 


Bl  * 

(X2) (X3-X3 ' ) 

B2  - 

(Xl)(X3)+(X3’)(Xr-l) 

B3  - 

(X2)  (Xl-Xl  '  )-(X2M 

m 

X2 

X3 

Bl(m) 

Pl(m> 

m 

D 

m 

XI 

X3 

B3(m)  P3(a) 

i 

0 

0 

0 

(X2 ') (X3' ) 

B 

0 

0 

0 

(X1-XX3) 

1 

0 

II 

-1 

(XI  *XX2') 

2 

J 

1 

0 

(X2 ' ) (X3) 

2 

0 

1 

0 

(X2’)(X3) 

2 

0 

1 

-1 

(XI ') (X2) 

I 

0 

-1 

(X2)  (X3 1 ) 

3 

1 

0 

-1 

(XI) (X31) 

3 

l 

0 

-1 

(XI)(X2’) 

4 

1 

1 

1 

(X2) (X3) 

4 

L 

L 

1 

(X1HX3) 

4 

1 

1 

1 

(X1)(X2) 

Using  Eq.  (2.72),  the  arithmetic  canonical  forms  for  the  coefficients 
Bk  become 

Bl  =  (0)(X2')  (X3')4-(0)(X2,)(X3)  +  (-l)(X2)(X3')+(l)(X2)(X3)  =  (X2)(X3)-(X2)(X3') 

B2  =  (0)(Xl')(X3')  +  (0)(xr)(X3  )  +  (-l)  (XI)  (X3')+(l)  (XI)  (X3)  =  (X1)  (X3)-(X1)  (X3’) 
B3  =  (-l)(xr)(X2,)+(-l)(xr)(X2)+(-l)(Xl)(X2’)+(l)(Xl(X2) 

-  (XI)  (X2)-(X1)  (X2  ' ) -  (XI ' )(X2)  -(XI 1 )  (X2  ' )  .  V  (2.73) 


In  the  arithmetic  canonical  form  for  Bk(t)  given  by  Eq.(2.72),  the 
terms  involved  in  the  summation  are  either  nonnegative  or  nonpositive. 

Denote  the  sum  of  all  nonnegative  terms  by  Bk-1- (t)  and  the  sum  of  all  non¬ 
positive  terms  by  -Bk  (t) .  Then 

Bk(t)  =  Bk+(t)-Bk  (t) .  (2.74) 


Note  that 


Example  2.6  (cont.)  Determine  analytical  expressions  for  z^"(t)  and 

•  _ 

—  (t)  in  terms  of  the  switching  variable  for  the  switching  function 
given  by  Eq. (2.64) . 

By  direct  substitution  of  Eqs.(2.76)  into  Eqs.(2.79),  one  obtains 
z+  =  (X1+)(X2)(X3)  +  (Xl~) (X2) (X3 ' )  +  (X2+) (XI) (X3)+(X2_) (XI) (X3 ' ) 

+  (X3+)(X1)(X2)  +  (X3~) [ (XI ') (X2)  +  (XI) 02')  +  (XI  * ) (X2 ' ) ] . 
z'  =  (Xl+) (X2) (X3 ' )  +(X1~) (X2) (X3)  +  (X2+) (XI) (X3 ' )  +  (X2_) (XI) (X3) 

+  (X3+) [ (XI') (X2)  +  (XI) (X2 ' )  +  (X1')(X2')]+(X3_)(X1)(X2) 

V  (2.8 

In  summary,  it  is  necessary  to  perform  the  following  steps 
in  order  to  obtain  analytical  expressions  for  £  (t)  and  z_  (t)  in  terms 
of  the  switching  variables  and  their  time  derivatives: 

1)  Express  z,(t)  as  a  sum  of  nonoverlapping  implicants  by  means 
of  the  Boolean  orthogonalization  procedure. 

2)  Obtain  an  arithmetic  expression  for  z_(t)  by  means  of  the 
arithmetic  expression  theorem. 

3)  Differentiate  with  respect  to  time  the  arithmetic  expression 


obtained  in  step  (2)  . 


4)  Use  Eq.(2.62)  to  eliminate  the  time  derivatives  of  all 

•  • 

complemented  variables  and  collect  terms  in  Xk  to  yield  z^(t)  in  the 
form  of  Eq. (2.63) . 

5)  Express  each  coefficient  Bk  by  its  arithmetic  canonical  form 
as  given  in  Eq.(2.72). 

_  t  •  -j_  •  _ 

6)  Separate  Bk  into  Bk  and  Bk  and  Xk  into  Xk  and  Xk  . 

7)  Use  Eqs.(2.79)  to  obtain  ^+(t)  and  z_  (t) . 

•  •  — 

Having  determined  analytical  expressions  for  z_  (t)  and  -l_  (t) ,  it 

•  +  •  _ 

is  now  possible  to  obtain  expressions  for  E[_z  (t)]  and  E[z^  ( t ) ]  which 
are  needed  in  the  evaluation  of  Eq.(2.51).  For  any  process  Y_( t )  which 
is  mean-square  differentiable,  it  follows  that  [18,  Sec.  9.6] 

E[Y(t)]  -  -^E[Y(t)]  -  E[Y(t)].  (2.81 

Recall  that  jlk  does  not  contain  Xk  as  a  factor.  Assuming  all  of  the 

input  switching  varibles  to  be  statistically  independent,  Bk+(t)  is 

•  •  —  — 

statistically  independent  of  Xk  (t)  and  (t) •  Similarly,  Bk  (t)  is 

*  "4"  *  — 

also  statistically  independent  of  Xk  (t)  and  Xk  (t).  Using  this  inde¬ 
pendence  and  the  property  stated  in  Eq.(2.81),  the  expected  values  of 
Eqs.(2.79)  become 

E  [_z +  ( t )  ]  =  l  (E[Xk+(t)]E[Bk+(t)]  +  E[Xk~(t)]  E  [  Bk  ~  ( t )  ]  > 
k*  1 

n  # 

E[_z“(t)]  =  l  (E[Xk  (t)] E[Bk“(t) ]  +  E[Xk  (t)  ]E[Bk  (t)  ]  }  (2.82 

k=l 


Removal  of  expected  values  of  complemented  variables  in  Eq.(2.82)  can 
be  accomplished  by  using  Eq.(2.38)  which  states  that 


E[Xk'(t)]  =  1  -  E[Xk(t)j. 


(2.83) 


Example  2 ■ 6  (cont . )  Utilizing  Eqs.(2.82)  and  (2.83),  obtain  the 
•  +  •  — 

expected  values  of  z_  and  z.  specified  by  Eqs.(2.80).  Assume  all 
switching  variables  to  be  statistically  independent. 

By  inspection,  it  follows  that 

E(z+(t)]  =  E[)Q+(t)]E[X2(t)]EQG(t)]  +  E[JO~(t)]  E[X2  (t)  ]  (l-E[X3  (t)  ]  } 

+  E[X2+(t)]E[JQ(t)jE[X3(t)]  +  E[X2~(t)  ]E[X1  (t)  ]  { 1-E[*3  (t)  ]  ) 

+  E[X3+  (t)  ]E[Xl_(t)  ]E[X2  (t)]  +  E[X3“(t)]  {  E[X2(t)  ]  (l  -  Epa(t)]} 
+  E-|Xl(t)j{l  -  E[X2  (t)  ] }  +  (l  -  E[Xl(t)]}  (1  -  E[X2(t)  ]  }} 


E[z"(t)]  =  E[n+(t)]E[X2(t)]{l  -  E[X3(t)  ]  }  +  E[n_(t)  ]E[}C2 (t)  ]E[X3(t)  ] 
+  E[X2+(t)]E[n(t)]{l  -  E[X3(t)>  +  E[X2~(t)]E[5Q(t)]E[X3(t)] 

+  E[X3+(t)]{E[X2(.t)Hl  -  E[Xl(t)]}  +  E[ja(t)]{l  -  E[X2(t)]} 

+  (1  -  E[Xl_(t)]}{l  -  E[X2(t)]}}  +  E[X3(t)]E[2a(t)]E[X2(t)] 

V  (2.84) 


The  previous  discussion  has  supplied  all  of  the  ingredients  necess¬ 
ary  for  the  evaluation  of  Eq.(2.51)  which  is  repeated  below. 


E[Z(t)]  =  E[Z(0)]  + 


+  •  _ 

(E[£  (0)]*f  (0)  -  E  [  z  (9)  ]  *f  (0)}d9 

0  -r  “  -f 


(2.85) 
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As  a  final  remark,  recall  that  jZ(t)  is  a  0,1  binary  random  process 
and  that  it  can  be  decomposed  into  counting  processes  as  in  Eq.(2.45) 
It  follows  that 

z(t>  =  z+(t)  -  z"(t)  .  (2 

The  expectation  of  Eq.(2.86)  is  given  by 

E[Z(t)]  =  E[Z+(t)]  -  E[Z“(t)].  (2 

Utilizing  the  interchangeability  of  differentiation  and  expectation, 
as  stated  in  Eq.(2.81), 

E[Z(t)]  =  E[Z+(t)]  -  E[Z(t)].  (2 

Taking  the  time  derivative  of  Eq.(2.85),  there  results 

E[Z(t)]  *  E[£+(t)]*f  (t)  -  E[z'(t)]*  fT  (t^  .  (2 

— r  -k 


Because  E[Z  (t)]»  E[£  (t)  ] ,  EU  (t)],  E[z  (t)],  (t) ,  f ^  (t)  are 

— r  — f 

each  nonnegative,  it  follows  that 


Example  2 . 7  Consider  the  physical  NAND  gate  whose  model  and  inputs 
are  shown  in  Fig.  2.7,  The  inputs  are  deterministic  and  are  assumed 


Xl(t). 
X2< 


(t>-  v 

w—l  r 


z(t) 


Xl(t) 


0.5 


3.5  5.0  6.5 


f  (t  ) 
ir  r 


UT 


T~r. 


X2(t) 


(b) 


1.0  3.0  4.5 

(c) 


7.0 


Fig.  2.7(a)  Circuit  model  (b)  p.d.f.'s  of  T_  and  _T  (c)  input  signals 


to  be  zero  everywhere  except  for  the  pulses  shown.  In  particular 
XI (t)  =  X2(t)  =  0  for  t  £  0.  As  a  result,  z(t)  =  1  for  t  £  0  and 

E[Z(0)]  =  Z(0)  =  1  .  (2.91) 

It  is  desired  to  compute  the  output  expected  value  E[Z_(t)],  using 

Eq.(2.85).  Following  the  methods  of  this  chapter,  the  output  expected 

•  •  — 

value  is  computed  in  two  steps.  First,  E[z  (t)]  and  E[z  (t)]  are 
computed.  The  convolutions  are  then  evaluated.  Note  that 


< 


z(t)  =  (Xl(t)  *X2(t))'  =  Xl'(t)  v  X2 ' (t)  . 


(2.92) 


In  terms  of  nonoverlapping  implicants,  z(t)  becomes 


z(t)  =  Xl'(t)  v  Xl(t) *X2' (t) . 


(2.93) 


The  arithmetic  expression  for  z(t)  is  given  by 

z(t)  -  Xl'(t)  +  Xl(t)  X2 • (t) .  (2. 

Differentiation  with  respect  to  time  yields 
z(t)  =  Xl'(t)  +  Xl(t)X2'(t)  +  Xl(t)X2*(t) 

=  -  Xl(t)  +  Xl(t)X2' (t)  -  XI (t)X2 (t) 

=  -  XI (t)  ( 1  -  X2'(t»-  Xl(t)X2(t) 

=  -  Xl(t)X2(t)  -  XI (t)X2(t) .  (2. 

Since  the  coefficients  in  Eq.  (2.95)  are  in  canonical  form, 

Bl+  =  0,  Bl“  =  X2,  B2+  =  0,  B2~  =  XI.  (2. 

Using  eqs.  (2.79), 

z+(t)  =  Xl~(t)  X2(t)  +  X2~(t)Xl (t)  (2. 

z~(t)  =  Xl+(t)  X2(t)+X2+(t)  Xl(t). 

Because  the  input  signals  are  deterministic,  E[Xk(t)]  =  Xk(t) . 

Equations  (2.97)  are  shown  graphically  in  Fig.  2.8.  In  order 
to  obtain,  E[_Z(t)],  the  derivative  of  the  expectation  of  the  delayed 
output,  convolutions  are  performed  as  in  Eq.(2.89). 


E[z+(t)]  =  z+(t) 


E[z  (t)]  =  z  (t) 


0  1.0 


Fig.  2.8  Derivatives  of  the  counting  processes  of  z(t) 
The  two  terms  of  Eq.  (2.89)  are  shown  in  Fig.  2.9. 


E[z+(t)]*frr(t) 


3.5  4.5 


7.0  8.0 


E[z~(t)]*fTf(t) 


1.5  3.0 


Fig.  2.9  Sketches  of  the  two  terms  in  Eq.(2.99). 

Using  the  results  of  Fig.  2.9  in  Eq.(2.85),  the  integration  yields  the 
two  terms  sketched  in  Fig.  2.10.  The  inital  condition  is  specified  by 
Eq .  (2.91) .  Performing  the  summations  in  Eq.(2.85),  the  final  result 
for  E[ji(t)]  is  shown  in  Fig.  2.11.  As  expected,  E[Z_(t)]  is  nonnegative 
and  varies  between  0  and  L.  To  evaluate  system  performance,  it  is 


2 


2.25  4.0 


6.25  7.5 


Fig.  2.12  Z(t)  when  1^=1  and  =  1.25 


Pr  {error  j  Z^)  =  0}  =  EtZCt^] 


Pr  {error  |  Z(t2)  =1}  =  1  -  EfZ(t2)] 


(2.98) 


It  follows  that  Pr  {error  |  Z(t)}  for  the  results  shown  in  Figs. 2. 11  and 
2.12  is  the  waveform  given  in  Fig.  2.13.  Note  that  the  error  does  not 
exceed  0.5.  Errors  larger  than  0.5  will  be  encountered  when  values  dif¬ 
fering  from  Tr  =  1,  Tf  =  1 . 25  are  specif ied . 


Pr  {error|z(t)} 


1.5  2.2  5  3.  3.5  4.  4.5  5.5  6.2  5  7.  7.5  8. 


Fig.  2.13  Conditional  error  when  T  =  1,  =  1.25. 


2 . 4  Tree-like  networks 

Thus  far,  the  discussion  has  dealt  with  analysis  of  a  single  logic 
block.  Given  the  expected  value  of  each  input  signal  and  assuming  the 
inputs  to  be  statistically  independent  it  was  shown  how  to  determine  the 
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expected  value  of  the  output.  Now  consider  a  tree-like  combinational 
network  of  logic  blocks.  A  tree-like  network  was  defined  to  be  a  network 
in  which  no  more  than  one  path  exists  from  each  input  to  every  output. 

As  a  consequence,  there  is  no  more  than  one  path  from  each  input  to  any 
point  of  the  network.  Given  that  the  inputs  to  the  network  are  statis¬ 
tically  independent,  it  follows  that  the  inputs  to  each  logic  block 
within  the  network  will  also  be  statistically  independent.  Also,  the 
expected  values  of  the  outputs  of  preceding  logic  blocks  provide  the 
expected  values  of  the  inputs  to  the  following  logic  blocks.  Therefore, 
by  repeated  applications  of  the  methods  developed  in  this  chapter,  it 
is  possible  to  extend  the  analysis  of  a  single  logic  block  to  the  analysis 
of  an  entire  tree-like  combinational  network.  The  extension  is  demon¬ 
strated  in  the  following  examples. 

Example  2.8  A  tree-like  combinational  network,  whose  logic  blocks  are 
assumed  to  include  pure  delay  elements,  is  shown  in  Fig.  2.14  along  with 
p.d.f.'s  of  the  delays.  For  clarity,  an  upper  case  letter  is  used  as 
a  subscript  to  the  symbol  _T  of  the  delay  to  denote  the  gate  type.  In 
particular,  T.  is  the  pure  delay  associated  with  AND  gates,  T  is 
the  pure  delay  associated  with  INVERTERS,  and  is  the  pure  delay 
associated  with  OR  gates.  All  inputs  are  deterministic  and  depicted 
in  Fig.  2.15.  It  is  desired  to  evaluate  the  output  expected  value. 


Fig.  2.14  (a)  network  for  Ex.  2.8, (b)  p.d.f  of  _r  (c)  p.d.f.  of  X  , 

i.  A 

(d)  p.d.f  of  x_ 


Xl(t) 


X2(t) 

X3(t) 

4 

1. 

t  n 

‘ 

L." 

10  0 


Fig.  2.15  Deterministic  input  signals  to  the  network  in  Fig.  2.14(a) 

The  notation  used  for  the  signals  in  this  example  is  identical  to  that 
illustrated  in  Fig.  2.1.  Beginning  at  the  network  input  and  successively 
applying  the  methods  derived  in  Sec.  2.2  results  in  the  expected  value 


of  the  network  output.  The  ideal  logic  outputs  for  the  AND  and  INVERTER 


Fig.  2.16  Sketches  of  (a)  yl(t)  and  (b)  y2(t). 

Using  Eq . (2 . 6) ,  the  output  expected  values  of  the  AND  and  INVERTER  gates 
are 

E[Y1  (t)]  =  yl(t)*fiA(t) 

E[Y2(t)]  =  y2(t)*fTl(t).  (2.100) 

These  are  sketched  in  Fig.  2.17. 


(a)  (b) 


Fig.  2.17  Sketches  of  (a)  E[Yl(t)]  and  (b)  E[Y2(t) ] . 

The  OR  gate  is  considered  next.  Since  YJ_(t)  and  Y2(t)  are  independent, 
it  follows  that  the  OR  gate  ideal  logic  output  expected  value  is  given  by 

E[  z  (t)  ]  =  E[Yl_(t)  v  Y2(t)]  =  E[Y_l_(t)  ]  +  E[Y2(t)  ] -E[Y1  (t)  ]E[Y2  (t)  ]  . 


Fig.  2.20  The  predicted  signal  Z(t),  when  Xj  =  0.375,  ta  =  1.0,  and  xQ  - 

The  conditional  probability  of  error,  for  this  case  Pr {error  |  Z(t)}, 
is  sketched  in  Fig.  2.21. 

Pr (error | Z (t) } 

1 

0.5 
0 

Fig.  2.21  Sketch  of  the  conditional  error  when  X  =  0.375,  x  *=  1.0,  and 

X  A 

x0  =  2.0 

V 

Example  2.9.  A  tree-like  combinational  network  is  shown  in  Fig.  2.22. 

All  logic  gates  are  assumed  to  include  discriminating  delay  elements 
and  the  network  inputs  are  assumed  to  be  statistically  independent. 

Also  assume  the  initial  conditions  E[Y1_(0)],  E[Y2(0)],  and  E[Z^0)]  are 
known.  It  is  desired  to  evaluate  the  output  expected  value,  E[Z(t)], 
in  terms  of  the  input  expectation;,  and  the  p.d.f.’s  of  the  gate  prop¬ 
agation  delays. 


jU(t) 

X2(t) 

X3(t) 

X4(t) 

Fig.,  2.22.  Tree-like  combinational  circuit  for  Ex.  2.9 

Once  again,  the  notation  of  Fig.  2.1  is  employed.  Using  the 

procedure  of  Sec.  2.3,  the  derivatives  E[yl+(t)],  E[yl  (t)  ] ,  E[y2+(t) j , 

•  __ 

and  E[y2  (t)]  are  obtained  as  follows: 

E[£l+(t)]  =  E[xT(t)]{l-E[X2(t)]}  +  E[X2“(t)]{l  -  E[Xl_(t)]} 

(2.103) 

E[yl~(t)  ]  =  Epa+ (t)  ]  { 1  -  E[X2(t)  ]  }  +  E[X2+(t)]{l  -  E[XA  (t)  ] } 

and 

E[y2+(t) ]  =  E[X3+(t)]E[X4(t)]  +  E[X4+(t) ]E[X3(t) ]  (2.104) 

E[^2“(t)]  =  E[3C3_(t)  ]  E[JC4  (t)  ]  +  E[Xi~(t)  ]E[X3(t)  ]  . 

•  "f*  •  — 

Carrying  out  the  convolutions  in  Eq.(2.90),  E[YJ_  (t)],E[Yl_  (t)  ]  , 

E[Y2+(t)],  and  E[Y2~(t)]  are  obtained.  Additional  upper  case  letters 
are  added  to  the  subscripts  of  the  propagation  delays  to  denote  the  gate 
type:  N-for  NOR  gate,  A-for  AND  gate,  and  0-for  OR  gate.  Thus, 

E[Yl_+(t)]  =  E[_^d+(t)]*fTrN(t)  ;  E[Yl_~(t)  =  E[jrl'(t)]*fxfN  (t); 

(2.105) 
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E[Y2  (t)]-  (t)  l*f xrA(t)  ;  E[Y2  (t)  ]  -  Et^'U)  ]*f xfA(t)  . 


Application  of  Eq.(2.85)  yields 


E[Yl (t) ]  =  E[Y1(0)]  +  {E[Yl+(6)]  -  E[Yl-(0)]}  d0 

•*0 


E[Y2(t) ]  =  E[Y2(0)]  + 


t 

f  (E[Y2+(0)]  -  E[Y2'(0)]jd0. 
0 


(2.106) 


(2.107) 


Now,  for  independent  inputs  Xi.(t)>  X2(c) »  ^(t) »  and  X4(t)  and 

independent  delays  T  ,  T  ,  t  and  t£a,  observe  that  the  two  signals 
~tN  —fN  — i rA  — f  A  ° 

Yl(t)  and  Y2(t)  are  statistically  independent.  Therefore,  the  same  pro¬ 
cedure  .can  be  repeated  to  obtain  the  network  output.  Specifically, 

E[£+(t)]  =  E[Yl+(t)]{l  -  E[Y2(t)]}  +  E[Y2+(t)]{l  -  E[Yl_(t)]} 

(2.108) 

E  [  z  ~  ( t )  ]  =  E  [  YJ_~  ( t )  ]  { 1  -  E[Y2(t)]}  +  E[Y2~(t)  ]  { 1  -  E[Yl(t)]} 


and 


E[Z(t)]  =  E[Z(0)]  + 


t 

{  E[z+(0)]*fTrO(0)  -  E[z  (0)  ] *f T f 0 (0)' } dG  . 


(2.109) 


Eq.  (2.109)  is  the  required  result. 


V 


Examples  2.8  and  2.9  demonstrate  the  simplicity  of  the  systematic 
method  for  evaluation  of  the  output  expected  value  of  a  tree-like  net¬ 
work.  The  method  becomes  considerably  more  complicated  when  networks 
with  reconvergent  fanouts  are  considered.  This  is  discussed  in  the 
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next  chapter. 


3.  ANALYSIS  OF  NETWORKS  WITH  RECONVERGENT  FANOUTS 


3.1  Introduction 

As  described  in  Chapter  1,  a  combinational  network  is  said  to 
contain  reconvergent  fanout  if  more  than  one  path  exists  from  an 
input  to  an  output,  (i.e.,  a  network  has  reconvergent  fanout  when  two 
or  more  branches  originate  from  the  same  point  and  merge  back  together 
at  some  other  point  of  the  network).  Example  3.1  depicts  a  network  con¬ 
taining  reconvergent  fanout. 

Example  3.1  Observe  the  combinational  network  shown  in  Fig.  3.1. 


Fig.  3.1  A  network  with  reconvergent  fanout. 

Note  that  the  signal  path  for  XI  splits  at  Y1  and  reconverges  at  gate 
E  (i.e.,  two  paths  exist  from  XI  to  the  output  Z) .  Also, the  signal  path 
for  X2  splits  between  gates  A  and  B  and  reconverges  both  at  gate  D  and 
gate  E  (i.e.,  three  paths  exist  from  X2  to  the  output  Z) .  Observe  that 
X3  has  only  one  path  to  the  output .  V 

To  demonstrate  the  difficulties  that  arise  in  the  analysis  of 
networks  containing  reconvergent  fanout  consider  the  following  simple 
examp le . 


Example  3.2.  A  simple  combinational  network  is  shown  in  Fig.  3.2. 

The  inputs  XI  (t)  and  X2(t)  are  assumed  to  be  statistically  independent. 
Observe  that  the  signal  X2(t)  has  two  paths  to  the  output.  Thus. by 
definition, this  network  contains  reconvergent  fanout.  For  simplicity. 


XI  (t) 


X2(t) 


Fig.  3.2  Circait  for  Ex.  3.2 


suppose  all  the  gates  in  Fig.  3.2  include  pure  delays  (i.e  x^  =  x^  =  x) . 
Following  the  notation  of  Ch.2,  denote  the  pure  delay  of  the  AND  gate 
by  and  the  pure  delay  of  the  OR  gate  by  x^.  The  switching  expression 
for  the  output  is  given  by: 


Z(t)  =  XlCt-T^  v  XlCt-^J-MCt-l^) 


(3.1) 


To  obtain  the  sum  of  nonoverlapping  iraplicants  form  note  that  the 
random  variable  X2 (t-XQ)  differs  from  the  random  variable  X2 (t-T  -  x_  q) 
since  they  are  taken  at  different  time  instants.  As  a  result,  X2(t-X^) 
and  X2(t-x_  -x_)  must  be  treated  as  two  different  switching  variables. 
Therefore,  the  sum  of  nonoverlapping  implicants  is  given  by 


Z(t)  =  X2(t-x  )  v  Xl(t-x  -t  )-X2(t-T  -t  )*X2'(t-T^)  . 


(3.2) 


Application  of  the  probability  expression  theorem  to  Eq. (3.2),  using 
independence  of  XI (t)  and  X2(t),  yields: 


E[Z(t)]  =  E[X2(t  -T^]  +  E[n(t-TA-T0)]E[X2(t-tA-i<))x2-(t-i0)].  (3.3) 


Note  that  the  second  term  of  the  right  side  of  Eq.(3.3)  includes  the 
expectation  E[X2 (t-TA-TQ)X2 '  (t-T^)],  called  the  second  order  joint 
moment  of  X2(t)  and _X2'(t).  The  above  second  order  joint  moment  can  be  fur¬ 
ther  simplified  into  E[X2  (t-T^-Tg)  ]  — E [_X2^ (t-T^-T^)X2  (t-T^)  ]  .  The  expecta¬ 
tion  EfX2(t-x  -T„)X2(t-T  ) ]  is  referred  to  as  a  second  order  moment 
—  —A  —0  —  —0 

of  X2(t)  .  V 

In  general,  the  8,^  order  moment  of  Xk(t)  is  defined  to  be  [18, p.297] 


^Xk  xk^  1  *  *  *  *  * ^8,^  ~  E [ Xk (t  ^  )  » • . Xk ( t ]  , 


(3.4) 


while  by  definition,  the  £  order  joint  moment  of  the  m  different  processes 
XI (t) , . . . ,Xm(t)  is 


^Xl.-.Xl  X2...Xm(tl’*’*,t:£l’  t£l+I,'*”t£2,'”’t£) 


=  E[2Q(t1)...n(t£1))C2(t£1+1)...2C2(tjl2)...Xm(tJl)],l  <_£1<£2<...<£m=£. 


(3.5) 


In  Example  3.2,  it  was  observed  that  due  to  the  different  delays,  the 


two  paths  between  X2  and  the  output  resulted  in  an  expression  for  the 
expected  value  of  the  output  which  includes  a  term  involving  the  second 


'->'3 
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order  moment  of  X2.  .The  order  of  reconvergence  for  a  signal  yk  is  defined 
to  be  the  number  of  paths  having  different  delays  from  the  node  ^k  to  the 
output  node  Z.  It  Is  now  shown  that  the  order  for  moments  of  required 
in  E[Z]  may  equal  the  order  of  reconvergence  for  yk. 

Assume  the  last  gate  preceding  the  output  node  Z_  is  an  AND  gate 
containing  m+q-1  delayed  inputs  as  shown  in  Fig.  3.3.  The  (m-1)  inputs 
Yi(t-T  11),  i  2,...  ,m,  have  only  one  path  from  the  node  Yi  to  the  output 
node  Z,  where  the  delay  associated  with  the  path  from  node  YJL  to  the  input 
of  the  AND  gate  is  denoted  by  jr^.  The  signal  at  node  Yl,  on  the  other 


Yl_(t-Tu) 


I2(t-T21) 


Ym(t-i  ) 
—  —ml 


Z(t) 


Fig.  3.3.  AND  gate  containing  nrt-q-1  delayed  inputs. 


hand,  has  reconvergence  of  the  order  q  (i.e.,  there  are  q  paths  with 
different  delays  from  the  node  H  to  the  output  node  Z) .  The  delays 
from  the  node  U_  to  the  corresponding  inputs  of  the  AND  gate  are  denoted 
by  In  , . . .  ,x_iq.  All  signals  jH.(t) .... ,Ym(t)  are  assumed  to  be  statis¬ 
tically  independent.  The  output  I(t)is  given  by 


Z(t)  =  YJ_  ( t  -j_i  i )  a  ...  a  Yl  (t-Tiq)AY2(t-T21)A  .  .  .AYm(t- tml) 


(3.6) 


Observe  that  Eq.(3.6)  is  in  the  form  of  a  sum  of  nonoverlapping  implicants 
where  the  sum  contains  a  single  term.  Applying  the  probability  expression 

theorem, 

E[Z(t)]=  E[Yl(t-T11)...Yl(t-Tlq)Y2(t-T21)...Ym(t-Tml)].  (3.7) 

Utilizing  the  statistical  independence  among  YJ_(t)  ,...  ,Ym(t)  ,  Eq.(3.7) 
is  written  as: 

E[Z(t)]  =  E[Yl(t-2n)...Vl(t-Tlq)]  E[Y2(t-T2i))...EtYm(t-TTOl)]. 

(3.8) 

It  is  observed  that  the  output  expected  value,  E[Z(t)],  includes 
til 

the  q  order  moment  of  YJ^t)  .  But  q  is  also  the  order  of  reconvergence 
of  the  signal  Y_l(t).  Thus,  in  this  case,  the  order  required  for  the 
moments  of  YJ_(t)  in  the  expected  value  of  the  output  signal  Z^(t)  equals 
the  order  of  reconvergence  of  Y1 (t)  . 

However,  the  order  required  for  the  moment  of  Yk(t)  does  not  necess¬ 
arily  equal  the  corresponding  order  of  reconvergence.  It  may  be  lower  if 
some  of  the  paths  are  maintained  at  a  constant  value.  This  is  illustrated 
by  the  following  example. 

Example  3.3.  Consider  the  network  shown  in  Fig.  3.4. 


Fig.  3.4.  Network  with  one  input  held  at  1. 


Observe  that  there  are  two  different  paths  between  the  input 
node  and  the  output  node  Z.  However,  since  one  of  the  inputs  to 
the  OR  gate  is  maintained  at  the  value  1,  the  output  Yl_(t)  is  also  held 
at  the  value  1.  Therefore,  the  output  j2(t)  is  affected  by  only  one  of  the 
paths  Involving  XI  (t).  In  particular,  denoting  the  delays  of  gates  A1 
and  A2  by  and  ,  respectively,  the  output  £(t)  is  given  by 

Z(t)  =  XI X^t-T^-T^).  (3.9) 

Assuming  the  inputs  XI (t)  and  X2(t)  are  statistically  independent 
and  using  the  conditional  expectation  theorem,  the  expected  value  of  Z^t) 
is  given  by 


oo  oo 

* 

Eli<t,]  ‘  l  Jo  EXi(t)lH(t-TA1-TA2)) 

(3.10) 

EX2  (t)  (t_TArTA2)  ■* 

fT_Al(TAl)fTA2(TA2)dTAl  dlA2 

Observe  that  only  first  order  moments  arise  in  Eq.(3.10)  even  though 

the  order  of  reconvergence  for  XKt)  is  two.  V 

When  discriminating  delay  elements  in  the  network  are  considered, 

the  signals  must  be  decomposed  into  their  counting  signals.  As  a 

result,  four  different  signals  may  be  encountered  for  every  signal 

•4*  •  — 

Yk(t),  in  evaluation  of  the  derivatives  z_  (t)  and  _z  (t)  .  These  are 
4*  *  4*  •  — 

Yk  (t) ,  Yk  (t) ,  Yk  (t) ,  and  Yk  (t) .  If  the  order  of  reconvergence  for 

Yk  is  q,  then  the  analysis  for  E[Z,(t)]  may  require  joint  moments  of  orders 

4  —  *  4*  •  — 

up  to  q  involving  Yk  (t) ,  Yk  (t) ,  Yk  (t)  and  Yk  (t) .  These  terms  can 
get  quite  large  in  number  since  the  number  of  different  ways  to  place 


4  objects  in  ■,  cells  with  repetition  is  given  by  [22, p. 909] 


( i+3) ! 
y!  3! 


(3.11) 


Example  3.4.  For  ,=2  there  are  — — —  =  10  possible  second  order 

joint  moments  (correlation  functions)  involving  the  counting 
signals  of  Yk(tj  and  their  derivatives. 

Specifically,  the  ten  possible  correlations  function  for  Yk(t) 
are  listed  below: 


RYk+Yk+(tl,t2)  ’  N'kSfk  ^l*^’  RYk+Yk+(t 1* ’V  ’  RYk+Yk  (ti*C2)»  Sk'Yk  (tl,t:2^ 


RYk_Yk+(trt2)  ,RYk~Yk  (tl’t2),RYk'fYk+(t:l,t2)  ’  ^k^k'^l  ,t2)  ’  ^k^Yk  (ti*t2) 


(3.12)  V 


Methods  for  reducing  the  number  of  needed  higher  order  moments  are 
discussed  in  Sec.  3.3.,  where  discriminating  delay  elements  are  considered. 
As  in  Chapter  2,  the  simpler  case  of  pure  delay  is  discussed  first. 

3.2.  Pure  Delay 

A  typical  logic  block  with  a  pure  delay  element,  i.e.  Tr  =  =  T> 

is  illustrated  in  Fig.  3.5  (same  as  Fig. 2. 3). 
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Xl(t). 


xaXt), 


Ideal 


Logic 


Circuit 


Delay  element 


Fig.  3.5  A  typical  logic  block  with  a  pure  delay  element. 


Again,  for  ease  of  discussion,  a  single  output  Z^(t)is  considered. 
The  inputs  Xl(t),...,Xn(t)  are  assumed  to  be  stochastic  0,1  binary 
processes  which  are  not  necessarily  statistically  independent. 

When  the  output  signal,  jl(t)  is  used  as  an  input  to  a  second  net¬ 
work  having  reconvergent  fanout,  the  higher  order  moments  of  Z^(t) 
are  required  for  evaluation  of  the  output  expected  value  of  the  second 
network.  Therefore,  in  this  section  it  is  shown  how  to  obtain  higher 
order  statistics  of  ,Z(t)  in  terms  of  the  higher  order  statistics  of  the 
inputs  Xl_ (t) , .  .  .  ,Xn(t)  and  the  p.d.f.  f^(t). 

For  the  pure  delay  case  the  sample  functions  of  Z(t)  are  delayed 
replicas  of  the  sample  functions  of  z(t).  Specifically, 


Z(t)  =  z  (t  -  t) . 


Therefore,  the  ■  order  moment  of  Z(t)  is  given  by 


E[Z(t ,)...  Z(t„) ]  =  E[z(t  -I)...  z(t  -  l) 

L  X.  i  —  —  K  “* 


(3.13) 


(3.14) 


69 


The  expectation  in  Eq.(3.14)  is  ts'-;..  with  respect  to  both  z(t) 
and  T .  Employing  the  conditional  expectation  theorem  [18,  p.208]  and 
using  the  notation  of  Sec.  2.2,  one  obtains 


RZ. . .Z(tl5t2’ • ’ ‘ 'H ) 


=  E[Z(t1)...Z(t£)]  =  ET[E2U(t1-T)...£(t£  -  T )  T.  =  T  ]] 


E  [z(t  -  T)...z(t  -  T)  T  =  T]  f  (t)  dx 

Z  —  i.  —  —  Z  —  —  T 


E  [  2  (t  -T)  ...  z(t  -  T)]f  (T)dT, 
—  1  —  k.  T 


(3.15 


The  following  simple  example  illustrates  this  evaluation. 


Example  3.5  Suppose  z(t)  is  deterministic  and  consists  of  a  single 


pulse  between  0  and  2.  Let  the  delay  T  be  uniformly  distributed  between 


0.5  and  1.  From  Eq.  (3.15),  the  output  autocorrelation  is  given  by 


Rzz(trt2)  =  £[z(ti  -T)z(t  -T)lfT(T)  dT  • 


(3.16 


Because  z(t)  is  deterministic  and  T  is  the  variable  of  integration 


E[z(t1-T)z(t  -  t)]  =  z(t  -  T)z(t2  -  T)  . 


(3-17 


Hence , 


RZZ(C1’C2)"  J  z(tl  "  T^z(t:2  "  T)fT(T>dT 


z(ti  -  x)z(t2  -  i)dx. 


(3-18 
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z(t),  fT(t),  the  input  autocorrelation  function  and  the  output 
autocorrelation  function  are  sketched  in  Fig.  3.6. 


Fig.  3.6.  (a)  signal  z(t)  (b)  p.d.f.  f^Cx)  (c)  Input  autocorrelation 
function  and  (d)  output  autocorrelation  function. 

The  output  autocorrelation  function  is  seen  to  be  a  shifted  and  distorted 
version  of  the  input  autocorrelation  function  due  to  the  random  delay.  V 
In  general,  knowledge  of  E[z(t  -  x)...£(t£-  x)]  is  necessary  in  order 
to  evaluate  F.q .  (3.15).  For  this  purpose,  it  is  convenient  to  derive 
EU(t  j)  .  .  .£(t^)  ]  in  terms  of  the  higher  order  moments  of  the  input 
switching  variable  Xl (t) , . . . ,Xn(t) .  It  is  now  shown  that  the  arithmetic 
expression  theorem  from  Ch.  2  is  also  applicable  here. 


Fig.  3.7  AND  gate  with  inputs  £(t  ^ ,z_(t  ^)  . 

Let  each  switching  function  _z(t^)  »  k  =  be  expressed  as  a  sum  of 

nonoverlapping  implicants.  Observe  that  performing  an  AND  operation  between 
a  set  of  nonoverlapping  implicants  and  other  switching  variables  generate 
another  set  of  nonoverlapping  implicants.  It  follows  that  application 
of  the  distributive  law  to  the  output  _z(tpA  £(t^)A  •  •  .^(tp  )  yields  a  sum 
of  nonoverlapping  implicants.  Use  of  the  arithmetic  expression  theorem 
from  Ch.  2  then  results  in  an  arithmetic  expression  for  the  output.  The 
output  expected  value  E[z_(t  ^ )  .  .  ._z  (t^,  )  ]  is  next  obtained  by  recognizing 
that  the  expectation  of  a  sum  equals  the  sum  of  the  expectations.  Finally, 
the  expression  can  be  simplified  by  taking  advantage  of  statistical  inde¬ 
pendence  where  applicable. 

A  useful  tool  in  the  evaluation  of  Ef z (t , ) . . . z (t  ) ]  is  to  recognize 
an  equivalent  approach  which  is  to  replace  each  switching  function  z(t  ), 

K 

k  =  1 . 7. ,  by  its  arithmetic  expression  and  then  apply  the  distributive 


law  on  the  arithmetic  expression. 


Example  3.6  Determine  the  third  order  moment  of  z^(t) ,  E[z(ti)z(t2)£(t3) ] , 
for  z^(t)  =  lil(t)  v  X2(t)  .  The  inputs  XI  (t)  and  X2(t)  are  assumed  to  be 
statistically  independent. 

An  arithmetic  expression  for  z^(t)  is  given  by 

£(t)  »  5Q(t)  +  XI ' (t)X2 (t) .  (3.19) 

Consequently, 


z(t1)z(t2)z(t3)  =  [xi(t1)+ n'(t1)X2(t1)]{x^(t2) 

+  n'(t2)x2(t2)] [Xi(t3)  +  xL'(t3)x2(t3)]  # 


(3.20) 


Application  of  the  distributive  law  results  in 
z.(t1)z(t2)z(t3)  =  [Xl(t1)Xl(t2)Xl(t3)+  Xl(t1)Xl.,(t2)Xl(t3)X2(t2) 

+  XV  (tj)Xl_  (t2)  XI (t 3)  X2(tx)+  XI’  (t1)Xl'  (t2)Xl(t3)X2(t1)X2(t2) 
+Xl(t1)Xl(t2)Xl'(t3)X2(t3)+  Xl(t1)n,(t2)Xl’(t3)5(2(t2)2C2(t3) 

+  Xl,(ti)Xl(t2)Xl'(t3)X2(ti)X2(t3) 

+  Xl'(ti)Xl,(t2)Xl'(t3)X2(ti)X2(t2)X2(t3)  (3.21) 


Taking  the  expectation  of  Eq.(3.21)  and  using  statistical  independence 
when  possible 


E [— ( c i )— (t2)— ( c 3)  ]  ~  E[xi(t1)iu(t2)x]_(t3)  ]+  EtxUt^jU 1  (t2)xi(t3)  ]E[x2(t2)  ] 
+  E[n'(t1)2y.(t2m(t3)]E[X2(t1)]  +  E[}Q’(t1)n,(t2m  (t3)Rx2x2(tl,t2) 


+  E[n(t1)2a(t2)2a' (t3)]E[x2(t3)]  +  e[xi  (t^xi '  (t2)xr ’  (t3>]Rx2X2(t2,t3) 

+  E[X1'  (t1)}Q(t2)2a'  (t3)]Rx2  X2(tl,t3) 

+  E[Xl'(t1)Xl'(t2m'(t3)]E[X2(t1)X2(t2)2(2(t3)].  (3.22) 

As  a  final  point,  note  that  complemented  variables  in  Eq .  (3.22)  could 
have  been  eliminated  by  using  in  Eq.(3.20)  the  property  that 


Xk'  (t)  =  1  -  Xk(t) .  V  (3.23) 

In  general,  networks  containing  reconvergent  fanout  also  contain 
tree-like  subnetworks.  The  tree-like  subnetworks  can  be  analyzed  using 
the  procedures  described  in  Sec.  2.2,  The  higher  order  moments  which 
arise  due  to  the  reconvergences  are  obtained  as  explained  i.  this 
section.  The  procedure  is  illustrated  in  the  following  network  example. 

Example  3.7  Consider  the  combinational  network  shown  in  Figure  3.8.  # 

All  gates  are  assumed  to  contain  pure  delay  elements.  Also,  the  inputs 

Xk(t),  k=l,2,3  are  assumed  to  be  statistically  independent.  Since  there 

are  two  paths  to  the  output  for  X2(t)  and  X3(t) ,  the  network  contains  # 

a  reconvergent  fanout.  It  is  desired  to  evaluate  the  expected  value 

of  the  output  E  [ Z  ( t )  ]  . 
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Fig.  3.8.  Combinational  network  with  a  reconvergent  fanout. 


The  model  for  the  physical  network  in  Fig.  3.8  is  illustrated 
in  Fig.  3.9.  Observe  the  ideal  logic  output  signals  ^_l_(t)  ,  ^(t),  and 
z(t)  and  the  pure  delay  elements  with  random  delays  and  x^. 


X3(t) 

Fig.  3.9.  Model  for  network  in  Fig.  3.8. 

The  expectation  of  the  ideal  logic  output  ^2(t)  is  given  by 


E[^2(t)]  =  E[X2(t)]E[X3(t)] .  (3.24) 

Using  Eq.(2.6)  the  expected  value  of  the  gate  output  Y2(t)  is 

E  [  Y_2  (t)  ]  =  E[£2(t)]*fTA(t)  =  {E[X2(t)]E[X3(t)l}*fxA(t).  (3.25) 


The  expected  value  of  YHt)  is  also  readily  obtained.  Noting  that 
XI  (t)  and  Y2(t)  are  statistically  independent  the  expected  value  of 
^l_(t)  is  given  by 


E[yl (t) ]  -  E[Xl(t)]  +  E[Y2(t)]-E[Xl(t)]E[Y2(t)]. 


(3.26) 


Once  agin,  by  means  of  £q.(2.6). 


E[Yl(t) ]  =  E[2l(t)]*f  (t) 


E[Xl(t)]*f  (t)  +  E[Y2(t)]*f  (t)  -  (E[Xl(t)]E[Y2(t)]}*f  (t) 


(3.27) 


with  respect  to  the  last  AND  gate,  the  expected  value  of  z_(t)  is 
given  by 


E [ z  (t) ]  =  E[Y1 (t)Y2(t) ] . 


(3.28) 


Note  that  _Yl_(t)  and  Y2(t)  are  not  statistically  independent.  The  arith¬ 
metic  expression  for  Yl(t)  is  given  by 


Yl(t)  =  2LL(t“I0)+  Ii(t  -I0)-2y_(t-T<))Y2(t-T0).  (3.29) 


Substituting  Eq.(3.29)  into  Eq.(3.28)  and  using  the  independence  of 


XI  (t)  and  Y2(t)  one  obtains 


E[z(t)]=  E[Xl(t-T<))Y2(t)+Y2(t-T0)Y2(t)  -  XI  (t-^Wi (t-l^m* (t)  ] 


=  E[X_l(t-T0)]E[Y2(t)]  +  R^2  Y2(t,t-T0)-E[Xl(t-x0)]RY2  Y2(t,t-T0) 


(3.30) 
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Using  the  conditional  expectation  theorem  [18]  to  average  over  the 
delay  Jq,  Eq.(3.30)  becomes 

OO 

E[z(t)]=  {E[Xl(t)]*fx0(t)}E[Y2(t)]  +  Ry2  (t  ,t-x)  fx()(x)dx 

OO 

* 

E[}Q(t-T)RY2  Y2(t,t-T)fT0(T)dT  . 


Observe  that  evaluation  of  Eq.(3.31)  requires  knowledge  of  the  auto¬ 
correlation  function  for  Y7(t) .  By  definition  the  autocorrelation 
function  of  the  ideal  logic  output  ^(t)  is 


RZ2  l2(t,£'T)  =  E[y2(t)l2(t-T)]  =  E[X2(t)X3(t)X2(t-x)X3(t-x)] 
=  ^2  X2(t,t-T)RX3  X3(t,t-T)  ' 


Making  use  of  Eq.(3.l5),  the  autocorrelation  function  of  Y2(t)  is 


(3.31) 


(3.32) 


&Y2  Y2(t,t-x)  =  j  ^2(t-a,t-x-a)fTA(a)da 


I  *S(2  X2  x3(t-a,t-x-a)fTA(a)da. 


(3.33) 


Substituting  into  Eq.(3.31),  Eqs.(3.25)  and  (3.33)  for  E[Y2(t) ] 

and  Ry2  Y2(t,t-x),  respectively ,  an  expression  for  E[:z(t)]  is  obtained 

in  terms  of  the  inputs  ^(t) ,  X^t)  »  and  X3(t)  • 

Finally,  from  Eq.(2.6)  E[^(t)]  is  given  by 


E[Z(t)]  =  E[z(t)]*fxA(t)  . 


(3.34) 


*.  •.  -V  * 
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Note  that  in  Eq.(3.34)  first  and  second  order  moments  for  the  recon¬ 
verging  signals  X2(t)  and  X3(t)  are  needed  while  only  the  mean  of  XI (t) 
is  required.  V 

In  summary,  the  procedure  for  obtaining  the  ouput  expected  value 
of  a  combinational  network  with  pure  delay  elements  is  performed  by 
progressing  from  the  inputs  to  the  output,  logic  block  by  logic  block. 

For  the  tree-like  portion,  the  expectation  of  a  block  output  is  evalu¬ 
ated  in  terms  of  the  expectations  of  the  block  inputs  and  the  p.d.f.  of 
the  corresponding  delay.  When  a  reconvergent  fanout  is  encountered, 
the  ideal  logic  output  signal  must  be  expressed  explicitly  in  terms  of 
the  reconverging  signals  using  the  arithmetic  expression  theorem.  Taking 
the  expected  value  of  the  arithmetic  expression  results  in  higher  order 
moments  as  illustrated  in  Eq.(3.30).  The  higher  order  moments  can  always 
be  expressed  in  terms  of  the  higher  order  moments  of  the  primary  inputs, 
as  illustrated  in  Eq.(3.33).  As  was  done  in  Eq.(3.31),  it  is  often  use¬ 
ful  to  employ  the  conditional  expectation  theorem.  Finally,  the  output 
expected  value,  E[Z^(t)],  is  obtained  by  convolving  E[z^t)]  with  the  p.d.f. 
of  the  output  delay,  as  in  Eq.(3.34). 


3.3.  Discriminating  Delay 

The  model  of  a  typical  logic  block  with  a  single  output  and  a  dis¬ 
criminating  delayCx^  f  t^)  is  illustrated  in  Fig.  3.10  (same  as  Fig.  2.5). 

Xl(t) 


Xn(t) 

Fig.  3.10.  The  model  for  a  typical  logic  block  with  T^_  4  i  ^ . 


It  was  pointed  out  in  Ch.l  that  in  such  a  case  the  signal  Z^(t) 
must  be  decomposed  into  the  difference  of  two  counting  signals  as 
given  by 

Z(t)  =  Z+(t)  -  Z"(t).  (3.35) 

However,  a  direct  derivation  for  Z_+( t)  and  jZ  (t)  in  terms  of  the 

input  signals  XI (t) , . . . ,Xn(t)  and/or  their  counting  signals  is  not 

known.  As  in  Ch.2,  the  difficulty  is  overcome  by  introducing  the 
•  4-  •  — 

time  derivatives  Z_  (t)  and  Z_  (t) .  When  reconvergence  is  encountered, 

4*  —  •  ^  •  _ 

the  higher  order  joint  moments  of  Z^  (t) ,  Z^  (t),  Z  (t)  and  Z  (t) 
are  required.  Their  derivation  is  presented  in  this  section.  Because 
a  very  large  number  of  such  moments  may  arise,  an  efficient  procedure 
for  performing  the  network  analysis  is  also  developed. 

Suppose  the  output  Z^t)  of  the  logic  block  shown  in  Fig.  3.10  is 
used  as  an  input  to  a  second  network  containing  reconvergent  fanout . 

Analysis  of  the  second  network  requires  higher  order  joint  moments  of 

4  —  *4"  •  —  til 

Z  (t) ,  (t),  Z  (t)  and  Z_  (t) .  Denote  the  £  order  joint  moment  by 

E[_Z  {ty)...Z  (t^j)_Z  (t£l+1)  •  •  -Z  ^*"£2+1  ^  (C£3^-  ^£3+1^  " ^£^  ; 

0  <  £1  _<  £2  _<  £3  <  £. 

Recall  that  Z + (t)  is  the  rise  transition  counting  process  and 

•  +  - 

Z  (t)  is  its  derivative.  Also  Z_  (t)  is  the  fall  transition  counting 
•  — 

process  and  Z  (t)  is  its  derivative.  As  discussed  in  Ch.  2,  these 
processes  can  be  expressed  in  terms  of  the  delay  element  input  processes 


£  (t) ,  £  (t) ,  z  (t) ,  and  £  (t)  as  follows: 

Z+ (t)  =  £+  (t-Tr)  ,  Z+ (t)  =  z+ (t-Tr) 

(t)  =  £  (t-x_f)  ,  Z  (t)  =  z  (t-Tf) .  (3.36) 

Substituting  Eqs.(3.36)  into  the  order  joint  moment  yields 

e[z  (t^-.-X  (t£1)z  (t£1+1)...z  (t£2)i  (t£2+1>...z  (t£3)z(tn+1)...|  (t£)] 

E[£  (t^Xp)  ••  -z.  (t 2.2+1“  -r^ 

**•-  ^W+l'-f^"*-  (3.37) 


Application  of  the  conditional  expectation  theorem  to  the  right  side 
of  Eq.(3.37)  results  in 


E[_Z  (t^.-.Z,  (t£l)X  (t£1+1)  •  • -X  (t£2)—  ('tZ2+l')"'~  ^23^-  ^23+1^*"'-  (t2^ 


J0  J0 


Et  £+< q-Tr>  •  •  ■£  ('li-'rli  KlM-'ft  •  •  ■£"«C£2-Tf)i+  <H2+rV 


i  (tt3'Tt>^tl3+r'I£)---i!trT£)1£TrTf(Tr>Tf)dTrdT£ 


(3.38) 


In  general,  whenever  the  binary  signal  Z^(t)  or  its  derivative  ,Z(t)  are 
directly  involved  in  an  expression,  they  can  always  be  decomposed  into 


Z(t)  =  Z  (t)  -  Z  (t)  or  Z(t)  -  Z.  (t)  -  _Z  (t),  respectively. 
Chapter  2  discusses  how  to  obtain  the  derivatives  of  the 


counting  processes  at  the  output  of  an  ideal  logic  circuit.  Therefore, 


terms  in  Eq.(3.38)  involving  _z  (t)  and  z_  (t)  can  be  obtained  by  inte¬ 
nt 


grating  with  respect  to  time  the  derivatives  _z  (t)  and  z_  (t)  .  The 
integrals  are  assumed  to  exist  in  a  mean  square  sense.  (This  simply 
requires  that  the  respective  autocorrelation  functions  be  integrable  over 


the  region  of  integration.)  Therefore,  selecting  t  =  0  to  be  the  time 

+  , 


reference,  z  (t^  -  x^)  and  £  (t^  -  x_^)  may  be  expressed  as 


z+(t.  -  x  )  =  z+(0)  + 
—  i  — r  — 


V^r 


z+(0i)d0i 


and 


z"(t.-Tf)  =  z'(0)  +  ] 


z  (0.)de. 


(3.39: 


0 


Without  loss  of  generality,  it  is  convenient  to  define 

z_(0)  -  0.  (3.40; 

This  implies  that  z_  (d  )  =  1  and  that 

z+(0)  =  z (0)  (3.41 


is  a  0,1  binary  random  variable.  Substituting  Eqs(3.39)  -  (3.41)  into 


Eq.(3.38)  results  in 


E[Z.  (tp..._Z  (tg^)Z  2+P  '  ‘ 3^—  ,  34-1 )  •  •  (t.,-)] 

^(O^dOj)  ...{z(0)  + 


E[  { z (0)  + 


0  J0 


t,-T 

1  r 


t  n  -T 
>A  r 


M+l  'f 


-+(02i)d':'..i}  J  Q 


tSL2  Tf 


l{di2>d\2Z  (tW+r‘lr)*”i  (t!l3~Tr)-~(t:£3+r'lf)  ' 


...z  (t;-T,)]f  ,  (T  ,i,)di  d.  . 

-  3  f  Tr  If  r  f  r  f 


(3.42 


After  carrying  out  the  products  in  Eq.(3.42),  the  final  result  can 
be  simplified  further  by  interchanging  the  order  of  integration  and 
expectation.  This  latter  step  is  permissable  due  to  the  linearity  of 
the  integration  operation. 

Example  3.8  Evaluate 


E[Z+(t1)2+(t2)Z  (t3)Z+(t4)Z  (t5)] 


t  .  -T 
1  r 


t2  1  r 


t3_lf 


E[  {_z  (0)  +  j  _z+(0^)dfij}{z(0)  +  + 


0  0 


0 


j  l  ■> 2)d-2 }  f  i  <-y d  3 

0  '  0 


_z  (VTr)z  (t5-if)]f..r  Tf  (  r,lf)d;rdTf. 


(3.43 


Carrying  out  the  products  in  Eq.(3.43)  and  using  the  idempotency  property 
of  the  binary  variable  £(()) (  i  . e . ,^(0)^(0)  =  _z(0))yields 


E(Z+(t1)Z+(t2)z‘(t3)Z+(t4)r(t5)l 


u 


k 


b 


t3_Tf 


E[z(0)  i'(63)d63  z  (t4-xr)z  (t5'Xf)]f  lr  Tf(Tr*Tf)dTrdTf 


0  0 


E  [  z (0) 


t0~T 
2  r 

i+(02)d02 


C3  Tf 


A  (03)  (t5_If)]fxr  if  (Tr,Tf)dTrdTf 


0  0 


E  [  z (0) 


crTr 

i+(Q1)d61 


VTf 


i"(03)d03  i+(VTr)i"(t5'Tf)]fxr  xf(VTf)dTr  dTf 


t «  “  T  t  * “  T  ^  o  ^-r 

oooo  lr  z  r  J  i 


0  J0 


E[ 


A+<Vd0i 


Jo 


l+(62)d62 


z“(63)d63  z  (t4-xr)z  (t5-xf)] 


f  C(T  ,1  )dx  di  . 
irif  r  f  r  f 


(3.44) 


Interchanging  the  order  of  expectations  and  integrations  results  in 


E[Z+(t1)Z+(t2)Z~(t3)z'T(t4)Z  (t5)]  = 


co  o° 


0  Jo 


E[z_(0)i  (03)i+(t4-xr)£  (t5-xf ) ] f Tr  Tf(Tr>Tf)d  3dTrdTf 


t0-L  t„-t' 
co  co  2  r  3  f 


E[z(O)i+(02)i  (03)i+(t4-xr)z  (t5-Tf)]fTr.f(ir,:f)d'3d  2dtrd:. 


OO  pOO 


t , — 1  t  -l, 
1  r  3  f 


J0J0J0  ;0 

t,-|  X.,- T  t,  -l, 


ElzCO^C-pi  (n3)i+(tA-xr)i  (t5-  i);i  3d  ldV  f 


1  r  ,2  r 


r  f  f 


3  f 


0^0  4)  J0 


0 


Elh",)A  c  U4"C^  lt5"  f11 


/'/ Y'V'v 


(  i  .is) 
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Note  in  Eq. (3.45  )  that  only  binary  variables  (i.e.  £(0)  and  derivatives 

•  •  __ 

of  counting  processes  (i.e.  £  (t)  and  z  (t))  occur  in  the  higher  order 

joint  moments  of  the  integrands.  V 

Evaluation  of  higher  order  joint  moments  in  terms  of  the  input  switching 

variables,  XI (t) , . . . ,Xn(t) ,  is  considered  next. 

A  seven  step  procedure  to  be  used  in  determining  analytical  expressions 
*'4~  •« 

for £  (t)  and- 2^  (t)  in  terns  of  the  input  switching  variables  and  the  time 

derivatives  of  their  counting  signals  is  discussed  in  Sec. 2. 3.  This 

•  •  — 

procedure  is  applicable  to  finding  expressions  for  jz  (•)  and  jz  (•)  such 
as  occur  in  Eq.(3.45).  However,  now  z  (•)  and  £  (•)  are  evaluated  with 
various  arguments  (e.g.,  0^  Q^,  t^-T^.,  In  order  to  perform 

arithmetic  operations  like  addition  and  multiplication  on  terms  involv¬ 
ing  jz(0),  z^(O)  must  be  expressed  in  the  form  of  an  arithmetic  expression. 
This  can  be  done  by  first  obtaining  the  switching  expression  for  z(0) 
in  terms  of  the  input  switching  variables,  XI (t) , . . . ,Xn(t) .  Next  the 
switching  expression  is  modified  into  the  corresponding  arithmetic  ex¬ 
pression  utilizing  the  arithmetic  expression  theorem  of  Ch.2.  The 
output  higher  order  joint  moments  can  then  be  obtained,  as  illustrated 
in  the  following  example. 

Example  3.9  Consider  the  physical  NAND  gate  whose  model  is  shown  in 
Fig.  3.11.  The  inputs  Xl_(t)  and  X2(t)  are  assumed  to  be  statistically 
independent.  It  is  desired  to  evaluate  the  output  second  order  joint 
moment  =  E[^(t1>Z+(t2)  ] . 


Fig.  3.11.  Circuit  model  for  a  physical  NAND  gate. 


The  first  step  is  to  express  the  desired  function,  R^j+Ct^.t^,  in  terms 

of  the  delay  element  input  signal  z(t)  and  the  time  derivatives  (t)  and 

•  ^  _ 

£  (t).  Decomposition  of  Z(t)  into  £  (t)-Z  (t)  yields 


RZZ+(tl’t2)  =  E[{Z+(t1)-Z"(t1)}Z+(t2)] 


E[Z+(t1)Z+(t2)]  -  E[z‘(t1)Z+(t2)] 


Transforming  in  Eq.(3.46)  the  counting  processes  and  their  derivatives 
in  terms  of  those  at  the  delay  element  input  results  in 


RZZ"Ktl,t2^  “  E[£+(t1-Tr)z+(t2-Tr)]  ~  EtzfUj-  If)z+(t2~  T^)] .  (3 


Using  Eq.(3.39)  -  (3.41),  z  (t.-  t  )  and  z  (t,-  T <.)  are  expressed  in 
integral  form.  Substitution  into  Eq.  (3.47)  results  in 

R„  •+(!., t.)  =  E[(z(0)+  z+(8)d0}z*(t,-  T  )]-E[  z"(0)d9  z+(t,-  x  )] 

—  —  1  2  -  JQ-  “  2  -r  Jo  -2-r 


fl"  % 


E[z(0)z 


i — t*  r  i  —1 

z+(t„-  T  )]  +  E[  z+(0)d0z+(t„-  x  )]-E[  r(0)d0z+(t 


Utilizing  the  conditional  expectation  theorem  [18]  yields 


E[z(0)zT(t2-Tr)]ftr(xr)dTr 


t.  -x 
1  r 


E[ 


z+(0)d0  ‘z+(t„-  x  )]f  (x  )dx 
—  —  2  r  xr  r  r 


00  00  Tf 


E[  £  i"(0)d0z+(t2-Xr)]fTr  Tf(Tr»Tf)dTrdTf-  (3-49) 


Because  integration  is  a  linear  operation,  the  order  of  integration  and 
expectation  can  be  interchanged.  Therefore, 


00  t,-x 
1  r, 


Rzz+(ti*t2)  "jj  Eiz(0)i+(t2-xr)]fTr(xr)dxr  +  jojQ  E[|f(e)z+(t2-xr)]fTr(xr)dedx 


oc  oo  tl-  Cf 


•0 


0J0 


EU"(0)_Z  (t2-Xr)]fTrTf(xr,xf)d0dxrdxf . 


(3.50) 


Note  that  the  integrands  in  Eq.  (3.50)  consist  only  of  the  binary  signal 

•f  •+  •  — 

z^(O)  ■  z^  (0)  and  the  derivatives  of  the  counting  processes  z  (t)  and  z  ^t). 

The  next  step  is  to  express  £+(tj,t2)  terms  t'1e  *-nPut 

switching  signals,  5U(t)  and  X2(t) ,  and  the  derivatives  of  their  counting 
signals. 

The  NAND  switching  operation  is  given  by 


Rzz+(tl,t2)  =  E[{l-Xl_(0)X2(0)}{2y"(t  -T  )X2(t„-x  )+X2~(t  -x  )Xl(t„-T  )}]f  (t  )dT 

—  JQ  1  ^  r  t  x  lx  Trr 


00  t.-x 
r  r  1  r 


’  f  1  r- 

E[{Xl_(0)X2(0)+X2"(0)n(0)}{Xl"(to-T  )X2(t0-x  ) 

Jo  Jo  z  r  —  1  r 


+  xf(t2-Tr)^(t2_Tr)}]fTr  <Tr)d0dTr 


t .  -T  _ 

GO  CO  \  f 


-iff  E[{iL+(0)X2(0)  +  X2+(0)ja(0)}{XI_(t,-T  )X2(t„-X 
JqJ0Jo  i.  r  t 


+  X2_(t0-x  )X1  (t  -x  )}]f  - (x  ,xjd9dx  dx£  . 

—  2  r  —  2  r  xrxf  r  f  r  f 


(3.54) 


Using  the  distributive  law  and  the  statistical  independence  of  the  inputs 
Xl (t)  and  X2(t)  results  in 


RZZ+(tI,t2)  =  j0  <E[X1  (t2-Xr)]EtX2(t2-Xr)]  +  E[X2  (t^)  ]E[3Q  (t^)  ] 


'"RXlxr(0,t2~Tr)RX2X2(0,t2~Tr)~RXlXl  (0,t2~Tr)RX2X2~(0,t2  "Tr)}fTr(VdTr 


■°°  ft  i 

+  .  0.  o  {Rxrxrf0,t2~Tr)RX2X2(9,t2~Tr)'1'  Rxm~(t2~Tr,6)RX2Xf(9,t2~V 


+  +  Rxlxl  (e,t2-Tr)RX2_X2-^e,t2_Tr^f xr(Tr)d('dTt 


00  00  t,-Xr 
f  r  r  1  f 


[  [  fjRXl+Xl-<e't2-Tr>RX2X2^*t2-Tr>  +  ^IXl*  (t2_Tr,0)RX2Xf^,t:2~V 


+  RXlxr(e,t2"Tr)RX2X2+(t2"'rr,0) 


+  *X1X1 <0  ’  VTr)RX2+X2-  (9  * W  }f  trxf  (Tr  *Tf  >d0dTrdTf  ■ 


(3.55) 


Eq. (2.81)  is  applied  to  Eq.(3.55)  to  interchange  order  of  differentiation 

•  __  •  _ 

and  expectation,  i.e.  E[Xk  (6)]  =  E[Xk  (9)].  The  same  argument  can  be 
used  in  the  correlation  functions  involving  derivatives.  For  example. 


Rxkik+(t1,t2)  -  EiHcUpa&t,,)]-  ■—  E[Xjt(t1)Xk+(t2)] 


3t2  tRXkXk+(tl,t:2^ 


(3.56) 


The  interchange  can  also  be  made  when  the  correlation  function  involves 
two  derivatives.  For  example. 


RXk+Xk“(tl*t2)  "  E[Mc+(t1)^~(t2)]  -  -£~  E[Xk+(t1)Xk  (t2)] 


9  r  9  „r.„  + 


j~  [--  E[)&;(t1)Xk  (t2)]]  =  -3tp-t^  E[Xk  (t1)3Ck  (t2)] 


3t13t2  [ISac+xk-(tl,t2)] 


(3.57) 


Substituting  Eqs.(2.81),  (3.56),  and  (3.57)  into  Eq.(3.55),  one 
finally  obtains 


RZZ+(tl,t2)  “  J0{^1—  (t2-VlE[^(t2_Tr)]  +  E[—  (t2-Tr)]El^-(t2*T 


-x  )] 
r 


Btz  ^RXlXl~(0,t2~'rr^RX2X2^0,t2’rTr)~RXlXl(0,t2  V  3t2  [RX2X2~(0>t2~Tr) ^ 


f  (x  )dx 
xr  r  r 


ntl-Tr  a2  a 

+  J0  J0  *  303t2  ^RXl~Xl~(9,t2~Tr^RX2X2(9,t2~Tr^  +  90  ^^jia"(t2'Xr,0)  ] 


dt2  [RX2X2~(9,t2  Tr^ 


+  3t2  [RX1X1~(9,C2  V]30  tRX2X2~(t2~V  9) 


+  R3an(0,t2"Tr)  303t2  [RX2  X£  (0,t2  V]  }fxr (Tr)d9dTr 


rrrtl_T5  .a2 


l0J0J0  {  303t2  (e,t2~Tr)]RX2X2(9,t2~Tr) 


+  96[RXlXl4~(G2'Tr,9)]  3t2  ^RX2  X2“(9,t2  Tr^+  3t2  tRXlxr(e,t2~Tr^ 


30[RX2X2+  (t2_Tr,9^+RXlXl(9,t2"Xr)  303t„[RX2+X2-(0,t2  Tr)1}fmf(Tr’VdedT 


In  the  evaluation  of  R^+ft^  §t2^  ic  is  r^<lu±C'ed  to  know  all 
the  expected  values  and  correlation  functions  in  Eq.(3.58).  However* 
all  correlation  functions  of  the  signal  Xk(t)  and  its  counting  signals 
Xk+ (t)  and  Xk  (t)  can  be  expressed  in  terms  of  three  basic  correlation 
functions.  A  procedure  for  determining  any  correlation  function  in 
terms  of  the  basic  correlation  functions  is  discussed  next.  V 

There  are  10  distinct  correlation  functions  which  may  be  defined  for 
the  signal  Y_(t) »  its  complement  Y’(t)  and  its  counting  signals  ]f+(t) 
and  Y_  (t) .  These  correlation  function  are  listed  below. 

®YY ^1 » ^2^  *  ^YY 1  ^1’  t2^  ’  * C2 ^  1  * ^”2^  *  * Y  *  ^1  * ^2^  * 


*®y»y~^1  *^2^  * Ry+y+  ^ i  *®y"*y—  ^1  *^”2^  *  ^Y— Y—^l  ’^2  ^ 

However,  using  the  decomposition  Y(t)  =  Y+(t)  -  Y  (t) ,  the 
first  seven  correlation  functions  on  the  above  list  may  be  expressed 
in  terms  of  Ry+Y+(t1,t2)  ,RY+y-(t1  ,t2)  ,  Ry-y-^  ,t2> ,  E[Y+(t)], 
and  E[Y  (t)].  The  expected  values  E[lf+(t)]  and  E[Tf  (t)]  arise  in 
correlation  functions  involving  the  complement  Y?(t)  =  1  -  Y(t)  . 

The  following  example  illustrates  the  procedure. 

Example  3.10.  Consider  the  correlation  function  Ry ,y+(t^ ,t2) .  Note 
that  Y/  (t)  =  l-]f(t)  =  1  -  Y_+(t)  +  Y  (t)  .  Therefore, 

Ry,y+(t1,t2)  =  E[Y’(ti)Y+(t2)]  =  E[{l-Y+(t1)+Y"(t1)}Y+(t2)]  v  .. 


Using  the  distributive  law 


Ry.Y+O^tj)  =  E[Y+(t2)]  -  E[Y+(t1)Y+(t2)]  +  E[Y"(tl)Y+(t2) ] .  (3.60) 

By  the  definition  of  the  correlation  function  (Eq.(3.5)) 

Vy+  <t1»t2)  “  E(Y+(t2)3  -  Ry+y+  (t1,t2)  +  Ry+Y-  (tj.tj).  (3.61) 

Note  that  the  right  side  of  Eq.(3.61)  contains  the  expected  value 
E[_Y+(t2)]  and  two  of  the  three  basic  correlation  functions.  Expressions 
similar  to  Eq.(3.61)  may  be  obtained  for  the  other  six  correlation  func¬ 
tions.  V 

Observe  that  Eq.(3.61)  consists  of  three  terms  on  the  right  hand 
side.  In  general,  expressing  correlation  functions  in  terms  of 

the  basic  three,  RY+Y+(ti  *C2>  *  VV^l  **2^  ’  ^“y'^I  ,t:2^  *  and  the  exPected 
values  of  the  counting  signals,  E[Y+(t)]  and  E[Y^  (t)),  result  in  an 
increase  in  the  number  of  terms  involved  in  expressions  similar  to  Eq.(3.58) 
For  example,  Eq.(3.58)  consists  of  the  sum  of  12  terms.  Expressing  the 
various  correlation  functions  in  terms  of  the  three  basic  ones  and  the 
expected  values  yields  an  expression  consisting  of  the  sum  of  52  terms. 
Obviously,  for  efficiency  of  computation,  summation  in  the  intergrands 
should  be  performed  before  integration.  For  example,  if  this  is  done  in 
Eq.(3.58)  evaluation  of  Rz^+(t^,t2)  requires  one  single  integration, 
one  double  integraion,  and  one  triple  integration. 

It  is  next  shown,  by  way  of  a  simple  example,  how  the  basic 


By  definition 


RX+X+(tl’t2)  =  E[X+(t1)X+(t2)].  (3.62) 

Recall  that  3(+(t)  is  an  integer  valued  signal.  It  follows  that 
X+(tl)-X+(t2)  is  also  integer  valued.  Therefore, 

OO 

E[X+(t1)X+(t2)]  =  l  k  Pr{X+(t1)X+(t2)  =  k}  .  (3.63) 

k=-°° 

In  this  case,  )(+  (t)  takes  on  only  the  values  0  and  1  .  Hence  the 
summation  in  Eq.(3.63)  is  reduced  to 

E[X+(t1)X+(t2)]  =  0-Pr{X+(t1)X+(t2)  =  0}  +  l-Pr{x+(tl)X+(t2)  =  1} 

=  Pr{X+(t1)X+(t2)  =  U.  (3.64) 

Note  that  the  product  X+(-.j)X+(t2)  takes  on  the  value  1  only  if  both 
JC  (t^)  =  1  and  =  *  ’  Thus  Eq.(3.64)  may  be  written  as 

E[X+(t1)X+(t2)]  =  Pr{X+(t1)  =  1,  X+(t2)  =  1}.  (3.65) 

Observing  Fig.  3.13(a),  one  can  see  that  X+(t)  =  1  if  and  only 
if  u^  <  t.  It  follows  that 
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Combining  Eqs.(3.62)  and  (3.66)  results  in 


RX+X+^tl,t2^ 


The  results  in  Eq.(3.67)  are  sketched  in  Fig.  3.14. 


Fig.  3.14  The  autocorrelation  Rv+v+(t. ,t„) 


The  methods  discussed  so  far  for  a  single  logic  block  can  be 
utilized  in  analyzing  combinational  networks  containing  reconvergent 
fanout.  As  pointed  out  in  Ch.  1,  the  output  expected  value  E[Z(t)] 
can  be  used  in  evaluation  of  system  performance.  It  was  observed  in 
Sec.  3.1  that  higher  order  moments  may  arise  when  output  expected  values 
of  networks  containing  reconvergent  fanouts  are  considered.  It  was 
also  shown  and  illustrated  in  this  section  how  output  higher  order 
joint  moments  of  a  logic  block  can  be  obtained  in  terms  of  input  higher 
order  moments  and  the  p.d.f's  of  the  delays.  An  efficient  method  for 
evaluating  the  output  expectation,  E[_Z(t)],  of  a  general  combinational 
network  is  described  below. 

In  networks  containing  reconvergent  fanouts  it  is  possible  that 
entire  tree-like  subnetworks  may  exist.  Entire  tree-like  subnetworks  are 
defined  to  be  those  in  which  there  is  only  one  path  from  any  primary 
input  Xk  to  an  output  Zj[  from  the  last  network  logic  level.  Entire 
tree-like  subnetwork  output  expectations  are  evaluated  by  progressing 
from  the  primary  inputs  to  the  outputs,  logic  block  by  logic  block, 
as  described  in  Sec.  2.3.  For  the  subnetworks  including  reconvergent 
fanouts,  the  output  expectations  are  determined  in  terms  of  the 
higher  order  joint  moments  of  the  counting  signals  associated  with  the 
reconverging  signals.  However,  only  some  joint  moments  of  a  given 
order  may  be  necessary.  In  order  to  avoid  unnecessary  computations 
and  storage  of  higher  order  joint  moments,  it  is  desirable  to  first 
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determine  those  moments  which  will  be  needed.  This  is  done  by  1) 

•  +  •  w 

deriving  the  analytical  expressions  for  E[£-(t)]  and  E[£  (t)]  in  terms 
of  the  reconverging  signals  and  2)  expressing  the  resulting  higher 
order  moments  of  the  reconverging  signals  in  terms  of  those  of  the 
network  inputs  XI (t) , . . . ,Xn(t)  and  their  counting  signals.  The 
procedure  is  illustrated  in  the  following  example. 

Example  3.12.  Consider  the  combinational  network  shown  in  Pig.  3.17? 
All  gates  are  assumed  to  contain  discriminating  delay  elements.  Also, 
the  inputs  Xl_(t)  and  X2(t)  are  assumed  to  be  statistically  Independent. 
Since  there  are  two  paths  to  the  output  for  X2(t) ,  the  network  contains 
a  reconvergent  fanout.  It  is  desired  to  evaluate  the  expected  value 
of  the  output  E  £_Z  ( t)  1 . 


Z(t) 


Fig.  3.17.  Combinational  network  with  a  reconvergent  fanout. 

The  model  for  the  physical  network  in  Fig.  3.17  is  illustrated 

in  Fig.  3.18.  Observe  the  ideal  logic  output  signals  £(t)  and  £(t) 

and  the  discriminating  delay  elements  with  the  random  delays 

and  Tc..  Also,  note  that  the  network  does  not  contain  an  entire  tree-like 
— f  A 


subnetwork. 


Fig.  3.18.  Model  for  network  in  Fig.  3.17. 

Because  of  the  discriminating  delays,  the  derivation  for  E [ Z^(t) ] 
requires  that  EfZ+Ct)]  and  K[Z  (t)]  be  determined  separately.  Then, 
by  applying  the  expectation  toEq.(2.45),  E[j2(t)]  is  given  by 

E[Z(t)]  =  E[Z+(t)]  -E[z"(t)].  (3. 

Observe  that  the  signal  X2(t)  splits  and  then  reconverges  at  the  node  z^. 

To  determine  which  higher  order  joint  moments  of  X2(t)  are  needed,  the 

•  •  — 

arithmetic  expressions  for  (t)  and  (t)  in  terms  of  the  reconverging 
signal  X2(t)  are  obtained  first.  Note  that  £(t)  is  the  output  of  the 
ideal  logic  circuit  performing  the  AND  operation.  Therefore, 

z(t)  =  X2(t)-Y(t).  (3. 

In  Eq.  (3.73)  the  right  side  can  be  considered  to  be  in  the  form  of 
a  sum  of  nonoverlapping  implicants.  Differentiation  of  the  corresponding 
arithmetic  expression  with  respect  to  time  yields 


£(t)  -  X2  (t)Y(t)  +  Y(t)X2(t).  (3.74) 

In  order  to  place  Eq.(3.74)  in  the  form  of  Eq.(2.63),  define 
Wl(t)  -  X2(t) 

W2(t)  -  Y(t) 

Then  Eq.(3.74)  can  be  written  as 
2 

l<t)  =  l  Wk(t)Bk(t) , 
k=l 

where 

B1 (t)  -  Y(t) 

B2(t)  -  X£(t). 

Note  that  Bl^(t)  and  IJ2(t)  are  in  the  arithmetic  canonical  form  required 
by  Eq.(2.72).  With  reference  to  Eq.(2.74),  it  follows  that 

Bl+(t)  -  Y(t),  Bl,_(t)  -  0,  B2+(t)  -  X2(t),  B2~(t)  «  0.  (3.75) 

Using  Eq.(2.79),  one  obtains 

£+(t)  -  WJ*(t)Y(t)  +  W2+(t)X2(t)  -  Xgtt)Y(t)  +  Y+(t)X2(t) 

i'(t)  “  W~(t)Y(t)+  W2_(t)X2(t)  -  X2"(t)Y(t)  +  Y~(t)X2(t) .  (3.76) 


Recognizing  that  ]f(t)  ■  Y_  (t)  -  Y  (t)  ,  Eq.(3.76)  become 
l+(t)  =  X2+(tXl+(t)  -Y"(t)}+  Y_+  (t)X2(t) 

£"(t)  =  X2“(t)(Y+(t)  - Y~(t)}  +  Y-(t)X2(t).  (3.77) 


Each  of  the  equations  in  Eqs.(3.77)  are  now  developed  separately. 


The  derivative  of  the  rise  counting  signal,  z  (t) ,  is  considered 


first.  Note  that  Y+(t)  =  j£+(t-  T^)  and  Y_  (t)  =  Z  (t-T^,).  Consequently 


z+(t) 


x2+(t)  (t-  t^)  -  £  (t-xf0) }  + 


(3.78) 


The  counting  signals  Z  (t“ir0)  and  Z  are  exPressed  in  terms 

of  their  derivatives  by  means  of  Eqs.(3.39)  -  (3.41).  Hence, 


t-x  n 
— ru 


t_— f  0 


,  — ru  (■  -tu 

£  (t)  -  X2+(t)  (Z(0)  +  i  (0)d0  -  i  (0)de}+f+(t-x  )X2(t) 

0  Jo 


t-  T 


=  X2tt)y(0)  +  X2+(t) 


Z+(6)d0  -  X2+(t) 


t"ifO 


y  (0)d0  +ZT(t-Tr0)X2(t) 


=  X2(t)y(0)  ^ 


t--r0 


X2+(t)y+(0)dG  - 


"“IfO 


X2+(t)x  (9)d0  +  fV-T^X^t) 


With  respect  to  ^(t)  only  2,(0) ,  _y(t),  and  y  (t)  are  involved  in 
Eq.(3.79).  These  signals  can  now  be  written  in  terms  of  the  primary 
input  signals.  In  particular,  2(0  =  5C1^C t >  v  X2 (t)  ,  and  a  corresponding 
arithmetic  expression  is  2(0  =  5U(t)  +  XI '  (t)X2(t) .  Therefore, 

2(0)  =  XI (0)  +  XI ’ (0)X2 (0)  (3.80) 

and 

£(t)  =  XI (t)  +Xl'(t)X2(t)  +  Xl'(t)X2(t).  (3.81) 

Using  Eq. (2.62)  in  Eq.  (3.81)  yields 

y(0  =  Xl_(t)  -  XHt)X2(t)  +X2(t)ja'(t) 

-  XI (t)  (l  -X2(t)}  +X2(t)Xl’(t) 

-  Xl(t)X2’(t)  +  X2(t)Xl/(t).  (3.82) 

Developing  Eq.(3.82)  into  the  form  of  Eq.(2.78),  it  follows  that 

Bl+(t)  -  X2'(t),  Bl"(t)  •"  0,  B2+(t)  -  XI.' (t),  B2*(t)  -  0.  (3.83) 

The  use  of  Eqs.(3.83)  in  Eqs.(2.79)  results  in 


Substitution  of  Eqs.(3.80)  and  (3.84)  into  Eq.  (3.79)  yields 


:t_^o 


i+(t)  -  x2+(t)xi (0)  +  x2+(t)2y.' 1  (0)ic2 (o)  +  I  x2+(t)xf(e)x2’(e)  de 

jo 

+  X2  (t)X2  (e)ja' (0)d6 

Jo 


f^o 


-cu.  +  •  -  f  ^u*  +  •  - 

X2  (t)Xl  (6)X2’(d)de  -  X2  (t)X2  (8)5a’ 

0  Jo 


(0)d8 


+  XI  (t-Tr0)X2’ (t-Tr0)X2(t) 


(3.85) 


The  right  side  of  Eq.(3.85)  consists  only  of  primary  input  signals. 

The  expected  value  of  (t)  is  now  readily  obtained.  Making  use 

of  the  statistical  independence  of  the  inputs  XI (t)  and  X2(t) ,  Et£<(t)] 

becomes 

t_ir0 

E[i+(t)]  =  E[X2+(t)]E[Xl(0)]  +  E[Xl_,(0)]Rx2'2+(0,t)+  E[Xl+(0)  ]Ry2  ,£2+(0  ,t)d6 

-  Jo  - 

t-T  n  t— T  rn 

■  — r0  r  — f0, 

+  E[Xl'(0)]R^2+^2+  (0,t)d0  -  j  E[X1“(8)]R^2+X2.  Ct,8)d0 


t“lfo 

E[2a'(0)]R^2+.2-(t,0)de 


+  E[ja  (t-  Ir0)]Ry2X2,(t,t_-r0)  +  EIXL'^Itq)!^  x2+^ ,t--^r0^  ’  (3.86) 


E 


Next  the  order  of  differentiation  and  expectation  is  interchanged  and 
the  conditional  expectation  theorem  is  used  to  obtain 


E[z+(t)]  -  E[xfl‘(t)]E[Xl(0)]  +  E[n*  (0)]  ^  [RX2X2+r0,t)] 


fooft~Tr0.  .  9 

+  E[X1>)]^  [Rx2,x2+(e,t)]flrO(TrO)d0dTrO 

0  0 

00  t-T  2 

+  f  |  'Elxi'ce)]-^  [VVKe.tOlf^  (Tr0>d6dTr0 

Jo  Jo  -  — 


00  t-T, 


f  f  .  -  a 

j  (0)1~9^~  tRx2'x2+^e,t^fTfO^TfO^d0dTfO 


°°  t-T, 


ru 

E[X1  '  (0)  ]  g£gQ  t  ^ C  ’  9 ^  ^  f  T f  0  ^ Tf  0 ^  d ^f  0 

J0  0  -  — 

00 

+  (t~TrO)]RX2X2*  (t,t~Tr0^fTr0(Tr0^dTr0 


+  J0El— ' (t~TrO)]  80  [RX2X2+  (t,0_TrO)1  0_t  f  T_rO(TrO)dTrO‘ 


(3.87) 


Finally,  the  derivative  of  the  expected  value  of  the  output  rise  counting 
signal,  E[Z_+(t)],  is  determined  using  Eq.(2.90).  Specifically, 


E[Z+(t)]  =  E[z+(t)]*fTrA(t). 


A  similar  derivation  for  E[z  (t)]  yields 


(3.88) 
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E[.z‘<t)]  -  E[X2~(t)  ]E[X1  (0)  ]+E[Xl '  (0)]  [^"(O.t)  ] 

t_Tr04-  » 

+  JJo  EtxiJ +(0)  £  (e.t)]flr0(rr0)dedrr0 


CO  t-T 


■  ru 


°°ft-Tf0 


E[X1_  (0)]  [Ex2tX2  (9*t)]  fTfO(TfO)d  Tf0 


f  f0  a2 

e[xi '  (0) ]  ^  iR^-xz^'^WW^fo 
Jo  Jo  —  —  ~ 

oo 

+  ^  EtXl^  (t-Tf0)  ]Rx2x2,(t,t-Tf0)fTf0  (xf0)dTf0 


r°°  a 

E[xr(t-Tf0)]  £  [R^  x2-  Ct  e-Tf0)]  /xf0^f0)dTf0- 

J0  —  -  0=t  ~ 


(3.89) 


Once  again,  using  Eq.  (2.90) 

E[Z"(t)]  =  E[z_(t)]*fTfA(t) . 


(3.90) 


Finally,  substituting  Eqs.(3.88)  and  (3.90)  into  Eq.(2.88)  and 
performing  the  integration  with  respect  to  time  yield 


E[  Z  (t)  ]  =  E[Z(0)  ]  +  E[Z  (0)]d0  -  E  [  Z  ~  ( 0 )]  d  0 . 


(3.91) 
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Because  of  Eqs.(3.87)  and  (3.89),  note  that  E[£(t)]  depends  on  the 
first  and  second  order  moments  involving  the  reconverging  signal 
X2(t)  and  its  counting  signals.  However,  only  the  expected  value  of 
JQ(t)  and  its  counting  signals  are  required. 

In  summary,  the  procedure  for  evaluating  the  output  expected 
value  of  a  combinational  network  with  discriminating  delay  elements 
includes  the  following  steps.  First,  the  network  is  analyzed  to 
identify  reconvergent  fanouts.  The  network  is  then  separated  into 
entire  tree-like  subnetworks  and  subnetworks  consisting  of  the 
reconvergent  fanouts.  Entire  tree-like  subnetworks  are  defined 
to  be  those  in  which  there  is  only  one  path  from  any  primary  input 
Xk  to  an  output  Zj_  from  the  last  network  logic  level.  Entire  tree¬ 
like  subnetwork  output  expectations  are  evaluated  by  progressing  from 
the  primary  inputs  to  the  outputs,  logic  block  by  logic  block,  as 
described  in  Sec.  2.3.  For  the  subnetworks  including  reconvergent 
fanouts,  the  output  expectations  are  first  determined  in  terms  of 
the  higher  order  joint  moments  of  the  counting  signals  associated  with 
the  reconverging  signals.  This  step  was  illustrated  in  Ex.  3.12. 

Next  the  higher  order  joint  moments  obtained  for  the  reconverging 
signals  are  evaluated  in  terms  of  the  higher  order  joint  moments  of 
the  counting  signals  associated  with  the  primary  inputs.  This  step 
was  illustrated  in  Ex.  3.9.  Finally,  the  results  of  the  last  step 
are  substituted  into  the  expression  for  the  output  expected  value, 
to  yield  an  expression  depending  only  on  expectations  and  higher 
order  joint  moments  of  signals  associated  with  primary  inputs. 


3.4  Approximations  for  Higher  Order  Moments  of  0,1  Binary  Processes 


In  their  very  significant  contributions  [15,16],  Bass  and  Grundmann 
employed  an  interesting  and  easy  to  use  approximation  for  the  auto¬ 
correlation  of  a  0, 1  binary  process  without  discussing  the  conditions 
under  which  their  approximation  is  valid.  This  approximation  is  devel¬ 
oped  in  this  section.  In  addition  to  discussing  the  limitations,  it 
is  also  shown  how  the  approximation  can  be  generalized  to  higher  order 
moments . 


Consider  a  0,1  binary  stochastic  signal  X(t)  whose  autocorrelation 
function,  R^Ct^.t^),  is  required.  In  reference  to  Fig.  1.5,  suppose 
every  sample  of  X(t)  has  a  minimum  gap  duration,  denoted  by  u  ,  such 
that  i)  _<_  'j  ;  i=l,2,,..  Also  suppose  every  sample  of  X(t)  has  a 
minimum  pulse  duration,  denoted  by  6^,  such  that  6^  ;  j  *  1,2,... 

Let  t^  =  t  and  =  t  +  A,  where  A  is  to  be  considered  as  a  fixed 
constant  during  the  discussion.  A  new  0,1  binary  stochastic  process 


is  defined  to  be 


Y^t)  =  X(t)X(t+A) 


(3.92) 


The  autocorrelation  function  of  X(t)  is  simply  the  expectation  of 
Y^(t) .  Specifically, 


R^t.t+A)  =  E[X(t)X(t+A)]  =  E[YA(t> } 


(3.93) 
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Consider  a  sample  function 


YA(t)  =  X(t)X(t+A) . 


A  rise  transition  in  Y^(t)  exists  only  when  a  rise  transition 
occurs  in  one  of  the  factors  of  Y^(t)  while  the  other  factor  takes 


.th 


on  the  value  1.  This  is  illustrated  in  Fig.  3.19.  If  |A|  <  u^,  the  i 
rise  transition  of  Y^(t)  occurs  at  the  i^  rise  transition  of  X(t)  when 
A>0  and  at  the  i*"*1  rise  transition  of  X(t+A)  when  A<0.  In  effect,  the 
i^  rise  transition  in  YA(t)  occurs  at  the  i*"*1  rise  transition  of  either 
X(t)  or  X(t+A),  depending  on  which  comes  last.  It  follows,  for  every 
sample  function  Y.(t),  that  the  rise  counting  signal  is  given  by 


^A (t)  = 


X  (t),  A>0 


1- 


=  min{X+(t),  X+(t+A)}. 


(t+A),  A<0 


(3.94) 


Eq.  (3. 

94)  is 

illustrated  in 

Fig.  3.19  for  0 
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Fig.  3.19.  Sample  functions  of  a)  21(0,  b)  21(t+A),  c)  Y.(t) 

f  |  ““A 

d)  X  (i)  ,  e)_X  (t+-A)  and  f)  _Ya  (t)  with  the  restriction  that  0<  A 
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It  is  Important  to  emphasize  that  Eq.(3.94)  does  not  hold  when 
| A |  >  Uq.  This  is  illustrated  in  Fig.  3.20. 
t  X(t) 


Fig.  3.20.  Sample  functions  of  a)  X(t),  b)  2i(t+A),  C)Y^(C)  under  the 
condition  that  u^<  A. 

Even  though  A>0,  the  second  rise  transition  in  (t)  occurs  at  the  second 
rise  transition  of  X(t+A)  as  opposed  to  the  second  rise  transition  of 
X(t) .  This  is  in  contradiction  to  Eq.(3.94).  The  difficulty  arises 
because  A  >  u^. 

Similar  reasoning  holds  for  the  fall  transitions  of  Y^(t).  The 
binary  signal  Y^(t)  has  a  fall  transition  if  one  of  the  factors,  X(t)  or 
X(t+A),  has  a  fall  transition  while  the  other  takes  on  the  value  1. 

This  is  illustrated  in  Fig.  3.21. 

If  j A ]  <  the  fall  transition  of  Y^(t)  occurs  at  the  jth  fall 

transition  of  X(t+A)  when  A  >  0  and  at  the  jth  fall  transition  of  X(t)  when 
A  <  0.  In  effect,  the  jt!'  fall  transition  of  Y^(t)  occurs  at  the  jth 
fall  transition  of  either  X(t)  or  X(t+A),  depending  on  which  comes 
first.  It  follows,  for  every  sample  function  Y^(t),  that  the  fall 


counting  signal  is  given  by 

X'(t  +A)  ,  A  >  0 

y'a  (t)  - 


max  {X  (t),  X  (t+A)} 


(t)  ,  A  <  0 

Eq.  (3.95)  is  illustrated  in  Fig.  3.21  for  0  <  A  <  6_. 


YA(t)=X(t)X(t+A) 


X  (t+A) 


Y.(t)  _  _ 

*'  max-.  X  (t )  ,X  ( t+V) 


Fig.  3.21.  Sample  functions  of  a)  X(t),  b)  X(t-t  A) ,  c)Y.  (t), 


d)  X  (t),  e)  X  (t+A)  and  f)  Y_  .(t)  with  the  restriction  that 
0  <  A  <  Q. 

Eq.  (3.95)  is  not  valid  when  •  .  This  is  illustrated  in  Fig.  3.22 


1  X(t) 

U1  dl 


X(t+A) 


Uj-A  dj-A  u^-A  d 

1  V*>  _ 


2  2 


Fig.  3.22.  Sample  functions  of  a)  X(t),  b)X  (t+A),  c)Y,  (t)  under  the 
condition  that  <  A  . 


Even  though  A  >  0,  a  fall  transition  does  not  occur  in 

Y^(t)  when  the  first  fall  transition  occurs  in  X(t+A) .  Eq.(3.95)  is 

violated  because  A  >  6^. 

Combining  the  piecewise  results  of  Eqs.(3.94)  and  (3.95)  in  accor¬ 
dance  with  Eq . ( 1 . 5 ) ,  an  expression  for  Y^(t)  is  given  by 


YA(f)  =  Ya  (t)-  Y"  (t)  = 


X  (t)  -  X  (t+A)  A>0 


X+(t+A)-X~(t)  ,  A<0  ,  !A|<min{u0,60}.(3.96) 


An  equivalent  expression,  using  the  min(*)  and  max  (•)  operations  from 
Eqs(3.94)  and  (3.95),  is 


YA(t)  =  min(X+(t),  X+(t+A)}  -  max  (X~(t) ,X~ (t+A) } 


| A l  <  min{u0,S0). 


(3.97) 


Observe  that  the  choice  of  the  upper  or  lower  equation  in  Eq.  (3.96) 
is  determined  only  by  the  sign  of  A.  This  choice  holds  for  each  time 
instant  t  and  for  all  sample  functions  of  the  process  Y^t).  It  follows 
that  Eqs.(3.96)  and  (3.97)  hold  for  the  stochastic  process  YA(t) . 
Specifically, 


Y , (t)  = 


X  (t)  -  X~(t+A)  ,  A>0 


X  (t+A)  -  X_(t)  ,  A<0 


=  nin{X+(t) ,X+(t+A) }-  max{X  (t) ,X  (t+A)}, 


| A |  <  min{Up,6Q} 


(3.98) 


\  .N  “•  .*•  * 


where  the  min{*}  operation  is  applied  to  each  of  the  sample  function 
pairs  from  X+(t)  and  X+ (t+A)  to  generate  the  random  process  Y^+  (t) 

and  Y  ^(t)  is  obtained  in  a  similar  manner  using  the  max  { • }  operation 
on  sample  function  pairs  from  X  (t)  and  31  (t+A) .  The  expected  value 
E[Y^(t)],  is  given  by  taking  the  expectation  of  the  first  part  in 
Eq.(3.98)  as 


EtY^t)]  =  < 


E[X  (t)]  -  E[X_(t+A)],  A>0 

E[X+(t+A)]  -  E [ X~ ( t )  ] ,  A<0  ,  |A|<min{u0,60}.  (3.99) 


Note  that  E[3C  (t)]  and  E[X  (t)]  are  nondecreasing  time  functions  since 


their  respective  time  derivatives,  E[3(  (t)]  and  E[X  (t) ] ,  are  nonnegative 


[15,  p. 18]  - 


It  follows  that 

E[X+(t)]  <  E[X+(t+A)  ]  ,  A>0 

E[X+  (t+A)  ]  <  E[X+(t)]  ,  A<0 


(3.100) 


and 

E[X~(t+A)  ]  >E[X~(t)]  ,  A>0 

E[X~  (t)  ]  >  E[X_(t+A)  ]  ,  A<0  . 


(3.101) 


Consequently,  E[Y^(t)]  can  be  written  as 

E[Y^(t) ]  =  min{E[X+(t)l,E[X+(t+A)]}  -  max{E[X-(t) ] ,E[X_(t+A) ] ) , | A | <min{u0,6Q} 

(3.102) 


Finally,  by  referring  to  Eq.(3.92)  and  letting  t  ■  t.  and  t  +  A  =  t„» 


Eq.(3.102)  becomes 


E[X(t1)X(t2)]  =  =  min{E[X+(t1)],E[X+(t2)]} 

-  max{E[x"(t1)],E[X'(t2)]}, 

provided  | t^-t2 |<  min{u0,60} .  (3.103 

The  result  of  Eq.  (3.103)  is  readily  generalized  to  higher  order 
moments.  Consider  the  n-fold  product 

Y (t)  =  X(t)X(t+  A )  ...X(t+A  .)  (3.104 

1  n— 1 


where  X(t)  is  a  0,1  binary  signal  and  A  ,  k  =  1,  ...,  (n-1),  may  be 
either  positive  or  negative  constants.  A  typical  scatter  of  the  points 
t,  t+A.,...,  t+A  is  sketched  in  Fig.  3.23. 


”  •  H .  n-1 

1  1  1  *  * 

*  1  *  *  '  §  ^  + 

t+A. 

1 

t+A1  t  t+A^ 

t+A  t+A, 

n-1  k 

Fig.  3.23.  Scatter 

Denote 

A^  =  max  {A^» 

of  the  points  t,  t+Aj,..., 

0} 

t+An_lt  0  <_  j,  k  £  (n-1) 

=  m*n  ^^k’ 

0}  ;  k  =  1 , . . . , (n-1) . 

(3.105) 

The  duration  of  the- time  interval  over  which  the  points  t,  t^^ , . . .  ,t+AQ 
are  scattered  is  specified  by 

T  ■  A^  -  A  •  (3.106) 

Note  that  Y(t)  equals  1  only  if  all  n  factors,  X(t)  ,X(t+Ap  , . . .  ,X(t+A^_j 
take  on  the  value  1.  Given  a  set  of  time  shifts  {A^},  such  that 

T  < min{uQ,6Q) ,  (3.107) 

the  ith  rise  transition  of  Y(t)  occurs  at  the  ith  rise  transition  of 
either  X(t),  X(t+A1) , . . .X(t+An-1) ,  depending  upon  which  one  occurs 
last.  Similarly,  the  jth  fall  transition  of  Y(t)  occurs  at  the  j 
fall  transition  of  either  X(t),  XU-W^) , . . .  ,X(t4An-1)  ,  depending  upon 
which  one  occurs  first.  It  follows  that 

Y+(t)  =  min{X+(t),X+(t+A1),...,X+(t+An_1)) 

Y_(t)  =  max{X~(t),X-(t-fA1),...,X-(t-tAn_1)},  (3.108) 

provided  T  <  mintu^.S^}. 

Consequently, 

Y(t)  =  Y+(t)  -  Y~(t)  =  min{X+(t),X+(t+A1) . X+(t+Anl» 

-  max{X  (t),X  (t+A^),...,X  (t+A^_^)}  ;  T <  mintu^.f^} .  (3.109) 
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Extending  che  same  reasoning  to  each  sample  function  of  the 


stochastic  process 


Y(t)  =  X(t)X(t+A, )  ...  X(t+A  .) 
—  —  —  1  —  n-x 


(3.110) 


where  JC(t )  is  a  stochastic  0,1  binary  process  and  A^'s  are  constants, 
it  follows  that 

Y  (t)  =  min  {X+(t),X+(t+A. ),..., X+(t+A  ,  )  }-  max(X-(t)  ,X-(t4A,  )  , . . .  ,X-(t-*A  ,)} 

i  n- 1  —  —  1  —  n — l 

(3. 111. 


T  <  min  {Uq,6q}. 


The  min  (*}and  max  { • }  operations  yield  Y+(t)  and  Y_  (t),  respectively, 
as  indicated  in  the  discussion  following  Eq.(3.98).  Because  the  A^'s 
are  constants  and  the  counting  signals  are  nondecreasing  waveforms, 
it  follows  that 


Y+(t)  =  min{X+(t),X+(t+A, ),..., X+(t+A  .)}  -  X+(t+A  ), 
—  —  —  i  n—  i  —  m 


Y  (t)  =  max  {X  (t)  .X^t+Aj)  ,  . . .  .X^t+A^)}  =  iMt+A^  .  (3.112) 

Eq.  (3.111)  becomes 

Y(t)  =  X+(t+Am)  -  X~(t+A^) .  (3.113) 

The  expected  value  of  £(t)  is  then  given  by 

E[Y(t)]  =  E[X+(t+Am)]  -  E[X_(t+AM)].  (3.114) 

Recognizing  that 

E[X+(t+Am)]  £  E[X+(t+Ak)], 

E[X_(t+AM)]  >  E[x"(t+Ak)]  k  =  l,...,(n-l)  (3.115) 


Eq.  (3.114)  can  be  written  as 

E[Y(t)]  -  min{E[X+(t)],E[X+(t+A1)l.-*-.E[^f(t+An_1)i} 

-  max{E[X'(t),E[X"(t+A1)],...,E[X“(t+An_1)]}jT<min  (uQ,60).  (3.116) 


Letting  t  -  t.  ,  t+A,  -  t„ . t  +  A„  ,  -  t  ,  it  follows  that 

i  l  j.  n- 1  n 

T  *  max{ t^}  -  min  {tk>  ;  k-l,...,n  (3.117) 

and 

E[X(t.)...X(t  )]  *  min{E[X+(t.)] . E[X+(t  )]}  -  max{E[X'(t, ) ] . E[X~(t  )]} 

i  n  —  i  —  n  l  —  n 

T  <  min{u0,60}.  (3.118) 


Eq.  (3.118)  is  the  generalization  of  Eq. (3.103)  to  nCl1  order  moments. 

The  result  in  Eq. (3.103)  is  the  approximation  used  by  Grundmann  in 
[15],  The  derivation  and  limitations  of  this  approximation  have  been 
presented  in  this  section.  Unfortunately,  the  restriction  T <  min{uQ,6g} 
cannot  always  be  met.  In  fact,  in  general,  the  existence  of  nonzero 
values  of  and/or  6^  is  not  guaranteed.  In  such  cases  use  of  the  approx¬ 
imation  results  in  some  error.  If  the  error  is  unacceptable,  the 
higher  order  statistics  introduced  in  this  chapter  must  be  used. 


4.  APPROXIMATE  MODELS  FOR  LARGE  LOGIC  BLOCKS 


4.1  Introduction 

As  was  seen  in  the  previous  chapters,  calculations  required 
for  the  determination  of  output  expected  values  may  be  very  long  and 
complicated,  even  for  simple  combinational  circuits.  One  way  for 
reducing  -the  amount  of  computations  and  memory  required  for  evalua¬ 
tion  of  complicated  combinational  circuits  is  to  subdivide  the  circuit 
into  large  logic  blocks.  Each  large  logic  block  is  modeled  as  a  single 
entity  and  the  entire  circuit  is  analyzed  by  taking  into  account  the 
interconnections  of  the  large  logic  blocks.  In  general,  a  large  logic 
block  will  contain  many  gates  involving  several  levels  of  logic.  In 
addition,  the  paths  connecting  the  various  inputs  and  outputs  may  differ 
significantly  in  terms  of  the  number  and  types  of  gates  encountered. 

For  example,  one  input  may  be  tied  directly  to  the  last  gate  of  the 
logic  block  whereas  other  inputs  may  go  through  several  levels  of  logic. 
As  a  result,  insertion  of  delays  into  a  simplified  model  of  a  large 
logic  block  is  extremely  difficult. 

Several  approaches  have  been  suggested  by  Grundmann  [15,  Ch.  5], 

The  most  accurate  appears  to  be  one  in  which  a  different  delay  is 
assigned  at  each  input  for  every  output  affected  by  that  input.  Even 
though  this  approach  does  not  take  into  account  reconvergent  fanouts, 
it  is  still  extremely  complicated  because  a  large  logic  block  with  n 
inputs  and  m  outputs  requires  that  2(n)(m)  delays  be  specified  in  order 


to  account  for  the  delays  associated  with  the  rise  and  fall  transi¬ 
tions  of  all  possible  input-output  pairs.  Another  approach  is  to  sim¬ 
ply  place  a  single  delay  at  each  input  irrespective  of  the  number  of 
outputs  affected  by  that  input.  In  this  way  a  separate  delay  is 
attached  to  each  input.  However,  a  disadvantage  is  that  the  same 
delay  is  attributed  to  every  path  connected  to  a  particular  input  even 
though  the  paths  may  have  different  delays.  A  third  approach  is  to 
assign  delays  both  at  each  input  and  at  each  output.  This  requires  that 
2(n+m)  delays  be  assigned  for  a  large  logic  block  with  n  inputs  and  m 
outputs.  The  assignment  is  complicated  by  the  fact  that  the  delays  are 
likely  to  be  strongly  dependent.  The  last  approach  suggested  by  Grundmann 
is  to  simply  assign  a  delay  to  each  output  of  the  large  logic  block.  A 
disadvantage  of  this  approach  is  that  the  same  delay  is  assigned  to 
every  path  connected  to  a  particular  output  even  though  the  paths  may 
have  different  delays.  Although  Grundmann  proposed  the  various  approaches 
described  above,  he  did  not  discuss  any  strategies  for  characterizing  the 
delays , 

In  this  chapter,  attention  is  devoted  to  the  approach  where  a 
delay  is  assigned  to  each  output  of  a  large  logic  block.  This  model 
is  referred  to  as  the  simple  output  delay  model.  The  justification  for 
this  choice  is  discussed  next.  One  factor  motivating  this  choice  is 
its  relative  simplicity.  As  seen  in  the  examples  of  the  previous  chapters, 
the  complexity  of  the  analysis  increases  drastically  as  the  number  of 
delays  increases.  In  addition,  it  has  been  observed  in  computer  simulations 


[3]  and  in  experimental  investigations  [2,  23]  that  digital  circuits 
are  most  sensitive  to  interfering  signals  injected  into  their  outputs. 

In  general,  interference  in  integrated  circuits  is  picked  up  in  long 
wires  and  conductors.  When  two  logic  blocks  are  connected  by  a  long- 
wire,  the  output  gate  in  the  logic  block  whose  output  is  connected 
to  the  wire  suffers  the  greatest  delays  due  to  interference  on  the  wire 
(see  Fig.  4.1).  Consequently,  the  physical  situation  with  respect  to 


Fig.  4.1.  Layout  of  a  logic  circuit  subjected  to  EMI. 

EMI  is  best  modeled  by  associating  the  delay  with  the  output.  In  addition, 
the  technique  developed  in  chapters  2  and  3  are  readily  applied  to  the 
simple  output  delay  model. 

Several  considerations  arise  when  subdividing  a  complicated  com¬ 
binational  circuit  into  large  logic  blocks.  Because  the  same  delay 
is  assigned  to  every  path  connected  to  a  particular  output,  it  is 
desirable  to  select  large  logic  blocks  such  that  the  number  of  gates 
encountered  in  each  path  from  any  input  to  any  output  is  approximately 


There  is  also  a  limit  to  how  large  a  logic  block  should  be. 

It  is  desirable  to  choose  the  logic  block  to  be  as  large  an  possible 
in  order  to  reduce  the  number  of  blocks  and  simplify  the  ar-alysis. 

However,  the  spreads  in  the  delays  increase  with  the  number  of  logic 
levels.  Once  the  spreads  exceed  the  minimum  time  duration  between 
transitions  in  the  input  signals,  some  pulses  or  gaps  that  actually 
occur  in  the  physical  circuit  are  not  predicted  by  the  ideal  logic  portion 
of  the  model.  In  such  cases,  adding  delay  to  the  output  of  the  ideal 
logic  block  results  in  errors.  The  phenomenon  is  illustrated  in  Example  1.1. 
Because  the  difference  in  the  delays  assigned  to  gates  01  and  02  exceeds 
the  time  duration  between  the  transitions  of  the  input  signals,  an  output 
pulse  is  generated.  If  the  ideal  logic,  without  delay,  is  performed 
first,  the  output  pulse  is  not  generated  and  adding  delay  to  the  output 
of  the  ideal  logic  yields  an  error.  To  avoid  such  errors,  it  is  recom¬ 
mended  that  the  size  of  large  logic  blocks  be  chosen  such  that  the  dif¬ 
ference  between  the  maximum  and  minimum  delays  for  all  paths  within  the 
block  be  less  than  the  minimum  time  duration  between  transitions  in  the 
input  signals , where  durations  from  one  input  to  another  are  considered  as 
well  as  durations  within  each  input.  This  is  a  conservative  recommendation 
in  the  sense  that  some  error  may  be  acceptable  in  order  to  simplify  the 
analysis  by  increasing  block  size. 

Obviously,  error  is  encountered  when  large  logic  blocks  are  modeled 
by  assigning  single  delays  to  each  output.  Can  various  strategies  for 
characterizing  the  delays  be  employed  to  trade  between  simplicity  and  the 
degree  of  error?  This  is  discussed  next. 


4.2  Output  Delay  Characterization 


A  difficult  problem  that  arises  in  the  simple  output  delay  model 
deals  with  characterizing  the  output  delays.  Focusing  attention  on 
a  particular  output,  there  is  usually  more  than  one  block  input  lead¬ 
ing  to  that  output.  Also,  more  than  one  path  may  exist  between  an 
input  and  the  output.  (In  such  a  case  the  network  is  said  to  have 
reconvergent  fanout).  Since  different  paths  are  likely  to  have  differ¬ 
ent  delays,  any  strategy  for  characterizing  the  output  delay  pair  {x 
results  in  an  approximation  to  the  large  logic  block.  The  problem  is 
further  complicated  because  delay  observed  at  an  output,  in  general, 
depends  on  values  of  the  inputs.  This  latter  point  is  illustrated  in 
Example  4.1. 

Example  4 ■ 1  Consider  the  circuit  shown  in  Fig.  4.2.  This  is  a  2 


Fig.  4.2.  A  2  out  of  3  majority  logic  circuit. 

out  of  3  majority  logic  circuit.  It  will  be  investigated  throughout  this 
chapter,  including  the  computer  simulations.  In  this  example,  attention 


is  focused  on  a  rise  transition  propagating  from  the  input  node  X3 
to  the  output  node  Let  XI (tj)  -1 ,  X2(t^)=0.  In  this  case  a  rise 
transition  in  X3(t)  which  occurs  at  t®^  propagates  through  the  gates 
01,  A1 ,  and  02,  and  characterization  of  the  output  delay  is  determined 
by  the  p.d.f.'s  associated  with  the  rise  propagation  delays  of  these 
three  gates.  Next,  let  XJ_(t^)=0  and  5(2(t^)  =  l.  Now  a  rise  transition 
in  X3(t)  at  t=t^  propagates  through  the  gates  A2  and  02,  and  character¬ 
ization  of  the  output  delay  is  determined  by  p.d.f's  associated  with 
the  rise  propagation  delays  of  these  two  gates.  It  is  seen  that  the 
delay  characteristics  observed  at  the  output  can  depend  on  the  input 
values.  V, 

The  difficulty  illustrated  by  Ex.  4.1  can  be  overcome  by  assigning 
p.d.f's  to  the  output  which  are  dependent  on  the  inputs.  However,  this 
becomes  very  cumbersome.  Consequently,  effort  in  this  chapter  is 
devoted  to  developing  strategies  for  assigning  p.d.f.'s  at  the  outputs 
which  are  not  input  dependent. 

Four  strategies  for  characterizing  the  output  random  delays 
by  assigning  inpi  indendent  p.d.f.'s  are  considered  next. 

1)  Longest  path  delay  In  this  strategy  the  longest  path  is  used 
to  characterize  the  output  rise  and  fall  random  delays.  By  definition, 
the  longest  path  is  that  which  contains  the  largest  number  of  gates. 

If  two  or  more  paths  contain  the  largest  number  of  gates,  the  longest 
path  may  be  selected  by  applying  a  suitable  criterion  such  as  the  sum 
of  mean  delays  or  the  sum  of  maximum  delays  along  each  path. 


The  characterization  of  the  output  rise  and  fall  transition  delays 
must  be  handled  separately.  First,  the  logic  block  inputs  are  selected 
such  that  transition  propagates  from  the  path  input  to  output.  As 
the  transition  propagates,  it  is  noted  whether  a  rise  or  fall  transition 
occurs  at  the  output  of  each  gate  along  the  path.  The  total  path  delay 
for  the  rise  or  fall  at  the  output  equals  the  sum  of  the  rise  and  fall 
propagation  delays  experienced  along  the  path.  Assuming  statistical  inde¬ 
pendence  from  gate  to  gate,  the  p.d.f.  of  the  total  delay  is  obtained 
by  convolving  the  p.d.f's  of  the  delays  in  the  sum.  This  procedure 
is  repeated  so  that  p.d.f's  are  obtained  for  both  rise  and  fall  trans¬ 
ition  delays  at  the  path  output. 

Example  4.2.  Consider  the  circuit  shown  in  Fig.  4.2.  The  longest 
path  in  the  circuit  includes  the  gates  01,  A1  and  02.  It  is  desired  to 
obtain  p.d.f. 's  for  both  rise  and  fall  transition  delays  at  the  longest 
path  output. 

Assume  the  logic  block  inputs  have  been  selected  such  that  a  rise 
transition  is  generated  at  the  output  of  the  gate  01  and  propagates  to 
the  output  ^(t)  along  the  longest  path.  By  inspection,  a  rise  transi¬ 
tion  at  the  input  to  gate  A1  will  create  a  rise  transition  at  its  output. 
This  rise  transition  will,  in  turn,  produce  a  rise  transition  at  the  outpu 
of  gate  02.  Denote  the  rise  transition  delays  for  gates  01,  A1  and  02 
and  for  the  longest  path  by  ,  t^q.,  and  J.r^»  respectively.  It 

follows  that 


-rT  — rOl  +  — rAl  +  -r02  ' 


(4.1) 


Similarly,  let  the  fall  transition  delays  be  denoted  by  ,  — fAl * 
t_£Q2  and  XfT*  an  anal°S°us  manner,  it  is  readily  shown  that 


-fT  ^01  +  —fAl  +  — f 02  ' 


(4.2) 


Let  the  delays  of  the  same  type  of  gates  be  identically  distributed. 
Assuming  statistical  independence  among  delays  of  different  gates,  the 


>.d.f.'s  of  the  total  delays  are  given  by 


fxrX(T>  =  fxr0(T)*fTrA<T)*fxr0(T)* 


fxfT(l)  fxf0('T)*fxfA(T)*fxf0(T)  * 


(4.3) 


If  the  entire  circuit  of  Fig.  4.2  is  considered  as  a  large  logic 
block,  the  internal  delays  are  ignored  and  a  single  discriminating 
delay  element  is  assigned  to  the  output  node  Z.  This  delay  element  has 
X  and  x_£T  as  its  rise  and  fall  transition  delays,  respectively, 
with  p.d.f.’s  given  by  Eqs.(4.3).  V 


2)  Mean  path  delay.  In  this  strategy,  an  expression  is  obtained 
for  the  arithmetic  mean  delay  of  all  possible  propagation  paths  from 
the  inputs  to  the  output.  The  arithmetic  mean  delay  is  expressed  in 
terras  of  the  delays  associated  with  the  various  gates  in  the  circuit. 
It  is  obtained  by  computing  the  total  delay  for  every  path  from  any 
input  to  the  large  block  output  (as  was  done  in  strategy  no.  1  for  the 
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longest  path) ,  summing  all  the  total  delays  of  the  paths  and 

dividing  by  the  number  of  paths.  The  following  notation  is  used. 

th 

For  an  output  rise  transition  associated  with  the  j  path  beginning 

at  the  input  node  Xk,  the  total  propagation  delay  is  denoted  by  T 

— rkj 

Similarly,  for  an  output  fall  transition  associated  with  the  j1"*1  path 
beginning  at  input  node  Xk,  the  total  propagation  delay  is  denoted  by 
— f kj '  Let  n  denote  t*ie  number  of  input  signals  and  let  n^  denote  the 
number  of  distinct  paths  from  the  input  node  Xk  to  the  output  node  Z. 
The  arithmetic  mean  delay  for  the  output  rise  transition  is  denoted  by 
and  is  given  by 


n  ^  n 

t  =  y  £  t  /  £  n  . 
k‘i  3ti  k-l  k 


(4.4) 


Similarly,  the  arithmetic  mean  delay  for  the  output  fall  transition, 
denoted  by  T_^,  is  given  by 


ni 

n  k 


-f  “  J-  J  -fkj/ J-  nk" 

k=l  j-1  J  k=l 


(4.5) 


Note  that  jr^j  and  JT^j  are  sums  of  rise  and  fall  propagation  delays 

£ll 

associated  with  gates  along  the  j  path  beginning  at  node  It  follows 

that  and  are  linear  combinations  of  the  rise  and  fall  propagation 
delays  of  the  various  network  gates.  Assuming  statistical  independence 


among  delays  of  different  gates,  the  p.d.f.'s  °f  J_r  and  can  be  obtained 
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using  the  following  rule.  Let  Y.  ■  £  a.  xi,  where  the  a.  s  are  constants 

i-1 

and  the  xi's  are  statistically  independent  random  variable  with  given 
p.d.f's  fx^(xi).  Extending  the  result  in  [18,  p.190]  to  several 
variables,  the  p.d.f  of  Y  is  given  by 


f y(Y)  = 


'n  a.1 
i-1 1 


f .  (I)*. 

xl 


.*f  (r-  ) 

xn  a 
—  n 


(4.6) 


The  method  is  illustrated  by  the  following  example. 


Example  4.3.  Consider  the  circuit  in  Fig.  4.2.  Note  there  are  five 
distinct  paths  from  the  inputs  to  the  output  node  Z.  It  is  desired 
to  obtain  in  terms  of  the  p.d.f's  of  the  gate  delays  the  p.d.f.'s  of  the 
mean  path  propagation  delay  for  rise  and  fall  transitions  at  the  output 


The  method  used  in  Ex.  4.2  to  determine  the  total  delay  for  the 

longest  path  is  used  here  to  evaluate  the  total  rise  transition  delay 

for  each  of  the  five  paths.  Recalling  that  Lj-^j  denotes  the  total  rise 

th 

transition  delay  for  the  j  path  beginning  at  input  Xk,  it  follows  that 


-T-rll  — rAl  +  -rOr 


-“r21  "  irOl  + -rAl  +  ^r02’ 


— r22  — tA2  +  — r02  ’ 


— r31  ”  — rOl  +  -rAl  +  -r02* 


— r32  =  ^rA2  +  -r02  * 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


(4.11) 


•  *'*  .*•  /•  .*  -  •  .*■  .*•  .*-*.‘*  >  •**  .*  .*  .**  •  *“  -**  ‘  -  ■  .**  /  . 
-  *  .*  •  *  ./o  V*  •  ■  .  •  .  *  .  -  •  *  -  -  .  •  *  *  •  *  .  •  ^  V*  ^  .  •  .  *  O  '  ■ 


Summing  Eqs. (4 . 7) — (4 -11)  and  dividing  by  the  number  of  paths  as  in 
Eq.(4.4),  yield 


T  =  (T  ..  +  I  „„  +  T  +  T  ..  +  X  +  T  +  T  +  T  +  X  +  T 

— r  — rAl  — r02  — rOl  — rAl  — r02  — rA2  — r02  — rOl  — rAl  — r02 


+  — rA2  +  ^r02)/5 


1  2  2 
=  —  T  +—  T  +—  T  +  T 

5  —rAl  5  —tA2  5  -±r01  -^02 


Application  of  Eq.  (4-6)  to  Eq.  (4.12)  results  in 


f?r  <£>  *  IT  <f>*£,r0  <f>‘£TfO<£) 


A  similar  derivation  for  results  in 


£Tf<T)  ■  TX  £TfA<T)',£T£A<f>*£TfO<f>*fTfO(T) 


(4.12) 


(4.13) 


(4.14) 


Although  Eq.(4.13)  and  (4.14)  are  of  the  same  form,  the  p.d.f’s  for 

and  may  differ  significantly  depending  on  the  gate  rise  and  fall 
delay. 


3.  Weighted  mean  path  delay.  As  in  the  case  of  mean  path  delay, 
all  paths  from  every  input  to  every  output  are  considered.  However, 
each  path  is  now  weighted  according  to  the  number  of  transitions  likely 
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to  occur  in  the  path  during  the  observation  interval.  Development  of 
the  weighting  scheme  is  considered  next.  As  in  Eq.  (2.79),  let  the 
derivative  of  the  counting  signals  associated  with  the  output  of  the 
ideal  logic  circuit  be  expressed  by 


£+(t)  -  l  (Xk+(t)Bk+(t)  +  Xk  (t)Bk"(t)} 
k=l 

and 

n  .  #  , 

’£ (t)  =  \  {Xk  (t)Bk“(t)  +  Xk“(t)Bk  (t)}  .  (4.15) 

k=! 


The  coefficients  Bk+ (t)  and  Bk  (t)  are  assumed  to  be  in  their  arithmetic 
canonical  form  as  presented  in  Eqs . (2 .72) , (2 .74) ,  and  (2.75).  Recall  that 
Pk(m)  contributes  to  either  Bk+  or  Bk  ,  but  not  both,  and  consists  of  a  dis¬ 
tinct  combination  of  input  variables  or  their  complements. 

When  Pk(m)  =  1,  either  Bk+(t)  or  Bk~(t)  equals  unity,  but  not  both. 

From  Eq . (4 . 15)  it  follows  that  a  transition  in  Xk(t)  will  produce  a  trans¬ 
ition  in  £(t) .  Consequently,  when  Pk(m)  =  1,  the  gates  in  the  network  are 
enabled  such  that  a  transition  in  Xk(t)  propagates  from  the  node  Xk  to  the 
output  node  .  Several  such  paths  may  be  activated  simultaneously  when 
_Pk(m)  =  1.  Denote  the  collection  of  such  paths  by  With  reference 

to  Eqs. (4. 15)  it  is  clear  that  a  rise  transition  occurs  at  the  output 
of  each  path  in  provided  either 

1)  Pk(m)  =  1  yields  Bk+(t)  =  1  and  Xk(t)  has  a  rise  transition 
or 

2)  _Pk(m)  =  1  yields  Bk~(t)  =  1  and  Xk(t)  bas  a  fall  transition. 


km 


Similarly,  a  fall  transition  occurs  at  the  output  of  each  path  in  tt 


provided  either 

1)  Pk(m)  =  1  yields  Bk+ (t)  =  1  and  Xk(t)  has  a  fall  transition 
or 

2)  ^k(m)  =  1  yields  Bk  (t)  =  1  and  Xk(t)  has  arise  transition. 

Denote  by  the  event  that  the  collection  of  paths  tt  are  activated 
such  that  transitions  in  Xk(t)  can  propagate  to  the  output  z(t).  This 
is  equivalent  to  the  event  that  Pk(m)  =  1.  Consequently,  the  probability 
that  the  collection  of  paths  ir^  are  activated  is  given  by 

Pr  {Akm^  =  =  E[Pk(m)].  (4.16) 

Attention  is  now  devoted  to  determining  the  expected  rate  of 
transitions  in  Xk(t) .  To  be  consistent  with  Eq.(3.40)  and  (3.41)  it  is 
assumed,  without  loss  of  generality,  that 

Xkf(O)  =  0 

2^+(°)  =  Xh(°) •  (4.17) 

Given  the  sample  function  Xk(t),  it  follows  for  t  >  0  that  Xk+(t  )  equals 

A  A 

the  number  of  rise  transitions  in  the  interval  [0,t.l  while  Xk_(t.) 

A  A 

equals  the  number  of  fall  transitions  in  [0,t  ].  The  rates  of  rise  and 

A 

fall  transitions  in  Xk(t)  at  some  time  t  e  [0,t  ]  are  given  by  Xk+(t) 

A 

•  _ 

and  Xk  (t) ,  respectively.  Taking  the  expected  values  over  the  ensembles 

♦  •  _ 

°f  and  Ilk  (t)  and  interchanging  the  order  of  expectation  and  dif¬ 

ferentiation,  as  in  Eq.  (2.81),  it  follows  that  E[in(.+  (t)]  and  E[Xk~(t)]  are 
the  expected  rates  of  rise  and  fall  transitions  in  Xk(t)  respectively. 


At  any  time  t  the  expected  rate  of  input  rise  and  fall  transitions 

which  propagate  through  the  collection  of  paths  tt,  are  given  by 

km 

E[Xk+(t) ]E[Pk(m) ]  and  E[Xk  (t) ]E[Pk(m) ] ,  respectively.  The  output 
transition  being  a  rise  ot  fall  depends  upon  whether  Pk(m)  ■  1  causes 
Bk+(t)  =  1  or  Bk  (t)  =  1.  Let  the  time  average  of  a  function  h.(t) 
over  an  interval  (tQ,t^)  be  given  by 


(< 

V _to  t 


h(t)dt , 


(4.18) 


Taking  time  averages  over  the  observation  interval,  denote  the  time 
averages  of  the  expected  rate  of  output  rise  and  fall  transitions  due 


to  transitions  propagated  through  the  paths  by  arkm  and  a  ,  respec¬ 


tively.  It  follows  that 


f  .  +  + 

<  Elxk  (t) ]E(Pk(m) ]  >  ,  Pk(m)  =  1  =>  Bk  -  1 

<  E[Xk~(t) ]E[Pk(m) ]  >  ,  Pk (m)  =  1  =>  Bk~  =  1 .  (4.19) 


and 


E[Xk  (t) ] E[Pk(m) ]  > 
E[Xk+ (t) ]E[Pk(m) ]  > 


,  Pk(m)  =  1  =>  Bk+ 
,  Pk (m)  =  1  =>  Bk 


-  1. 

-  1.  (4.20) 


a^n,  an<^  aficm  are*  in  a  sense»  measures  of  usage  for  the  paths  n  . 

To  determine  the  rise  and  fall  delays  associated  with  each  of  the  paths 
in  tt  ,  it  is  necessary  to  identify  the  various  gates  included  in  each 


path.  This  is  done  by  first  assigning  to  the  inputs  X j ,  j=l,...,n;  j  t  k. 


those  values  which  result  in  Pk(m)=l •  By  tracing  each  path  from  the 
input  node  Xk  to  the  output  node  £  ,  it  is  noted  which  gates 
experience  transitions  due  to  a  transition  in  Xk(t).  Paths  for 
which  transitions  do  not  propagate  completely  to  node  z^  are  ignored. 

Denote  by  D  ^  and  the  number  of  paths,  respectively,  which 

produce  rise  and  fall  transitions  in  z^(t)  due  to  a  transition  in 
Xk(t).  Using  the  same  method  as  in  strategies  1  and  2,  a  total  delay 
is  determined  for  each  one  of  the  0-^^  +  M^^)  paths.  Denote  by 
— rkm  an<*  — fkm  t*ie  arithmetic  mean  delays  for  the  Pr  and  paths, 

respectively . 

Let  Y^  denote  the  number  of  nonzero  coefficients  Bk(m)  in  the  summation 
for  Bk  given  by  Eq.(2.72).  Y^  equals  the  number  of  distinct  collections 
of  paths  Tt  from  the  input  node  Xk  to  the  output  node  z_.  Taking  into 


consideration  each  of  the  inputs  to  the  large  logic  block  and  performing 

n 

a  weighted  arithmetic  mean  over  the  \  Y,  path  delays,  the  weighted 

k=l  k 

mean  path  delays  for  the  output  rise  and  fall  transitions  in  z(t). 


<T  >  and  <T_> ,  are  defined  to  be 
— r  — f 


n  Yk 


<  V  =  l  l 


k= 1  m= 1 


3rkm  —rkm 


/n  Yk 

A  l 

k=l  m=l 


(4.21) 


n  \ 


<  "  l  l 


k=l  m=l 


afkm  —rkm 


/n  \ 

A  £ 

/k=l  m=l 


(4.22) 


Because  of  the  weights  employed  in  Eqs.  (4.21)  and  (4.22),  <T^_>  and 
<T_^>  reflect  the  usage  of  the  paths  t>^. 
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As  in  Eqs.(4.4)  and  (4.5),  <T^>  and  <t_£>  are  linear  combinations 
of  rise  and  fall  propagation  delays  associated  with  individual  gates. 
Assuming  statistical  independence  among  delays  of  different  gates,  the 
p.d.f.'s  f<Ttf>(x)  and  f<T£>(T)  can  be  obtained  by  application  of  Eq.(4.6). 
The  previous  discussion  is  illustrated  by  the  following  example. 

Example  4.4.  Consider  again  the  2  out  of  3  majority  logic  circuit  shown 
in  Fig.  .4.2.  The  switching  operation  performed  by  this  circuit  is 
Z  -  X1*X2  v  XI *X3  v  X2*X3.  It  is  desired  to  obtain  expressions  for 
f<^r>(T)  and  when  this  circuit  is  considered  as  a  single  large 

logic  block. 

The  large  logic  block  model  for  the  circuit  in  Fig.  4.2  is  depicted 
in  Fig.  4.3.  The  output  of  the  ideal  logic  circuit  is  given  by 


Fig.  4.3.  Large  logic  block  model  for  circuit  in  Fig.  4.2 

£(t)  =  JU(t)*X2(t)  v  XI_(t)  -X3(t)  v  X2(t)  *X3(t)  .  (4 

•  ^  •  •» 

First,  the  derivatives  z  (t)  and  z  (t)  are  obtained  in  their  arithmetic 
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canonical  form.  An  arithmetic  expression  for  z^(t)  in  Eq.(4.23)  is  given  by 


(4.24) 


£(t)  =  )Q(t)2C2(t)  +  n(t)X2' (t)X3(t)  +  XI '  (t)X2(t)X3Ft)  . 
Differentiation  of  Eq.  (4.24)  yields 


z(t)  =  XUt)X2(t)  +  X2(t)3U(t)  +  XHt)X2 '  (t)X3 (c>  +  X2'  (t)n(t)X3(t) 

+  X3  (t)Xl_  (t)X2.  •  (t)  +  XT  (t)X2(t)X3(t)  +  X2(t)2a'  (t)X3(t) 

•  ,  N  (4.25) 

+  X3(t)  XI’  (t)X2(t)  . 

Applying  Eq.(2.62)  to  Eq.(4.25)  and  collecting  terms  result  in 

£(t)  =  h( t)  {X2(t)  +  X2' (t)X3(t)-X2(t)X3(t)}  +  X2(t){Xl_(t)  -  Xl_(t)X3(t) 

+  XI'  (t)X3(t)}  +  2C3(t){ja(t)X2' (t)+  )Q'(t)X2(t)}  (4.26) 

With  reference  to  Eq.(2.63), 

Bl(t)  =  X2(t)  +  X2'  (t)  X3(t)  -  X2(t)X3(t) 

B2(t)  =  Xl(t)  -  Xi(t)X3(t)  +Xl’(t)X3(t) 

B3(t)  =  Xl(t)  X2'(t)  +  XI'  (t)X2(t)  .  (4.27) 


As  discussed  on  page  43,  the  arithmetic  canonical  forms  for  Bk(t) , 
k=i.2,3,  are  obtained  from  Table  4.1  shown  below: 

TABLE  4.1 


li 

m  X2  X3  Bl(m)  Pl(m) 

1  0  0  0  (X2’)0a') 

201  1  (X2’)(X3) 

310  1  (X2)(X3') 

4110  ( X2 ) (X3) 


m 

1 

2 

3 

4 


B2 

XI  X3  B2(m)  P2(m) 

0  0  0  ()U’)(X3') 

0  1  1  (2a')(X3) 

1  0  1  (XJJ  (X3 ') 

1  1  0  (XV)  (X3) 


1 
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t 

TABLE  4.1  (continued) 

-  -  -w 

J 

• 

_B3 

; 

a 

m  XI  X2  B3(m)  P3(m) 

• 

1  0  0  0  (X£')(X2;) 

2  0  1  1  oa  ’)(X2) 

K 

3  10  1  (5y^)  (2^' ) 

4  11  0  (XV)  (X2 ) 

• 

fcft 

With  reference  to  Table  4.1  the  arithmetic  canonical  forms 

k  =  1,2,3,  are 

for  Bk(t) , 

• 

Bl_(t)  =X2'(t)X3(t)  +  X2(t)X3'(t) 

• 

B2(t)  =  XT' (t)X3(t)  +X£(t)X3’(t) 

B3(t)  =  XI’ '  (t)X2(t)  +  Xl(t)X2'(t). 

(4.28) 

• 

a 

From  Eqs.(2.74)  and  (2.75)  it  follows  that 

-  ....  ^ 

» 

(t)  =  X2'(t)X3(t)  +  X2(t)X3’(t) 

a 

=  £1(2)  +  ££(3)  , 

B2+(t)  =  2LL*  (t)X3(t)  +  X£  (t)3C3 ' (t) 

-  ■  ■  -A 

• 

=  £2(2)  +  £2(3) , 

(4.29) 

• 

B3+(t)  =  X£*(t)X2(t)  +  XI  (t)  X2/(t) 

=  £3(2)  +  £3(3)  , 

Bl(t)  =  B2~(t)  =  B3"(t)  =  0. 

• 

• 

The  substitution  of  Eqs.(4.29)  into  X4.15)  results  in 

• 

z+(t)  =  l  Xk+(t)[Pk(2)  +  Pk(3) ] 
k=l 

* 

• 

3 

z(t)  =  l  Xk"(t)[Pk(2)  +  Pk(3) ] . 
k=l 

(4.30) 

• 

• 
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The  conditions  under  which  Pk(m)  =  1  for  k 
tabulated  in  Table  4.2. 


1,2,3  and  m  =  2,3  are 


TABLE  4 

.2 

PI  (2) 

PK3) 

P2  (2) 

P2(3) 

P3  (2) 

P3(3) 

(X2 ' ) (X3) 

(X2) (X3 ' ) 

(XI') (X3) 

(XI) (X3 ' ) 

(XI ' ) (X2) 

(XI) (X2 1 ) 

n 

- 

0 

1 

0 

1 

X2  0 

1 

- 

- 

1 

0 

X3  1 

0 

1 

0 

_ 

Referring  to  Eqs.(4.19)  and  (4.20)  and  utilizing  statistical  independence, 
the  expressions  for  the  weig.  a  a^^  and  afkm  are  given  by 

arl2  =  "  EfAL+(t)]  E  tX2’ (t) ]  E[X3(t) ]  > 

arl3  =  E(Xi+(t)J  E[X2( C)]  E(X3'(t)]  > 

ar22  =  <  E[X2+(t)]  E[Xl'(t)]E[X3(t)]  > 

ar23  =  <  E[X2+(t)]  E[Xj.(t)]  E[X3'(t)]  > 


ar32  KlXT'(t)]E[X2(t)] 

ar33  =  "  El?i3+Ct )  ]  E[XUt))E[X2’(t)] 
and  . 

afl2  =  <  E[Xl_~(t)]E[X2'(t)]E[X3(t)] 


> 


(4.31) 

> 


afl3  =  El2LL~(t>]E[X2(t)]E[X3’(t)]  > 


af22  =  "  E[X2'(t)]E[Xl_'(t)]E[X3(t)) 


af23  =  <  E[X2  (t)lE[Xl_(t)]E[X3’(t)} 


k 

m 

Gates  in 
Path  tt 

Kill 

PVkm 

— irkm  ''fkm 

— f  km 

n 

1 

2 

Al  ,02 

1 

— rAl 

+  T  1 

— r02 

T  +  T 

“f  A 1  —f02 

• 

1 

3 

A1 ,02 

1 

—rAl 

+  T  1 

— r02 

— f Al  +  — f  02 

a 

2 

2 

A2 ,02 

1 

-rA2 

+  T  1 

~r02 

— f A2  +  — f02 

• 

2 

3 

01 ,A1 ,02 

1 

— rOl 

+  AAl+ir02  1 

-f01+-fAl+-f02 

3 

2 

A2 ,02 

1 

— rA2 

+  A02  1 

T  +  T 

— f  A2  — f 02 

a 

3 

3 

01 ,A1 ,02 

1 

— rOl 

+  irAl%02  1 

-f  l+-fAl+-f02 

• 

The  next  step  is  to  identify  the  paths  tt  and  determine  the 
total  rise  and  fall  propagation  delays  associated  with  each  path. 

The  gates  in  each  path,  hrkm,  UfkuT  and  — fkm  are  listed  in  Table 

4.3.  To  illustrate  how  entries  in  the  Table  are  obtained,  focus  attention 

TABLE  4.3 


on  the  cases  k  =  2,  m  =  2  and  k  =  2,  m  =  3.  For  k  =  2,  m  =  2,  P2 (2) 

=  (XI ' ) (X3)  =  1  provided  =  0,  1^3  =  1.  For  these  values  of  XI (t)  and 
X3(t),  the  gate  A1  is  inhibited  while  the  gates  A2  and  02  are  enabled. 
Therefore,  a  transition  in  l(2(t)  propagates  through  the  path  9.,  consis¬ 
ting  of  gates  A2  and  02.  Observe  that  rise  and  fall  transitions  in 
}£2(t)  produce  rise  and  fall  transitions,  respectively,  at  the  outputs 
of  gates  A2  and  02.  For  k  =  2,  m  =  3,  P2^(3)  =  (X 1 )  (X3 ' )  =  1  provided 


-v.'-iktJk.V.  V.nV.  ■  -« 
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XI =1 ,  X3=0.  For  these  values  of  Xl^t)  and  X3(t) »  the  gate  A2  is 
inhibited  while  gates  01,  Al,  and  02  are  enabled.  Therefore,  a  trans¬ 
ition  in  X2(t)  propagates  through  the  path  consisting  of  gates  01, 

Al,  and  02.  As  before,  rise  and  fall  transitions  in  X2(t)  produce  rise 
and  fall  transitions,  respectively,  at  the  outputs  of  gates  01,  Al,  and 
02. 

Finally,  the  results  from  Eqs.(4.31),  (4.32),  and  Table  4.3  are  sub¬ 
stituted  into  Eqs.(4.21)  and  (4.22)  to  obtain  expressions  for  <x  >  and 

— r 

<T.f>.  Eq.(4.6)  is  then  employed  to  obtain  the  output  p.d.f's  x<t  >(i) 

and  f<Tf>(T)- 

4.  Assigned  Gaussian  delay.  The  assigned  Gaussian  delay  strategy 

was  motivated  by  Magnhagen  and  Flisberg[12]  who  assigned  Gaussian  p.d.f.' 

to  each  gate  delay  in  their  simulation  program  called  D1GSIM.  However, 

in  the  strategy  proposed  here,  the  large  logic  block  is  modeled  by  simply 

assigning  rise  and  fall  Gaussian  delays  at  each  block  output.  The  Gaussi 

parameters  are  obtained  as  discussed  below. 

The  mean  and  variance  for  each  Gaussian  p.d.f.  are  determined  from 

the  minimum  and  maximum  delays  experienced  bv  transitions  propagating  fro 

any  input  to  the  output  under  consideration.  Denote  these  rise  and  fall 

delays  by  T  .  ,  T  ,  x.  .  ,  and  i.  .  It  is  assumed  that  minimum  and 
rmin  rmax  fnun  fmax 

maximum  delays  are  known  for  each  gate  within  the  large  logic  block.  i 
is  determined  by  1)  identifying  the  paths  with  the  minimum  number  of  gate 
2)  evaluating  the  minimum  output  rise  propagation  delay  for  each  path  by 


summing  the  appropriate  minimum  delay  for  each  gate  along  the  path,  and 

3)  assigning  to  x^^  the  minimum  of  the  sums  evaluated  in  step  (2). 

Tfmin  ■'-S  determined  in  like  manner.  Tmax  is  determined  by  1)  identifying 

the  paths  with  the  maximum  number  of  gates,  2)  evaluating  the  maximum 

output  rise  propagation  delay  for  each  path  by  summing  the  appropriate 

maximum  delay  for  each  gate  along  the  path,  and  3)  assigning  to  Tmax 

the  maximum  of  the  sums  evaluated  in  step  (2) .  A  similar  procedure  is 

used  to  obtain  T, 

fmax 

Let  the  random  variables  for  the  rise  and  fall  delays  assigned  to 
the  output  be  denoted  by  and  _T^  ,  respectively.  By  definition, 
the  means  are  selected  to  be  the  averages  of  the  maximum  and  minimum  delays. 
Hence, 


E[x  _]=-=■  (t  +  t  .  ) 
— rG  2  rmax  nan 


E[-fG]  2  (Tfmax  +  Tfmin' 


(4.33) 

(4.34) 


The  variances  are  determined  by  assigning  the  3 a  points  of  the  Gaussian 
p.d.f.  to  the  maximum  and  minimum  delays.  Therefore, 


E[t  _]  -  3a  _  =  X  . 

— rG  xrG  rmin 


E[x  ]  +  3a  =  x 

rG  TrG  rmax 


(4.35) 


and 


Ef-fG1  ^TfG  Tfmin 


E  [  x  ]  +  3a  =  x 
— fG  xfG  fmax. 


(4.36) 
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Solution  of  Eqs.  (4.35)  and  (4.36)  for  the  standard  deviations  result 


O  =  -T  (T  -T  .  ) 

xrG  6  rmax  rmin 


°xfG  6  ^Tfmax  Tfmik 


(4.37) 


(4.38) 


The  approximated  Gaussian  p.d.f.'s  for  the  output  delays  are  given 

(  18[t  -  -kx  +x  .  )]2\ 

(  r  2  rmax  ruin  \ 


(t  -  x  .  )2 
rmax  rmin 


f  (x  )  =  - - - e 

^rG  r  /^T  (x  -x  .  ) 
rmax  rmin 


(4.39) 


f  (V 

xfG 


18[Tf-I  (lf„ax+Wl2 

(x -  -  X  )2 

fmax  fmin 


/27T^Tfmax  Tfmin^ 


(4.40) 


The  strategy  is  illustrated  by  the  following  example. 

Example  4.5.  Consider  again  the  circuit  in  Fig.  4.2.  Suppose  the  minimum 
and  maximum  values  for  the  rise  and  fall  propagation  delays  associated  with 
the  AND  and  OR  gates  are  known  to  be 


0. 1  <  x  <  0.6 
— rA 


0.2  <  x,.  <  0.5 
— f  A 


(4.41) 


0.2  <  x  .  <  0.9 

— rO 


0.3  <  x_fQ  <  1 .0 


(4.42) 


where  it  is  assumed  that  like  gates  have  the  same  minimum  and  maximum 
delays.  It  is  desired  to  assign  Gaussian  p.d.f.'s  to  the  output  delays 
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when  the  entire  circuit  of  Flg>  4.2  Is  modeled  as  a  large  logic  block* 
For  the  circuit  of  Fig.  4.2,  the  minimum  number  of  gates  in  a  path 
equals  two.  There  are  three  paths  with  two  gates:  1)  a  path  from  node 
XI  to  node  Z  via  gates  A1  and  02,  2)  a  path  from  node  X2  to  node  Z 
via  gates  A2  and  02,  and  3)  a  path  from  node  X3  to  node  Z  via  gates  A2 
and  02.  Observe  that  rise  and  fall  transitions  propagate  in  the  same 
manner  through  each  path.  Since  all  three  paths  contain  one  AND  gate 
and  one  OR  gate  and  like  gates  have  the  same  minimum  delays,  only  one 
of  the  three  paths  need  be  analyzed. 

It  follows  that 


T  .  =  min  { T  . }  +  min  { T  . } 

rmrn  — rA  — rO 


=  0.1  +  0.2  =  0.3  , 


Tfmin  min  f A^  +  min  {-f0} 


=  0.2  +  0.3  =  0.5 


(4.43) 


(4.44) 


The  maximum  number  of  gates  contained  in  a  path  equals  three. 

There  are  two  paths  with  three  gates:  1)  a  path  from  node  X2_  to  node 
via  gates  01,  A1  and  02,  and  2)  a  path  from  node  X3  to  node  Z^  via  gates 
01,  A1  and  02.  As  before,  rise  and  fall  transitions  propagate  in  the  same 
manner  through  each  path.  Since  both  paths  contain  the  same  gates,  it 
follows  that 


X  =  max  { T  „ }  +  max  {x  .}  +  max{x 
rmax  — rO  — rA  — rO 


=  0.9  +  0.6  +  0.9  =  2.4, 


(4.45) 
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Tfmax  =  maX  {ifO}  +  m3X  {^fA}  +  max  {ifO} 


=  1.0  +  0.5  +  1.0  =  2.5. 


Substitution  of  T  .  and  x  in  Eq.(4.39)  and 
rmxn  rmax 


Eq.  (4.40)  result  in 


18[xr-1.35r 


f^(TJ  =  - ^  e  (2 . 1)  2 


xrG  r 


2.1/2n 


=  1.14  e 


4.08(x  -1.35)‘ 


18 [x^-1 -5]' 


fxfG(xf) 


6  e 


2/2v 


=  1.20  e 


-  4.5(xf-1.5)‘ 


The  p.d.f.'s  for  X  n  and  x,  ,  are  sketched  in  Fig 

■ — ro  — iLj 


As  a  final  point,  it  is  noted  that  computer  simulation  programs,  such  as 


LDPS  by  A1 -Hussein  and  Dutton  [14],  may  be  employed  to  obtain  x  .  ,  t  , 

min  max 

X  _  ,  ,  and  xc 
fmin  fmax 

4.3  Comparison  of  Output  Delay  Characterizations 

The  four  strategies  discussed  in  Sec.  4.2  are  compared  in  this  section. 
Specific  p.d.f.'s  are  assigned  to  the  individual  gates  in  the  circuit  of 
Fig.  4.2.  The  circuit  is  modeled  as  a  large  logic  block  using  the  simple 
output  delay  model  shown  in  Fig.  3.10.  For  each  strategy,  the  p.d.f.'s 
of  T_r  and  are  determined  and  E[Z^(t)]  is  evaluated.  In  addition,  a  Monte 
Carlo  simulation  is  perfomed  to  provide  a  reference  waveform  for  E[Z(t)]. 

An  exact  analysis  of  the  circuit  was  not  perfomed  because  of  its 
complexity  (see  App.B).  Several  considerations  entered  into  the  design  of  the 
Monte  Carlo  simulation.  The  first  step  in  each  run  was  to  assign  numerical 
values  to  the  gate  rise  and  fall  propagation  delays.  Random  number  generators 
which  reflected  the  assigned  delay  p.d.f.'s  were  used  for  this  purpose.  A 
new  set  of  gate  delays  was  randomly  and  independently  selected  for  each 
run  and  the  output  signal  obtained.  Let  the  output  of  the  run  be 

denoted  by  Z^(t).  Also,  let  N  denote  the  total  number  of  runs.  The  reference 
waveform  derived  from  the  Monte  Carlo  simulation  was  defined  to  be 

1  N 

ZM(t)=N  l  Zi(t)-  (4’49) 

i=l 

The  accuracy  of  the  Monte  Carlo  simulation  is  now  discussed.  Consider 
a  particular  time  instant,  t^,  in  the  observation  time  interval.  Denote 


by  Pj  the  probability  that  2^(t^)  ■  1.  Specifically, 

P1  =  pr  -  1}  =  EfZCtp].  (4.50) 

Since  Z(t  )  is  a  binary  random  variable  with  values  0  or  1,  the  Monte 
Carlo  simulation  for  t  =  t^  can  be  interpreted  as  repeated  Bernoulli  trials 
with  probability  P^[18,  Sec. 3-2]  assuming  the  runs  of  the  Monte  Carlo  simula¬ 
tion  to  be  statistically  independent.  Let  denote  the  number  of  runs  in  which 

_Z(tj)  =  1.  Because  _Z(t^)  is  a  0,1  binary  random  variable  and  Z^(tj)  denotes 
,  N 

the  i  sample,  N  =  £  Z.(t  ).  By  the  law  of  large  numbers  [18], 

i=l  1  1 

N 

£im  _ 1_  _  £im  J_  £  Z.(t.)  =  P.  .  (4.51) 

N-*°°  u  N-x«  N  i=l  ^  1 


Comparison  of  Eqs.  (4.49)  and  (4.51)  reveals  that  Z^(t^)  is  a  consistent 
estimator  of  P^ .  The  second  order  moment  of  2^(tp  is  given  by 

E[Z2(t1)]  =  (D2Pr{Z(t1)  =  1}  +  (0)2pr(Z(t1)  =  0}  =  Pr  (4.52) 

This  result  can  also  be  obtained  using  the  idempotency  property  of  the  0  >1 

2 

binary  random  variable.  Since  Z^  (t^)  =  ,Z(t^),  it  follows  that 


E[r(t:)]  =  E[Z(tl)]  =  Pj  . 


(4.53) 


The  variance,  o  ,  .,  is  given  by 

^  t  ^  ) 


=  E[Z2(t1)]^E[Z(t1)]}2  =  Pl  -  PL2  =  Pj(l  -  Pj) 


(4.54) 
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When  the  Monte  Carlo  simulation  is  repeated  many  times,  the 
waveforms  in  Eq.(4.49)  can  be  viewed  as  random  processes.  Then 


w 


N  A  -i^P  ' 


i=l 


(4.55) 


The  expected  value  of  Z^  and  its  variance,  respectively  are  readily  shown 
to  be  [18,  p.246] 


Etwi=  P1 

2  Pl^ 

aZM(t.)  N 


(4.56) 

(4.57) 


Assuming  N  to  be  large  and  applying  the  central  limit  theorem  to  Eq.(4.55), 
can  be  approximated  as  a  Guassian  random  variable  with  mean  and 
variance  given  by  Eqs.(4.56)  and  (4.57). 

For  a  Guassian  p.d.f,  99%  of  the  area  is  located  in  the  symmetric 
interval  IQ  g9 ,  extending  2.6  standard  deviations  on  each  side  of  the  expected 
value.  It  follows  that 


Pr  tnzupi  -  2.6  <_Z„(tl)  <  EIZCt,)]  +  2.6  -  0.99 

(4.58) 

Therefore,  with  probability  0.99,  each  outcome  of  the  Monte  Carlo  simulation 
falls  in  the  interval  1^  Note  that  the  center  of  the  interval  is  located 

at  the  true  E[_Z(tj)].  The  width  of  the  interval  is  given  by  2e  where 


(4.59) 


1 

I 

3 


=  2.6 


l 


For  p  =  0.5,  P^(l-P^)  equals  Its  maximum  value  of  0.25.  With  N  ■  1000, 
the  interval  half  width  is  bounded  by 


e  < 


2.6 


.25 


1000 


0.041. 


(4.60) 


To  determine  a  reference  waveform,  a  single  Monte  Carlo  simualation  was  performed 
with  N  =  1000.  Eq.(4.60)  can  be  viewed  as  an  upper  bound  on  the  simulation 

error.  Values  of  C  and  e/E[Z^t^)]  are  tabulated  in  Table  4.4  for  various 
values  of  P^  =  E[Z(t^)].  It  is  seen  that  the  Monte  Carlo  simulation  yields 
a  reasonable  estimation  of  E[j£(t)]. 


TABLE  4.4 


P1 

io~3 

10-2  0.1  0.3 

0.5 

0.7 

0.9 

0.99 

0.999 

c 

.0026 

.0082  .0247  .0377 

.0411 

.0377 

.0247 

.0082 

.0026 

c 

2.6 

.82  .247  .1256 

.0822 

.0539 

.0247 

.0082 

.0026 

EfZUp] 

To 

reflect 

the  actual  behavior  of 

a  physical 

circuit , 

an  additional 

effect  was  taken 

into  consideration  in 

the  Monte 

Carlo  simulation. 

The 

Simula- 

tion  program  ignored  all  pulses  at  a  gate  input  whose  durations  <5^  were  less 
than  the  rise  propagation  delay,  x  ,  of  the  gate  and  all  gaps  whose  duration 
were  less  than  the  fall  propagation  delay,  x^,  of  the  gate.  This  effect 
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is  referred  to  as  blanking.  Indeed,  in  physical  gates,  and  typically 
approximate  the  response  times  of  the  circuit.  Therefore,  the  circuit  does 
not  respond  to  pulses  and  gaps  of  shorter  durations.  The  effect  of  blanking 
in  the  probabilistic  analysis  was  ignored  because  it  was  not  possible  to  develop 

an  analytical  procedure  for  determining  the  statistics  of  the  pulse  and  gap 
durations  as  the  signals  propagated  through  the  various  gates  in  the  circuit. 

The  results  of  this  chapter  indicate  that  blanking  in  the  large  block  can  be 
ignored  when  the  pulse  and  gap  durations  at  the  output  of  the  ideal  logic 
circuit  are  longer  than  the  maximum  delays  assigned  to  the  discriminating 
delay  element. 

As  mentioned  at  the  outset  of  this  section,  the  2  out  of  3 
majority  logic  circuit  in  Fig.  4.2  was  used  as  a  vehicle  for  comparing 
the  four  strategies.  For  convenience,  the  circuit  is  shown  again  in 
Fig.  4.5.  The  p.d.f's  of  the  delays  assumed  for  the  gates  of  this  circuit 


Yl(t) 


Y2(t) 


Fig.  4.5.  2  our  of  3  majority  logic  circuit. 


ire  sketched  in  Fig.  4.6,  and  the  deterministic  input  signals,  XI (t),  X2(t) 


and  X3(t)  are  depicted  in  Fig.  4.7. 


Fig.  4.6.  Sketch  of  p.d.f's  of  assumed  delays  for  gates  in  Fig.  4.5. 

The  circuit  is  intended  to  perform  majority  logic  on  its  three 
inputs.  In  other  words,  the  output  should  take  on  the  value  1  if  two  or 
more  inputs  equal  1  and  the  value  0  if  two  or  more  inputs  are  zero. 
Clearly,  in  the  large  logic  block  model  of  the  circuit  in  Fig.  4.5,  the 
ideal  logic  portion  will  perform  majority  logic.  z(t),  the  output  of  the 
ideal  logic  portion  is  also  sketched  in  Fig.  4.7.  Note  that  pulse  Az 
is  due  to  the  overlapping  of  pulse  XI,  of  XI (t)  and  pulse  X3,  of  X3(t). 


Similarly,  pulse  Bz  is  due  to  the  overlapping  between  pulses  X2^ 
and  X3^,  and  Cz  is  due  to  the  overlapping  between  pulses  Xl^  and  X2Z . 
Finally,  pulse  Dz  is  due  to  the  overlapping  between  1)  pulses  Xl^  and 
X2^  2)  pulses  Xl^,  X2^  and  and  3)  pulses  Xl^  and  X39  . 

Because  the  p.d.f's  associated  with  each  gate  have  finite  domain, 
the  times  of  occurrence  of  each  rise  and  fall  transition  in  the  physical 
output  Z(t)  have  finite  spreads .  Let  ,  ufiZ,  u^,  u^  and  d^,  d^z, 
d^z»  denote  the  rise  and  fall  transition  times  respectively,  of  the 

pulses  in  Z(t)  corresponding  to  the  pulses  Az,  Bz,  Cz,  and  Dz  in  z(t). 
Assuming  the  circuit  performs  as  intended,  the  spreads  for  the  random 

variables  u,„  and  d,„  are  discussed  next. 

— AZ  — AZ 

The  rise  of  pulse  Az  occurs  when  XI (t)  =  1,  X2(t)  =  0  and  X3(t) 
rises.  Because  X2(t)  =  0,  gate  A2  is  inhibited.  On  the  other  hand  gate 
AL  is  enabled  since  XI (t)  =  1.  Therefore,  the  rise  in  X3(t)  propagates 
through  gates  01,  Al,  and  02.  The  minimum  rise  propagation  delay  of 
the  AND  gates  is  observed  from  Fig.  4.6  to  be  0.1.  Similarly,  the 
minimum  rise  propagation  delay  for  the  OR  gates  is  observed  to  be  0.2. 
Hence,  the  minimum  total  rise  propagation  delay  in  this  path  is 
0.2  +  0.1  +  0.2  =  0.5.  In  a  like  manner,  the  maximum  value  of  the  total 
rise  propagation  delay  is  determined  to  be  0.9  +  0.6  +  0.9  =  2.4.  Since 
the  rise  in  X3(t)  responsible  for  the  rise  of  pulse  Az  occurs  at  t  =  2, 
follows  that  the  spread  of  u  is  given  by 


2.5  <  u . „  <  4.4. 


The  fall  of  pulse  Az  is  due  to  a  fall  in  XI (t)  while  X2(t)  ■  0 
and  X3(t)  =  1.  This  fall  transition  propagates  through  the  path  con¬ 
taining  gates  A1  and  02.  The  minimum  delay  for  fall  transitions  in  the 
AND  gates  is  seen  from  Fig.  4.6  to  be  0.2  and  the  minimum  delay  for 
fall  transitions  in  the  OR  gates  is  0.3.  Thus,  the  minimum  total 
delay  for  the  fall  transition  of  pulse  Az  is  0.2  +  0.3  =  0.5.  Similarly, 
the  maximum  total  delay  for  this  fall  transition  is  0.5  +1.0=  1.5. 
Because  the  fall  inXl(t)  responsible  for  the  fall  of  pulse  Az  occurs  at 
t  =  2.5,  it  follows  that 

3.0  '  d  <  4.0  .  (4.62) 

—AZ 

The  ranges  of  transition  times  in  Z^(t)  corresponding  to  the 
pulses  Bz ,  Cz,  and  Dz  are  obtained  in  like  manner.  The  results  are  tabu¬ 
lated  in  Table  4.5. 


TABLE  4.5 


Transit  ion 

— AZ 

-AZ 

-BZ 

-BZ 

u  „ 

- CZ 

^CZ 

^DZ 

^DZ 

min 

2.5 

3.0 

4.8 

5.5 

9.5 

11.5 

13.3 

19.5 

max 

4.4 

4.0 

6.0 

6.5 

11.4 

12.5 

14.5 

20.5 

These  spreads  are  depicted  in  Fig.  4.8. 

The  sample  functions  in  the  physical  circuit  output  Z(t)  differ, 
one  from  the  other,  depending  upon  the  actual  gate  delays  within  the 
circuit.  In  particular,  some  pulses  expected  to  appear  may  not  occur 
and  some  pulses  not  expected  to  appear  may  occur.  To  illustrate  this 


2.5  4*  4.8  6 


9.5  11.4  13.3  14.5 


-AZ  -BZ 


3.0  4  5.5  6.5 


11.5  12.5 


19.5  2  0.5 


Fig.  4.8  Spreads  of  transition  times  in  Table  4.8. 


point,  assume  the  gate  delays  are  given  by  the  average  values 

=  0.35,  =  0.35,  =  0.55,  t ^  =  0.65.  Also,  assume  each  gate 

is  modeled  as  an  ideal  logic  circuit  followed  by  a  discriminating  delay 

AAA  A  A  /S 

element.  The  resulting  waveforms  w(t) ,  W(t) ,  yl(t),  Yl(t),  y2(t),  Y2(t), 
z(t),  Z(t)  are  sketched  in  Figure  4.9.  With  reference  to  Fig.  4.7(d),  a 
pulse  AZ  corresponding  to  pulse  Az  is  expected  to  occur  in  Z(t).  However 
from  Fig.  4.8,  the  rise  transition  for  AZ  must  occur  in  the  interval 
(2.5,  4.4)  while  the  fall  transition  must  occur  in  the  interval  (3,4). 

A  A. 

Due  to  the  average  gate  delays,  the  pulse  AZ  does  not  occur  in  Z(t),  as 
can  be  observed  by  examining  the  waveforms  in  Fig.  4.9.  The  pulses  BZ, 

CZ,  and  DZ  in  Z(t)  correspond,  as  expected,  to  the  pulses  Bz ,  Cz,  and 

A  A 

Dz  of  z(t).  However,  the  pulses  EZ  and  FZ  deserve  special  attention. 

Their  rise  transitions  occur  at  8.9  and  20.5  while  their  fall  transitions 
occur  at  9.15  and  21.65.  These  transitions  times  are  not  predicted  by 
Fig.  4.8.  Consequently,  pulses  EZ  and  FZ  are  unintent ional  output  pulses 
which  occur  due  to  the  choice  of  gate  delays.  It  is  seen  that  gate 
propagation  delays  can  remove  intentional  pulses  and/or  add  unintentional 


pulses  to  the  output. 


Taking  into  account  the  spread  of  the  gate  delays,  it  is  readily 
shown  that  the  transition  times  for  pulse  EZ  and  FZ  in  Z(t)  must  satisfy 
the  inequalities 


8.3 

< 

—EZ 

< 

9.5 

8.3 

< 

-EZ 

< 

10.0 

20.3 

< 

-FZ 

< 

21.5 

20.8 

< 

-FZ 

< 

22.5 

(4.63) 


The  results  of  the  Monte  Carlo  simulation,  Z^(t),  is  plotted  in 
Fig.  4.10.  Noting  the  transition  times  it  is  possible  to  distinguish 
between  the  intentional  and  unintentional  pulses.  Observe  that  each 
pulse  spans  the  spread  predicted  by  Table  4.5  and  Eqs.  (4.63).  The  inten¬ 
tional  pulses  are  denoted  by  AZ,  BZ,  CZ,  DZ  while  the  unintentional  pulses 
are  denoted  by  EZ  and  FZ.  The  height  of  each  pulse  in  the  Monte  Carlo 
simulation  is  a  measure  of  its  probability  of  occurrence.  It  is  seen  that 
the  intentional  pulse  AZ  rarely  occurs.  This  is  also  true  of  the  uninten¬ 
tional  pulse  EZ.  Four  pulses  occur  more  than  50%  of  the  time.  These  are 
the  intentional  pulses  BZ,  CZ,  and  DZ,  and  the  unintentional  pulse  FZ. 

The  Monte  Carlo  simulation  serves  as  the  reference  waveform  for  evalu¬ 
ating  the  four  strategies.  The  evaluation  is  done  by  comparing  the  height, 
center  position,  and  width  of  each  of  the  four  intentional  pulses.  The 
pulse  width  is  defined  to  be  the  width  at  which  the  estimate  of  E[Z(t)] 
equals  0.5.  No  width  is  given  when  the  height  of  a  pulse  is  less  than  0.5. 
in  general,  the  height,  center  position,  and  width  are  denoted  by  C^, 

and  W  where  I  =  A,  B,  C,  or  D  depending  on  the  pulse  and  J  =  M,  1,  2,  3  or  4 


depending  on  the  procedure  used  to  determine  the  estimate  for  E[Z(t)]. 

M  stands  for  Monte  Carlo  simulation  while  1,  2,  3,  and  4  stand  for  strategies 
1,  2,  3,  and  4,  respectively.  Features  of  E[£(t)],  as  estimated  from  Z^Ct) 
in  Fig.  4.10,  are  tabulated  in  Table  4.6. 

TABLE  4.6 


Pulse  IZ 

AZ 

BZ 

EZ 

CZ 

DZ 

FZ 

• 

him 

0.008 

0.625 

0.086 

1.00 

1.00 

0.55 

WIM 

- 

0.4 

- 

1.6 

6.3 

0.4 

CIM 

3.10 

5.35 

9.0 

11.25 

17.0 

21.3 

• 

The  quality  of  each  strategy  for  characterizing  output  delays 
in  the  simple  output  delay  model  is  next  evaluated  with  reference  to  Z^(t)  • 
The  output  delays  for  each  strategy  were  assigned  and  the  output  expected 
value  of  the  large  logic  block  was  calculated  and  plotted  using  the  computer. 
Each  strategy  is  evaluated  by  comparing  the  features  of  each  pulse  in  the 
calculated  expected  value  with  those  given  by  Table  4.6.  Since  the  output 
of  the  simple  output  delay  model  is  a  delayed  version  of  z(t),  the  model 
cannot  predict  the  unintentional  pulses  EZ  and  FZ  no  matter  which  strategy 
is  employed.  Also,  because  pulses  Az  and  Bz  have  the  same  width  of  0.5, 
as  seen  in  Fig.  4.7(d),  pulses  AZ  and  BZ  have  identical  shapes  independent 
of  the  strategy.  These  inherent  errors  in  the  simple  output  delay  model  can 
get  very  large. 

The  computations  for  each  strategy  are  discussed  next  where  use  is 
made  of  the  results  from  the  examples  in  Sec.  4.2. 
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1)  The  longest  path  delay.  The  results  of  Eqs.(4.3)  are  used  with  the  par¬ 
ticular  p.d.f.'s  sketched  in  Fig.  4.6  .  The  resulting  p.d.f.'s,  fTrT(T)  and 
fTfT^T)  ’  anc*  e9t^matec^  expected  value  of  the  block  output  Z^(t),  denoted 

by  EjJZU)],  are  plotted  in  Fig.  4.11.  Features  of  E^^t)],  are  summarized 
in  Table  4.7. 


TABLE  4.7 


Pulse  IZ 

AZ 

BZ 

CZ 

DZ 

HI1 

0.73 

0.73 

1.00 

1.00 

WI1 

0.7 

0.7 

2.2 

6.2 

CI1 

3.9 

6.4 

11.6 

17.6 

As  mentioned  earlier,  because  pulses  Az  and  Bz  in  the  ideal  logic  block 
output,  z(t),  have  the  same  width,  pulses  AZ  and  BZ  have  identical 
shapes  in  Fig.  4.11(c)  This  is  an  inherent  error  of  the  simple  output 
delay  model  which  is  also  observed  in  the  other  three  strategies. 

The  figures  of  merit  used  in  evaluating  the  heights  and  widths  of 


the  pulses  obtained  by  the  various  strategies  for  the  estimates  of 
E[Z(t) ]  are  defined  to  be 


where  the  subscripts  I  and  J  have  the  same  interpretation  as  before. 


To  obtain  a  figure  of  merit  for  the  pulse  center  position,  let  C  denote 


sixb  x _ i '-ci 


EZ  1 


2.  50 


the  center  position  of  the  pulses  in  the  ideal  logic  output  z(t) 
where,  as  before,  I  ■  A,  B,  C,  or  D  depending  on  the  pulse.  The  figure 
of  merit  for  the  pulse  center  position  is  defined  to  be 


'IJ 


(CIM  "CIz  ) 


(4.< 


The  figures  of  merit  for  the  longest  path  delay  strategy  are  tabulated 
in  Table  4.8. 

TABLE  4.8 


Pulse  IZ 

AZ 

BZ 

CZ 

DZ 

hIl 

90.0 

0.168 

0. 

0. 

WI1 

- 

0.75 

0.375 

0.016 

CI1 

0.94 

0.95 

0.28 

0.6 

Note  that  the  figures  of  merit  are  decidedly  better  for  the  wide  pulses 
of  z(t).  Specifically,  the  figures  for  the  pulses  CZ  and  DZ,  are  much 
smaller  than  those  for  AZ  and  BZ.  In  general,  it  will  be  seen  for  all 
four  strategies  that  the  simple  output  delay  model  better  simulates 
the  physical  circuit  as  the  pulses  and  gaps  in  z(t)  become  longer. 

2)  The  mean  path  delay.  The  p.d.f.’s  f— f(r)  and  f— ^ (x )  are  deter¬ 
mined  by  carrying  out  the  convolutions  in  Eqs.  (4.13)  and  (4.14).  The 
results  are  used  to  obtain  an  estimate  for  the  expected  value  of  _Z(t), 
denoted  by  f^CZ^t)].  Fig.  4.12  shows  the  plots  of  f-^-  ("0,  f^fC1),  anc* 

E2 f-Z ( t )  ] .  Features  of  E2[Z(t)]  and  the  figures  of  merit  are  summarized  in 


Table  4.9. 
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TABLE  4.9 


Pulse  IZ 

AZ 

BZ 

CZ 

DZ 

:t 

HI2 

0.825 

0.825 

1.00 

1.00 

• 

WI2 

0.65 

0.65 

2.2 

6.1 

CI2 

3.55 

6.05 

11.3 

17.25 

a 

h  2  102.0 

0.32 

0. 

0. 

• 

WI2 

— 

0.625 

0.375 

0.032 

CI2 

0.53 

0.636 

0.04 

0.25 

r  ■ 

** 

3)  The  weighted 

mean 

path  delay. 

The  input  s 

ignals  XI (t),  X2(t)  and 

X3(t), 

• 

sketched  in  Fig. 

4.7, 

are  used  in 

Eqs.  (4.31) 

and  (4.32)  to  determine 

the 

■ 

• 

weights  a  ,  and 
°  rkra 

afkm" 

Because  the  inputs  are 

deterministic  signals. 

the 

• 

M 


expected  values  in  Eqs.  (4.31)  and  (4.32)  are  replaced  by  the  signals  them¬ 


selves.  For  example,  ^  =  <Xl+(t)X2' (t)X3(t)>  = 


24 


XI  (t)X2 ' (t)X3 (t)dt 


=  0,  since  no  rise  transition  in  XI (t)  occurs  when  X2(t)  =  0  and  X3(t)  =  1. 

24  24 

•  4-  1 

Also  a 


*r i3  =  <Xl+(t)X2(t)X3'(t)>  =  ~ 


L 


Xl+(t)X2(t)X3'(t)dt  =  -jr  |  <5  (t-13)dt  = 

0  Jo 


— ,  because  XI (t)  has  a  rise  at  t  =  13  while  X2(13)  =  1  and  X3 (13)  =  0. 
Similarly,  all  other  weights  are  obtained,  and  the  values  are  tabulated  in 
Table  4.10. 

TABLE  4.10 


•  •  km 

12 

13 

22 

23 

32 

33 

arkm 

0 

1/24 

1/24 

1  /24 

0 

liv¬ 

* 

f  km 

2/24 

1/24 

0 

0 

1/24 

'd 

• 
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The  results  of  Eqs.(4.68)  and  (4.69)  are  used  in  Eq.(4.6)  to 
derive  f<Tr>(T)  and  f<T£>(T)>  respectively.  These  p.d.f.'s  are  used 
to  obtain  an  estimate  for  the  expected  value  of  the  output  2^(t) ,  denoted 
by  E^[2_(t)].  Computer  results  for  f<Tr>(x)»  and  E^[jl(t)]  are 

plotted  in  Fig.  4.13.  Features  of  E^[Z^(t)]  and  the  figures  of  merit  are 
tabulated  in  Table  4.11. 


TABLE  4.11 


Pulse  IZ 

AZ 

BZ 

cz 

DZ 

HI3 

0.45 

0.45 

1 .00 

1 .00 

WI3 

- 

- 

1.80 

5.90 

CI3 

3.40 

5.85 

11.15 

17.15 

hI3 

55.0 

0.28 

0.0 

0.0 

WI3 

- 

- 

0.125 

0.063 

CI3 

0.35 

0.45 

0.08 

0.15 

4)  The  Assigned  Gaussian  delays.  Note  that  the  minimum  and  maximum 
values  of  rise  and  fall  propagation  delays  given  by  Eqs.  (4.41)  and 
(4.42)  are  the  same  as  the  minimum  and  maximum  values  for  the  p.d.f.'s 
of  the  delays  shown  in  Fig.  4.6.  Therefore,  the  results  of  Eqs. (4.47)  and 
(4.48)  for  fTrG(T)  and  f  fG(x),  respectively,  are  readily  employed  as  the 
large  logic  block  output  delay  p.d.f.'s.  These  two  p.d.f.'s  and  the 
estimate  of  the  output  expected  value,  denoted  by  E^[Z(t)],  are  plotted 
in  Fig.  4.14.  Features  of  E^.[Z^(t)]  and  the  figures  of  merit  are  tabulated 
in  Table  4.12. 

TABLE  4.12 


Pulse  IZ 

AZ 

BZ 

CZ 

DZ 

HI4 

0.656 

0.656 

1.00 

1.00 

WI4 

0.6 

0.6 

2.15 

6.1 

CI4 

3.75 

6.25 

11.50 

17.50 

hI4 

81 

0.05 

0. 

0. 

WI4 

- 

0.5 

0.34 

0.03 

°14 

0.765 

0.82 

0.20 

0.50 
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Summary .  The  simple  output  delay  model  is  seen  to  have  serious 
short  comings.  Because  no  delays  are  assumed  within  the  ideal  logic 
circuit,  phenomena  like  the  nonoccurrence  of  the  intentional  pulse  A Z 
and  generation  of  the  unintentional  pulses  EZ  and  FZ  cannot  be  predicted 
by  the  model.  Such  phenomena  occur  when  durations  between  transitions 
in  the  various  input  signals  are  shorter  than  the  spreads  of  the  output 
delays . 

On  the  other  hand,  when  the  durations  of  pulses  and  gaps  in  the 
ideal  logic  output,  z(t),  are  sufficiently  long,  the  simple  output  delay 
model  does  a  reasonable  job  of  estimating  E[Z(t)].  None  of  the  four 
strategies  yielded  decidedly  superior  performance.  The  weighted  mean 
delay  strategy  did  slightly  better  in  predicting  a  smaller  probability 
of  occurrence  for  pulse  AZ.  However,  it  also  did  slightly  worse  since 
it  predicted  a  smaller  probability  of  occurrence  for  pulse  BZ .  The 
results  of  this  example  tend  to  suggest  that  all  four  strategies  perform 
approximately  the  same.  Since  the  assigned  Gaussian  delay  strategy  is 
the  only  one  which  does  not  require  convolutions,  and  therefore,  is 
the  simplest  to  implement,  it  appears  to  be  preferable. 
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5.  SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 


5.1.  Summary 

This  dissertation  develops  a  probabilistic  analysis  for  combin¬ 
ational  circuits  with  random  delays.  A  logic  block  is  modeled  as 
an  ideal  logic  circuit,  which  performs  the  intended  switching  function 
on  the  block  inputs,  followed  by  a  delay  element,  which  accounts  for 
propagation  delays  experienced  by  the  signals.  This  model  is  used 
both  for  single  gates  and  in  the  approximation  of  large  logic  blocks. 

Two  kinds  of  delay  elements  are  discussed.  The  pure  delay  element, 
whose  output  is  a  delayed  replica  of  its  input,  is  considered  first.  The 
more  complicated  case  of  the  discriminating  delay  element,  in  which 
input  rise  and  fall  transistions  experience  different  delays,  is  treated 
next.  The  output  of  the  discriminating  delay  element  is  a  distorted  ver¬ 
sion  of  its  input.  The  amount  and  complexity  of  the  computations  and  the 
required  computer  storage  are  significantly  larger  when  logic  blocks 
include  discrimating  delay  elements. 

Two  types  of  networks  are  considered.  The  simpler  type,  referred 
to  as  tree-like  networks  have  the  property  that  all  gate  inputs  within 
the  network  are  statistically  independent  provided  the  primary  inputs  to 
the  network  are  statistically  independent.  In  this  case,  evaluation  of 
output  expected  values  requires  knowledge  only  of  the  first  order  input 
statistics.  In  the  more  complicated  type  of  network,  containing  recon- 
vergent  fanouts,  higher  order  moments  of  the  primary  input  signals  may 
be  required  for  evaluating  output  expected  values. 

In  order  to  simplify  computations  and  reduce  the  amount  of  required 
computer  memory  an  approximate  model  for  a  large  logic  block  is  suggested. 


Different  placements  of  the  delay  elements  in  the  large  logic  block, 
model  are  considered.  The  simple  output  delay  model  is  developed  and 
selected  as  the  one  most  appropriate  for  EMI  applications.  Four  strate¬ 
gies  for  assigning  p.d.f.'s  to  the  output  delay  element  are  proposed. 

A  computer  example  is  used  to  compare  the  effectiveness  of  the  four 
strategies.  The  reference  waveform  for  the  true  expected  value  of  the 
output  signal  is  obtained  by  means  of  a  Monte  Carlo  simulation. 

5.2  Recommendations  for  future  work. 

The  probabilistic  analysis  developed  in  this  work  results  in  extreme¬ 
ly  complicated  expressions,  even  for  simple  circuits.  The  computational 
efforts  required  to  produce  numerical  results  are  considerable.  Therefore, 
if  the  methods  described  here  are  to  be  applied  to  practical  circuits, 
it  will  be  necessary  to  develop  an  efficient  computer  code  to  assist  with 
the  analysis.  Employing  techniques  such  as  used  by  Debany  [10]  in  obtain¬ 
ing  his  probability  expressions,  it  may  be  possible  to  automatically  gen¬ 
erate  complicated  expressions  such  as  appear  in  appendix  B.  This  would 
facilitate  application  of  the  probabilistic  analysis  to  larger  logic 
circuits.  Assuming  computational  difficulties  can  be  overcome,  other 
problems  remain  to  be  studied. 

Perhaps  the  most  natural  extension  of  this  work  is  to  sequential 
circuits  and  digital  circuits  containing  feedback  loops.  As  pointed  out 
by  Grundmann  [15],  signal  dependencies  caused  by  either  reconvergent 
fanouts  or  feedback  are  somewhat  similar  in  nature.  Therefore,  techniques 
developed  for  combinational  circuits  with  reconverg  nt  fanouts  may  provide 
a  good  starting  point  for  the  analysis  of  sequential  circuits  and  circuits 


with  feedback.  Some  preliminary  work  may  be  found  in  references  [15], 
[17],  and  [24]. 

An  interesting  theoretical  problem  deals  with  the  introduction  of 
blanking  effects.  Additional  knowledge  about  the  signals,  such  as  the 
statistical  properties  of  the  pulse  and  gap  durations,  6^  and  u . ,  are 
needed.  In  particular,  evaluation  of  p.d.f.’s  for  and  at  a  gate 
output,  given  the  gate  switching  function  and  the  p.d.f.'s  of  6,  and  u 

‘  — i  — ] 

for  gate  inputs,  is  required.  Introducing  the  blanking  effect  into 
the  model  eliminates  another  problem.  Recall  that  rise  and  fall  transi¬ 
tions  may  be  delayed  differently.  Let  a  rise  be  delayed  by  x  ,  and  a  fall 

rl 

by  T^,  where  .  Consider  a  short  pulse  of  duration  6^  such 

that  -  x^.  In  this  situation,  the  output  of  the  delay  element 

has  the  deficiency  that  the  fall  transition  associated  with  the  pulse 

occurs  before  the  rise  transition.  It  was  pointed  out  by  Grundmann  [15] 

that  this  effect  may  cause  significant  error  in  evaluation  of  the  output 

expected  value.  A  dual  effect  may  occur  when  x  2  <  and  a  short  gap 

of  duration  U  ^  ,  such  that  u  <  Tf2~Tr2*  intr°duced  to  the  input  of 

the  delay  element.  In  this  case  the  rise  of  the  gap  at  the  delay  element 

output  appears  before  its  fall.  These  two  problems  are  automatically 

eliminated  by  including  a  blanking  circuit  in  the  model  which  rejects 

any  pulse  shorter  than  x^  (which  is  certainly  larger  than  x  ^-x^  for 

any  x^)  and  any  gap  shorter  than  x  (which  is  clearly  larger  than 

T  .  — t  for  any  t  ) . 
f2  r  r 

Finally,  experimental  work  is  needed  to  obtain  p.d.f's  for  the 
delays  of  the  various  gates.  Measurement  of  delays  under  different 
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environmental  and  load  conditions  need  to  be  performed.  Analytical 
expressions  for  the  p.d.f.'s  can  be  derived  either  from  approximations 
to  well  known  distributions  as  pointed  out  in  [1,  Sec.  9],  or  by 
Gaussian  approximations,  as  suggested  in  [12]. 


APPENDIX  A.  EXPECTATION  AND  CORRELATION  OF  RANDOM  PROCESSES  INVOLVING  IMPULSES 


Throughout  this  work  expected  values  and  correlation  functions  of  0,1 
binary  random  processes  and  their  time  derivatives  are  evaluated.  It  was 
pointed  out  in  Ch.2  that  derivatives  of  sample  functions  in  a  0,1  binary 
process  consist  entirely  of  positive  and  negative  unit  area  impulses. 

The  mathematical  meaning  of  expectations  and  correlation  functions  of 
processes  whose  sample  functions  consist  entirely  of  impulse  trains  is 
discussed  in  this  appendix. 

The  unit  area  impulse,  referred  to  as  the  Dirac  delta  function  6(t), 
is  characterized  by  its  sampling  property  [27,p.35],  In  other  words, 
if  g(t)  is  any  function  which  is  continuous  on  a  neighborhood  about  the 
origin,  then 

oo 

J  g ( t ) 6 (t)dt  =  g(0) .  (A. 1) 

— oo 


By  a  simple  change  of  variable,  this  property  can  be  readily  modified  to 
yield 


g(t)  '  (t-tQ)dt 


g(<;+to)^(;)df.  =  g(tQ) 


(A. 2) 


provided  g(t)  is  continuous  in  the  vicinity  of  t=tQ.  There  is  no  ordinary 

function  that  has  the  sampling  property.  However,  the  "delta  function"  can 

be  viewed  as  a  limit  of  a  sequence  of  functions  { (t)  -  that,  as  n^0,  tend 

n 


't  -i 


to  zero  everywhere  except  t=0  where  they  tend  to  infinity,  and  the  area 
under  each  Sn(t)  is  unity  [25»Sec.l.5] .  As  an  example  of  such  a  sequence, 
consider 


6  (t)  •  —  e  2 

^  ./Ott 


(A. 3) 


which  is  a  Guassian  function  with  zero  mean  and  standard  deviation  equal  to  — 
This  function  has  unit  area  and  exhibits  the  desired  behavior  as  n-*». 

In  the  theory  of  generalized  functions  [?6]»  sequences  { yn ( t ) }  are 
defined  such  that,  for  every  finite  n,  the  function  yn(t)  is  infinitely 
differentiable  with  respect  to  t  and 


£im  t^y  (t)  =  0,  £im  t^  y  (t)  =0,  k>0. 

£-K_00 


(A. 4) 


Note  that  the  sequence  {^^(t) }  in  Eq.(A.3)  satisfies  these  conditions. 

A  generalized  function  G(t)  is  defined  by  an  equivalent  sequence (yn(t) } 
such  that 


K 


(t)g(t)dt  =  G(t)g(t)dt. 


(A. 5) 


In  particular,  the  Dirac  delta  function  is  defined  by  any  equivalent  sequence 


having  the  sampling  property.  Specifically, 

00  oo 

[  <5 (t)g(t)dt  A  £im  I  <5  (t)g(t)dt  -  g(0). 
J  n-*»  J  n 


(A. 6) 


Two  additional  important  properties  of  the  Dirac  delta  function  are 


6(-t)  -  5(t) 


(A.  7) 


and 


i 


<5(t-t0)g(t)  -  5(t-t0)g(tQ)  (A. 8) 

for  every  g(t)  which  is  continous  in  the  viclntiy  of  t*t^. 

The  definition  of  Eq.(A.6)  can  be  generalized  to  two  or  more  variables 
[26, .Sec. 7. 6]  by 


S(0)6(t)g(0,t)d8dt  ■  Him  Him 

nr*®  n -*00 


j  |  6m(0)6n(t)g<8’t)d0dt  -  g(o,o) 


(A. 9) 


for  any  function  of  two  variables  g(0 ,t)  which  is  continuous  in  the  vicinity 
of  the  origin  (0,0). 

The  notation  S(t)  is  commonly  used  as  a  function  symbol.  However, 
it  has  a  meaning  only  as  a  factor  of  an  integrand.  Integrated  over  the 
variable  t.  Similarly,  the  "product"  5(8)6(t)  is  meaningful  only  as  a  factor 
of  an  integrand  under  a  double  integration  over  the  variables  0  and  t. 

Attention  is  now  focused  on  random  variables  and  stochastic  processes. 


Let  i  be  a  random  variable  with  p.d.f.  denoted  by  f  (x),  and  let  g(x)  be 

f°° 

a  function  of  x  such  that  the  integral  g(x)fT(T)dx  exists.  The  expected 

— CO 

value  of  g(x)  is  given  by  [18,p.l42J 

E[g(x)]  =.  J°°g(x)f.L(x)dx  (A,  10) 


(T)dx 


(A.  10) 


If  i  is  a  random  variable  of  a  discrete  type,  with  values  taken  from 
the  set  {x^}  ,  then 


E  [  g (x ) ]  -  l  g(x.)P  {  x  -  x  } 
i 


(A. 11) 


......  .*• 

*-'  •.*  *.* 


Suppose  X(t)  is  a  stochastic  process  such  that  XI  (t)  ■»  X(t,x). 

In  this  case,  the  source  of  randomness  is  the  random  variable  x  and 

for  the  specific  outcome  x_  =  x,  the  sample  function  X(t,x)  is  deterministic. 


Example  A.l  Let  5C(t)  =  A  sin((jdt+x) .  Specifying  a  known  value  for 
the  random  variable  x_,  i.e.,  setting  x_  =  x,  results  in  the  deterministic 
sample  function  X(t)  =  A  sin(ait+x) .  V 

As  in  Eq.(A.lO),  the  expected  value  of  X(t)  =  X(t,x_)  is  given  by 


OO 

E[X(t) ]  =  j  X(t,x)fT(x)dx. 

—  OO 


(A. 12) 


Note  that  EfX(t)]  is  a  deterministic  function  of  time.  If  x_  is  a  discrete 
random  variable,  it  follows  from  Eq.(A.il)  that 


E[X(t)]=  l  X(t,xi)Pr  {x  =  x±}  . 
i 


(A. 13) 


Let  X^Ct)  be  a  stochastic  process  whose  sample  functions  consist  of 
a  single  impulse  occurring  at  the  random  time  instant  X_.  Specifically, 


X^t)  =  <5(t-x)  . 


(A. 14) 


The  p.d.f  of  the  random  time  instant  x_  is  given  by  f^X).  Direct  sub¬ 
stitution  into  Eq.  (A. 12)  results  in 


EtXjU)]- 


6(t-x)f^(x)dx 


(A. 15) 


Changing  the  variable  of  integration  yields 


E[Xd(t)] 


00 

* 

5(X)fx(t-T)dT. 

00 


(A. 16) 
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Utilising  the  sampling  property  of  Eq.(A.l),  one  readily  obtains 


ElX^t)]  -  fx(t). 


(A. 17) 


However,  the  meaning  of  the  expectation  of  X^(t)  Is  better  understood 
using  the  sequence  approach,  as  in  Eq.(A.6).  Application  of  Eq.(A.6)  to 
Eq.(A. 16)  results  in 


Et.XjCt)] 


js(T)fx(t-T)dT 

—CO 


£im 

n-**> 


|  6n(x)fT(t-T)dT. 

-CO 


(A. 18) 


Note,  for  any  finite  n,  that  the  integrand  6n(x)fx(t-x)  is  a 
finite  function  of  T  for  every  t  and  the  meaning  of  the  integral  is  clear. 
Letting  n-*00  yields  the  result  of  Eq.(A.17).  If  T  is  a  random  variable 
of  a  discrete  type  it  follows  from  Eq.(A.ll)  that 

EfX^t)]  -  l  6(t-x1)Pr{r  -  T±}.  (A. 19) 

Even  though  X  (t)  consists  of  a  single  impulse,  the  expectation  obtained 
in  Eq.(A.19)  is  an  impulse  train  with  Impulses  occurring  at  the  time  instants 
specified  by  the  set  {x^}  and  having  areas  equal  to  the  corresponding  proba¬ 
bilities  of  occurrence.  In  this  work  all  expectations  evaluated  for  processes 
whose  sample  functions  contain  impulses  are  eventually  integrated.  Con¬ 
sequently,  ordinary  functions  are  always  obtained  from  expressions  such  as 
that  in  Eq.(A.19)> 

The  result  in  Eq.(A.17)  is  readily  generalized  to  processes  having  a 
finite  or  countable  number  of  impulses  in  their  sample  functions.  Let  the 
process  X*!(t)  have  sample  functions  consisting  of  N  impulses  which  occur 
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at  the  random  time  instants  t^, . » . ,t^. 

random  time  Instants  is  given  by  f  . 

Cl  •  • 

is  expressed  as 


The  joint  p.d.f.  of  the  N 
tN  The  process 


N 


l  5(t-0. 
k-1 


(A.  20) 


The  expected  value  of  X*!(t)  is  given  by 


oo  oo 


Et^Ct;)]  -  [...  £  6(t-tk)ftl  ^tN(t1,...,tN)dt1...dt1}.  (A. 21) 

j  '  k=*l  —  — 


— oo  — oo 


Changing  the  order  of  summation  and  integration  and  using  marginal 
p.d.f .s  determined  by 

ftl. • .tN(tl . tH)dtr  * "dtk-ldtk+l *  * ,dtN’ 

—OO  —CO 

(A. 22) 


ftk(tk> 


i 


result  in 


E[^(t)] 


N 


(A. 23) 


Note  that  Eq.(A.23)  is  a  direct  generalization  of  Eq.  (A. 17). 

The  autocorrelation  of  a  stochastic  process  X(t)  is  defined 
as 

Rj^Ctj.tj)  =  E[X(t1)X(t2)]  .  (A. 24) 


Once  again,  let  X^(t)  ■  5(t-_x)  be  the  stochastic  process  whose  sample 
functions  consist  of  a  single  unit  area  impulse  at  a  random  time  instant  t 


By  definition,  the  autocorrelation  of  X^Ct)  is  given  by 

00 

Rxd  Xd(tl,t2)  "  "  J6(t1-T)6(t2-T)fT(T)dT.  (A. 25) 

00 


Denote  6 (t2-x)f^(x)  by  g(x).  The  symmetry  property  of  Eq.(A.7)  yields 
S(tj-x)  -6(x-tp.  Application  of  the  sampling  property  in  Eq.(A.l)  results  in 


*Scd  Xd(tl*t2) 


00 

* 

<$(x-t.)g(x)dx  -  g(t.) 
1  1 

00 


<S (t2“ti)  fT(tP* 


(A. 26) 


Note  that  R^  xd  (t^)  is  an  at  t2*tl  w:*-t^  an  area  equal  to  f^Ct^. 

Note  that  application  of  the  property  in  Eq.(A.8)  yields 


<S<t2-ti)fx (ti>  *  <S(t2~t1)fT(t2). 


(A. 27) 


Hence,  as  expected,  Xd  is  symmetric  with  respect  to  t^  and 

t2 .  As  with  expectations,  all  correlation  functions  evaluated  for  pro¬ 


cesses  whose  sample  functions  contain  Impulses  are  eventually  integrated. 
Once  again,  ordinary  functions  are  obtained  (see  Ex. 3. 12). 


APPENDIX  B.  EXACT  EXPRESSION  OF  E[z (t) J 
FOR  CIRCUIT  USED  IN  THE  COMPUTER  SIMULATION 


To  demonstrate  the  complexity  which  may  occur  in  the  analytical 
expression  for  the  output  expected  value  of  even  a  relatively  simple 
circuit*  the  exact  expression  of  E[£  (t)]  for  the  circuit  used  in  the 
computer  simulation  of  Sec.  4.3  is  given  below.  This  is  expressed 
in  terms  of  the  circuit  inputs,  XI (t) ,  X2 (t) ,  and  X3(t) ,  their  counting 
signals  and  their  derivatives.  Statistical  independence  is  assumed 
among  all  three  inputs. 


E[z+(t)]  =  E[Xl+(t-  irA1)  3E[X2(0)]E[3C3(0)  ]  +  ElXl+U-T^)  ] 


rt_— r01~— rAl 

J  W  (0.e)Rx_3x3+(0’e)de+  <t-IrAI)1 
0 


t-T  n1-T 

— rOl  — rAl 


W(O*0)BX3X3'  (°*9)d9'  ^-IrAl” 


t  -fOr^rAI 


RX2X2' (O,0)RX3X3_(O,0^d0~  E[—  ^  ^rAl^ 


t-T  -T 

L  -fOl  -rAl 


■St2b-  <°-9>RX3X3'  (O'e)d0  + 


RX2X2' ^0,t--r0l“-rApRX3X3+  ^  0,t“-rOl”-rAl^ 


+  E[L1(t“IrA1)]Rx2X2+(0,t“^rOr-rAl)RX3X3,(0,t--rOr-rAl) 


t  -^rk2  -^rOl  -h-Al  ,  -^£01  -rAl 
+  J{E[Xl+(t-TrA1)Rx2x2(O,0)E[X3.  (0)] 

0  0  0  ~~~ 


+E[Xl+(t-_lrA1)]  Rx2*2+(O,0)E[X3(0)]  +  1^2X2' (0*O)RX3X3+  (0’9) 

+  E[X1  ^“lrAl^Rx2'x2+^0*6^SX3X3^0*e^  +  ^H('t~lrAl^RX2X2  1  >n^RX3+X3'H(9 ,n^ 

+  E[xi+(t-irA1)l^2^2+  (n.0)Rx3x|(0.n) 

+  E[X1^  ^ ~Ir A1 )  J Rx2xf* 9 * n ^  RX3'X3+^n*9^  +  ^fc~-^rAl^  ^RX2+X2+^0>ri^RX3X3 '  ^9,r^ 
-  E[X1  (t~lrA1)  ^Rx2X2'  ^9>^RX3^X3T^9>^  “  E[X1  (t~XrA1)  1Rx2,X2+^,9^RX3X3~^9*^ 
-  E[X1.  (t-lrA1)]Rx2_x£(e»C)R^3,^3-(C»9)-E[Xi+(t-JTrA1)]R^2+^2-(0>C)Rx3x3,(0>C) 

+E [ XI ( t-lrA1 ) 1 rx2x2 1 ( 0  * t_— ro l _-rAl ) RX3+X3+ ( 9 ’ t_-rO 1 "-rAi ) 

+E [XL(t-TrA1 )  ] ®X2 ' 'X2+  ( t-Ir0rIrAl  * 9)RX3X3+(9,t~^rOr^rAl) 


+E [ XI  ( t  T_r A1 )  ] ^2X2^ ( 6 ’ C  ~-r  0 1  “-rAl ) *Sc3 ' X3+  ( C  _-r 0 l_Ir A 1 ’ 9) 


+E  IrAl )  ]  RX2+X2+  ( 0  ’ c  -—rO  1  ~^rk  1  ^  RX3X3 ' ( 9 ’ C  ~ir 0 1  -^r A1  >  > d  ^dr)d  6 


ft_±fA2  ft~^rOrIrAl  ft_^fOr^rAl 


{E[Xl+(t-TrA1)]Rx2x2(O,0)EtX3"(9)] 


+E[Xl+(t-TTA1)]Rx2J2-(O,e)E[X3(0)]+E[)a+(t-TrA1)]Rx2x2,(8,O)Rx3—  (0,0) 


+E[xi+(e-TrAI)Rx2,i2-<o.e>Rx3x3(0.e)  +  E(xi(t-itA1)]R12H,(e,n)i^3+j3- 


-(n,0) 


«lEi+(t-Tt.A1)IR)t2,j2-(n.e)Rx3j3+(e,n)+EtxiT(t-Tr41)]Rx2-2+(e,n)i<23,i3-<ii,0) 


«  4- 

nTtrl  1  / 


+E[Xl+Ct-lrA1)lR^2+^2-(n»9)Rx3x3.  (9,n)-E[ja  (t~IrA1^RX2x2,^9,C^RX3  X3 


-E[X1  (t-xrA1)]Rx2,^2-(^,9)Rx3jj:3-(6.0-E[ja  (9,C)RX3'X3  (C* 


•  4- 

p,  r,n  1  t 


+E  [2LL<t-JLrA1 )  ]  ^2X2’  (  8 » t--r0 1  ~^r A1  ^  RX3+X3 ' "  ^  t_-r0 1  "^rAl , 9  * 


+E  [  XI  ( t-T_r  A1 )  3  ^2 '  x2-  ( t--r0 1  r  Al’ 9 )  RX3X3+  ( 9  * C  _-r  0 1  _-r  A 1 J 


+E[^(t-TrA1)]RX2x2+(0»t-TrO1-TrA1)Rx3,^3-(t-Troi-lrAi»6) 


4*. 


+  B[Xl(t--TrA1iHjfo-fjk+  ^t-2r01~JTrA1 . - )rX3X3 ,  (Q.t-TTQ1~lrA1))  dgdnde 


+  E[X1<0)  11^(1,  t-IIA2>*IB+<'-irA2>  1+Elil«»  '8X2x2+<°'t-IrA2)E[S<'-IrA2) 1 


*--4ai  e-4oi  9-ifoi 

+  |  I  1  {*‘*‘*<e»*m2<°-c-lrA2>*&*<t-lrA2> 

0  0  0 


+  E[Xl+(0)]Rx2^(O.t-TrA2)EtX3(t-lrA2)] 

+E[X1  (9)  1^2x2'  (t~^A2,0)RX3X3+(0,t~ITA2)'t'E[—  (6)  ]RX2  ,X2+(0,t:~-^TA2)RX3X3(0,t 
+  E[3Q  (6)  1^2x2'  (t~lrA2*n)RX3+X3+^t~  •^tA2’T1) 

+  E[)Q  (0>JRx^»x2+(O»t-TrA2)R^3^+(t::TrA2,n) 

+  E[X1  (9)Rx2x2+(t~^rA2>n)RX3,X3+(n,t:~^rA2^ 

+  E[X1  <9)]R^+^+(t-IrA2*r>)Rx3x3«(t-IrA20) 

-  E[XI  <6)RX2X2I (t~lrA2,^)RX3+X3~(t~^rA2,n) 


-  E[X1  (9)lRX2'x2+(?,t“irA2)RX3X3“(t"^A2,^^ 


-  E[X1  (9)iRX2X2~^t~^rA2»^RX3,X3't'^,t~-n:A2^ 

-  E[X1  (6)]Rx2+x2"<t-i1.A2»^>Rx3X3'  (t~IrA2*^) 
+E[X1(9)]R^2x2,  (t-X,.A2>6-  Iroi)Ri3+x3+(t'lrA2*0“-^rOl) 


+  E[Xl(0)RX2,^2+(9-  irOi»t‘IrA2)RX3X3+(t~-rA2,0_ 


+  E[X1  (0)lRX2X2+(t_-rA2,e“irOl)RX3'X3+(e_  ^r0rt~^rA2) 


+  E[2a(0)R^2+*2+(t  ITA2‘e-IToi)Rx3X3,^t-'^rA2,0--irOl^d^dnde 


C  -fA2  ,-0  -rA2  r  6  -fOl 


O  1  -3  {E[  —  (0)-|RX2X2(O,t  ^rA2)E[—  (t  -^rA^1 


+  E[3U  (0)  ^Rx2X2+^°’t  — rA2^E^— rA2^+E^—  ^0^RX2X2'^t  — rA2,0^RX3X^0,t  — rA2^ 


+  E[Ll~(8)]Bx2.fa+t°»t--J,A2>BX3X3<0*t0W 


Ef—  ^6^RX2X2'  ^  ■^2,n^RX3+X3+^t_-rA2,T1^ 


+  E[_X1  (°>]Rx2,X2+(ri,t~-rA2)RX3X3+(t~^rA2,rl) 
+  E[X1  (9)]Rx2X2+(t~^rA2,r|)RX3’X3'l'(n,t~-^rA2^ 
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+  E[X1  (9)  ]Rx2+X2+^t~-^rA2,T1^RX3X3'  ^t~-rA2>1'1^ 

-  E [XI  (9)  ]Rx2X2*^t~-^rA2 ,^^RX3+X3~^t~— 'fAZ 

-  E [Xl_  (B)]Rx2,^2+^»t-lrA2)RX3X3~^t~^TA2,^) 

-  E[n-(0)]Rx2b-(t-TrA2,ORx3,i3+(?,t--TrA2) 


E[X1  (9)Rx2~t~X2~^t~— rA2,^^RX3X3'  ^  -rA2’^ 


+  E[ja(6)]Rx2x2(t-T_rA2,9-  lfoi)RX3+3G_(t_-irA2,e''^f01) 

+  E[n(e)]Rx2,b+(9-If01,t-TrA2)Rx3^3-(t-TrA2>e-Tf01) 


+  E[jq(9)]Rx2j|:2-(t-TrA2,9  ifoi)Rx3'x3+(e-  ^£01  ,t_^rA2) 


+  E[2a(0)]R^+^-(t-TrA2,e-TfO1)FX3M,(t-TrA2,0-Ifol)}  dUndB. 


The  above  expression  contains  75  terms.  As  pointed  out  in  Sec.  3.3 
the  various  correlation  functions  can  be  determined  from  a  much  smaller 


set  of  basic  correlation  functions. 


REFERENCES 


[1]  A.  Ephrath  -  "EMC  modeling  and  analysis  -  a  probabilistic  approach"  - 

M. S.  Thesis  -  Syracuse  University,  1982  (Adviser  -  D.  D.  Weiner). 

[2]  J.  Gormadi,  D.  D.  Weiner,  et  al  -  "  Random  Susceptibility  of  an 

IC  7400  TTL  NAND  gate"  -  1983  IEEE  Symposium  on  EMC,  Washington,  D.C., 
Aug/23-25/1983. 

[3]  J.  Alkalay,  D.  D.  Weiner  -  Performance  degradation  of  a  7400  TTL 
NAND  gate  due  to  sinusoidal  interference  -  RADC-TR-80-257 , 

Aug.  1980. 

[4]  L.  Fratta,  U.  G.  Montanari  -  "A  Boolean  algebra  method  for  com¬ 
puting  the  terminal  reliability  in  a  communication  network"  - 
IEEE  Trans,  on  Circuits  and  Systems,  May  1973. 

[5]  K.  P.  Parker,  E.  J.  McCluskey-"  Probabilistic  treatment  of  general 
combinational  networks"  -  IEEE  Trans,  on  Computers,  June  1975. 

[6]  K.  P.  Parker,  E.  J.  McCluskey  -  "Analysis  of  logic  circuits  with 
faults  using  input  signal  probabilities"  -  IEEE  Trans,  on  Computers, 
Mav  1975. 

[7]  S.  K.  Kumar,  M.  A.  Breuer  -  "Probabilistic  aspects  of  Boolean 
switching  functions  via  new  transforms"  -  Jour,  of  the  Assoc,  for 
Computing  Machinery,  July  1981. 

[8]  J.  D.  Mur ch land  -  "Fundamental  concepts  and  relations  for  reliability 
analysis  of  multi-state  systems"  -  SIAM,  Reliability  and  fault  tree 
analysis,  Philadelphia  1975.  Editors:  R.E.  Barlow,  J.  B.  Fussell, 

N. B.  Singapurwalla. 

[9]  W.  G.  Schneeweiss  -  "Calculating  the  probability  of  Boolean  expression 
being  1"  -  IEEE  Trans,  on  Reliability,  April  1977. 

[10]  W.  H.  Debany  -  "Probability  expression  with  applications  to  fault 
testing  in  digital  networks"  -  M.S.  Thesis  -  Syracuse  University, 

1983  (Advisor  -  P.K.  Varshney) . 

[11]  G.  Coraluppi  -  "Binary  network  analysis"  -  Alta  Frequenza,  Feb. 1963. 

[12]  B.  Magnhagen,  R.  Flishberg  -  "A  high  performance  logic  simulator  for 
design  verifications"  -  Summer  Computer  Simulation  Conference, 
Washington  D.C.,  July  1976. 


[13]  B.  Magnhagen  -  "Practical  experiences  from  signal  probability 
simulation  of  digital  design"  -  Proc.  of  the  14th  Design  Auto¬ 
mation  Conference,  New  Orleans  LA,  June  1977. 

[14]  H.  K.  Al-Hussein,  R.  W.  Dutton  -  "Path  delay  computations  for 
integrated  systems"  -  1982  IEEE  Circuits  and  Computer  Conference, 

New  York  City,  NY,  Sept.  1982. 

[15]  J.  W.  Grundmann  -  "Probabilistic  analysis  of  digital  networks"  - 
Ph.D.  dissertation  -  Purdue  University,  1979  (Adviser  -  S.C.  Bass). 

[16]  S.  C.  Bass,  J.  W.  Grundmann  -  "Expected  value  analysis  of  combin¬ 
ational  logic  networks"  -  IEEE  Trans,  on  Circuits  and  Systems,  May 
1981. 

[17]  J.  W.  Grundmann,  S.  C.  Bass  -  "Expected  value  analysis  of  digital 
networks  with  memory"  -  IEEE  Trans,  on  Circuits  and  Systems, 

Sept.  1981. 

[18]  A.  Papoulis  -  Probability,  random  variables  and  stochastic  processes 
McGraw  Hill,  1965. 

[19]  Z.  Kohavi  -  Switching  and  finite  automata  theory  ,  2nd  Ed.  -  McGraw 
Hill,  1978. 

[20]  D.  D.  Givone  -  Introduction  to  switching  circuits  theory  - 
McGraw  Hill,  1970. 

[21]  N.  Balabanian  -  Digital  logic  and  sequential  machine  design  - 
Class  notes,  Syracuse  University,  1983. 

[22]  E.  Kreyszig  -  Advanced  engineering  mathematics  ,  5th  Ed.  -  John 
Wiley,  1983. 

[23]  "Integrated  circuit  electromagnetic  susceptibility  handbook" 

MDC  report  E  1929,  sponsored  by  USN  Surface  Weapon  Center, 

Contract  N60921-76-C-A030. 

[24]  K.  P.  Parker,  E.  J.  McCluskey  -  "Sequential  circuit  output  proba¬ 
bilities  from  regular  expressions."  IEEE  Irans.  on  Computers, 

March  1978. 

[25]  M.J.  Lighthill  -  Introduction  to  Fourier  analysis  and  generalised 
functions  Cambridge  University  Press,  1959. 


[26]  D.S.  Jones  -  Generalised  functions  -  McGraw  Hill,  1966. 

[27]  R.  F.  Hoskins  -  Generalised  functions  -  Ellis  Horwood(John  Wiley) 


V  < 

^  * 

\  MISSION  < 

$  Of  l 

^  Rome  Air  Development  Center 

RAVC  plans  and  executes  neseanch,  de.vzlopme.nt,  teAt  and 
-6 elected  acquisition  pn.ogn.amA  in  Support  o£  Command,  Contnol 
&  Communication!,  and  Intelligence  (C^ j )  activities.  Technical 
1%  and  engineening  Auppont  within  areas  technical  competence 
is  pnovided  to  EST>  Pnognam  O^ices  l  POa  )  and  othen  ESD 
ij*)  elements.  The  principal  technical  mission  aneas  one 
?  communications,  electromagnetic  guidance  and  contnol,  sur- 
^  veillance  o $  gnound  and  aerospace  objects,  intelligence  data 
\  collection  and  handling,  information  system  technology,  ' 

ionospheric  pnopagation,  solid  state  sciences,  microwave 
%  physics  and  electronic  neliability,  maintainability  and 
&  compatibility. 


END 

FILMED 

3-85 


DTIC 


